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Chapter 1 

Introduction 



The study of the parity violating electron scattering (PVES) has played a 
fundamental role in corroborating the Standard Model (SM) of electroweak 
(EW) interactions and understanding the detailed structure of the weak- 
neutral current. The first experiment of this kind was performed at the 
Stanford Linear Accelerator Center (SLAC) in the late 1970's [T] where the 
parity violating asymmetry 

AnL^'-^, (1.1) 

where a^i and are the cross sections for a right handed and left handed 
polarized beam respectively, was measured in the scattering of longitudinally 
polarized electrons off a deuteron target. Since the electromagnetic parity 
conserving (PC) interaction is insensitive to the polarization of the beam, 
while the electroweak parity violating (PV) interaction changes sign upon 
flipping the spin of the incoming beam, taking the difference between gr 
and ai isolates the PV contribution to the cross section, which ultimately 
depends on the weak mixing angle. In the SLAC experiment, the value of 
the weak mixing angle extracted from Am indicated an amount of parity 
violation which was consistent with the SM. Beside establishing the SM as 
one of the most promising frameworks to describe the electroweak interac- 
tion, this experiment also ruled out a number of possible models which were 
considered at the time as plausible. Subsequently, PV experiments involv- 
ing scattering of electrons off ^^C [3] and ^Be [1] were performed to further 
test the SM at low energy. In the last few years, because of the constant 
improvement of experimental techniques, measurements of Am at the level 
of a few percent, are conceivable. Such high precision allows one to measure 
radiative corrections to the SM and look for deviations from the SM, open- 
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ing a window on possible physics beyond the SM, and also to have better 
insights into the hadronic structure. A quantity that is being extensively 
studied to find deviations from the SM, is the weak mixing angle 9wi defined 
by the relation 

'2 

siii^ 0w = -^n-, — 2' (1-2) 
g +9 

and where g' and g are SUl{2) x ?7y(l) guage couplings. An important fea- 
ture of sin^ Owi is that it value "runs" (i.e. changes) with energy. They way 
it changes with energy is described by the renormalization group equation, 
which is a first order differential equation, therefore, by measuring the angle 
at a given energy, one could predict its value at any energy. The standard 
model (SM) of electro-weak interactions makes a precise prediction on the 
running of the electro-weak angle, and measuring it in high precision exper- 
iments in different kinematical regions, would be a powerful way to test the 
SM and to look for new physics. 

In this respect, a series of high precision low energy PV experiments have 
been recently performed, approved or proposed. Among them, it is worth 
to mention the M0ller experiment E-158 performed at SLAC, where the 
weak mixing angle is extracted by measuring Am in electron-electron scat- 
tering with an accuracy of (5(sin^ 9^) ~ .0007, and the Qweak experiment [7] 
approved at the Jefferson Laboratory (JLAB) facility, which will measure 
the weak mixing angle by measuring the on a proton target, aiming 
for a precision of (^(sin^ 9y/) ~ 0.0005. In both the experiments the average 

is 0.03 GeV^. As pointed out in Ref. [5], the sensitivity of an experi- 
ment to new physics depends on the kinematical region and on the target 
that is used. The experiments in Refs. [6l[7], in tandem with the proposed 
experiments of Refs. [HI [9], will form a complementary program which will 
explore a broad range of possibilities to extend the SM. In Refs. [HI [9] they 
propose a measurement of Am on a deuteron target in the deep inelastic 
region, with a momentum transfer squared of the order of 3 GeV^ in Ref. [9], 
and between 16 and 28 GeV^ in Ref. [8], and extract the weak mixing angle 
from the asymmetry, with a precision of 0.5% on sin^ 9y/- In particular, the 
latter experiment would explore a kinematical region far from the Z pole, 
but at a sufficient high so that, presumably, one should minimize the 
uncertainty due to the higher twist (HT) effects and also be able to test new 
physics. Following Ref. [5], one could describe the effects of new physics, by 
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introducing a four fermion contact interaction, whose Lagrangian is 



where A is the characteristic mass that sets the new physics energy scale, k'^ 
is the couphng of the new interaction, and hij are the hehcity-dependent cou- 
phngs (i , J denote the hehcity of the fermions). Such a form of interaction is 
suitable to describe different scenarios of new physics which could be tested 
by measuring the deep inelastic scattering (DIS) Arl on deuteron. For in- 
stance, super string inspired theories and some super symmetric theories 
admit the existence of an extra gauge boson, and the proposed experiment 
might set a lower bound on its mass. One could also test compositeness of 
fermions, which could be described by a contact interaction similar to the 
one in Eq. ()1.3p . and set limits on the mass scale A at which the composite- 
ness becomes manifest. Finally, Eq. (jl.Sp could also be applied to describe 
leptoquarks formation, i.e. the formation of an intermediate bound state of 
a lepton and a quark, in the process e + (7 — > LQ — > e + g where LQ is the 
leptoquark state. In this case the parameter A represents the mass of the 
leptoquark state, and the proposed experiment could be used to set a lower 
limit on it, given some reasonable assumptions on the scale of the couplings. 
A PV deep inelastic scattering (DIS) experiment might be proposed in the 
near future, whose goal is to measure HT effects [10]. In this experiment, 
which might eventually run at JLAB, a PV asymmetry is measured using 
a deuteron target at a of nearly 3 GeV^. My study on HT effect might 
provide a guidance in developing the details of the experiment, in particular 
at what kinematics and for which values of the Bjorken variable xb, the HT 
effects could be prominent. 

Because of the high precision of these experiments, it is absolutely nec- 
essary to know how the hadronic effects might impact the interpretation 
of AfiL. In this study, which was mostly conducted along with Michael 
Ramsey-Musolf and Shi-Lin Zhu, I consider a few of such effects which 
might be relevant to the interpretation of the different experiments. 




(1.3) 
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1.1 Complete Treatment of the Parity Violating 
Asymmetry at Leading Twist 



In the second Chapter, I consider an extensive study of the possible leading 
twist hadronic corrections which affect Afij^ in a deep inelastic scattering 
(DIS) process. The main motivation for studying such corrections is that in 
Refs. [HI [9], it has been proposed to use the deep inelastic polarized cross 
section to extract sin^ 6w with an accuracy of 0.5%. The corresponding 
precision required to measure Arl is therefore of the order of 1%, which is 
extremely challenging for this kind of experiment. Consequently, one has to 
make sure that all the possible theoretical uncertainties are below 1%. In the 
approximate case in which sea quarks are neglected, the mass of the target 
M is negligible compared to the momentum transfer Q, and the perturbative 
QCD (Quantum Chromo Dynamics) corrections are not included, Ajil for 
an isoscalar target like the deuteron assumes a particularly simple form [11] 



Gpq^ 9 

Arl = 



(l_^sin2^^) + (l_4sin2 0W')^^ 



(1.4) 

where Gp = 1.16639(1) x 10~^ Gev~^ is the Fermi constant, a = is the 
fine-structure constant, q is the momentum transfer and y = 1 — ^ with 
E [E') the energy of the incoming (outgoing) electron. It's worth noticing 
that, under the previous assumptions, the asymmetry does not depend on 
the structure of the target, but a more complete treatment, which includes 
the hadronic corrections, is necessary if one wants to reach a precision of 1% 
in the interpretation of the measurement of Arl ■ 



1.1.1 Sea Quark, Perturbative QCD and Target Mass Con- 
tributions to Arl 

In Chapter [2] I present a detailed computation of Arl for a deuteron target, 
including the sea quark (see for instance Ref. [2]), the perturbative QCD 
(PQCD) [12] and the target mass (TM) [131 E] corrections. Eq. has 
been obtained by using the quark parton model, where the cross section of 
the process of the lepton scattering off the target (deuteron, in this case), 
is viewed as an incoherent sum of cross sections of the lepton scattering 
off free pointlike particles (the partons). If one restricts the process to 
only valence quarks, Eq. ()1.4p is recovered. The sea quark corrections, 
consist in including all the possible quarks which could participate to the 
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Vq^Gq YG^qq YG^qq 

Figure 1.1: Order as contribution to eP — > e' X DIS process (only the 
hard part is shown) arising from antiquarks and gluons 

process, compatibly with the kinematical constrains. The PQCD corrections 
are corrections of first order in the strong coupling constant as and which 
cannot be reabsorbed in a redefinition of the parton distribution functions 
(PDFs). Such terms arise from corrections to the processes depicted in Fig. 
11.11 which go beyond the leading logarithmic approximation, and they are 
renormalization scheme dependent (see Sec. 12. 2p . Since these correction are 
expressed in terms of integrals of the leading twist PDFs, the are known once 
the latter are. In the PQCD correction, also the gluon PDF contributes. 
Finally the TM corrections are contributions which are suppressed by a 
factor of M^/Q^, where M is the target mass and Q the momentum transfer. 
When deriving the naive parton model, one assumes that the nucleon is made 
up of free partons that carry a faction x (ajxjl) of the total its momentum 
P^. If one considers also the possibility for the partons to interact among 
each other through the exchange of gluons, their momentum acquires also 
a component which is transverse the nucleon momentum. Because of this 
transverse component of the parton momentum, the PDF obtained in the 
naive parton model receive corrections which are proportional to M^/Q^ \12>\ 
I47j . In DIS processes, where the momentum transfer Q is normally much 
larger that the nucleon mass M, those corrections are small, but they might 
become important at the kinematics proposed in [8l |10], where = 2.8 
GeV^. Also in this case they can be computed if the PDFs are known, but 
the gluon PDF does not contribute. 

In Sec. 12.11 1 give an overview of the DIS experiments to better un- 
derstand the role they played in the last few decades in understanding the 
hadronic structure and how they lead to the formulation of the quantum 
chromo dynamics (QCD) gauge theory. In Sec. 12.1.11 1 introduce most of 
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the nomenclature involved in DIS, from the kinematical variables, to the 
hadronic and leptonic tensor, in which terms one can express the cross sec- 
tion. Finally, in Sec. 12.21 I consider Arl in the case of deuteron target, 
including the corrections which were neglected in Eq. ()1.4p . As mentioned 
before, all the corrections considered above can be expressed in terms of 
the PDFs. These quantities, at the moment, cannot be computed from first 
principle (although there are some attempts to obtain them from lattice cal- 
culation [E]), and one has to rely on model independent extractions. Two 
groups, in particular, the CTEQ [16] and the MRST [IT], have, in the last 
years, continuously updated their PDFs extracted from the experiments. By 
using the most recent sets of PDFs provided by the two groups, I compute 
the leading twist hadronic corrections to Aju, and take the difference be- 
tween the two predictions as a theoretical uncertainty arising form the PDFs. 
Although the two sets quite agree in describing the valence distribution, the 
sea quark and the gluon distributions might substantially differ. These dif- 
ferences might represent a source of uncertainty in the interpretation of the 
measurement of Arl, especially because of the gluon distribution, which is 
the least known, and appears in the PQCD corrections. 



1.2 Higher Twist Corrections to Aji^ 

Chapters [3] and U] are dedicated to the definition and estimate of the twist 
four correlation functions. Similarly to the TM corrections, the twist four 
corrections appear in the DIS cross section and asymmetry suppressed by 
a factor of ^ compared to the twist two contributions, where A is the 
QCD scale of the order of few hundreds MeV. Because of the small value 
of the momentum transfer which will be eventually used in the experiments 
proposed in [51 [10] , these corrections might give a substantial contribution to 
Arl ■ The main difference between TM and HT corrections is that the former 
can be computed in terms of PDFs, which are known from experiments 
. As far as the HT contributions are concerned, there exists only some 
model independent extractions from experimental data [181 [T9] which are 
not extremely precise. The problem is that it is very hard to disentangle 
and distinguish between the contributions arising from the HT and those due 
to the higher order PQCD corrections, since at small momentum transfer 
Q, the logarithmic behavior of the PQCD corrections resembles the power 
behavior of the HT. It is therefore very important, when considering high 
precision experiments in the DIS region, to have an estimate of the HT. 
Already, more than twenty years ago, Politzer [20] pointed out the ne- 
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cessity of a complete theoretical description of HT [21] in DIS processes. 
In his pioneering work he gave for the first time a formal definition of the 
twist four operator of the nucleon using the operator product expansion 
(OPE) [22]. His work was then extended by Jaffe and Soldate [23l and 
soon after Ellis et al. [13] gave a different but equivalent description which 
did not require the use of OPE. More than twenty years have passed since 
those authors built the theoretical framework to compute the HT, but not 
many attempts have been done to do so [25l [26l [271 [28l [29] . The main prob- 
lem is that in order to compute the HT one must solve the QCD equations 
which are currently intractable. Most likely, years from now, people will 
be able to compute the HT effects on the lattice but, at this time, even 
the computation of the twist two PDFs presents some problem, and only 
a small number of moments of the PDFs have been computed [15]. Until 
this time will come, one has to rely on phenomenological models to estimate 
such effects. 

In this part of the study I use the MIT bag model (MITBM) to de- 
scribe the nucleon wave function and estimate the twist four contributions 
to the different structure functions (SFs). The MITBM is a phenomeno- 
logical model first proposed by the authors of Ref. [30] to describe how 
the quarks are confined inside the nucleon. Since no free quarks have been 
experimentally observed, they must be tightly confined inside the hadrons. 
The main assumption of the model is that the quarks in the hadron reside in 
a region (the volume of the bag) of true vacuum (or perturbative vacuum), 
in which they behave as free particles. On the surface of the bag acts the 
pressure due to the QCD vacuum (a continuous creation and annihilation 
of quark-antiquark pairs and gluons) which keeps the quarks confined in 
the hadron. Such a model has been extensively used to estimate different 
physical quantities (hadron masses, magnetic moments, charge radii) and 
the results are within twenty to forty percent of accuracy [31]. Keeping this 
in mind, I would like to emphasize that the main goal of this part of the 
study is not to give an exact quantitative description of the HT, but instead 
a semi-quantitative one, in which the main concern is to obtain an order of 
magnitude estimate of such effects and see if they might represent a prob- 
lematic source of uncertainty in the extraction of the weak mixing angle or 
in looking for hints of physics beyond the SM. 

In Sec. 14.21 1 present the contributions to the nucleon SFs of the two 
fermion correlation functions up to twist four (see Fig. 13. Ih . In this part of 
the study I estimate the two fermion correlations functions by computing 
their moments in the MIT bag model (MITBM) as it was originally sug- 
gested by the authors of Ref. [30], and then reconstruct them by means of 
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the inverse Mellin transformation (IMT) [32]. Indeed, given a function f{t) 
sufficiently well behaved, one could compute the Melhn transformation 

POO 

cp{z)= / dte-^f{t). (1.5) 
Jo 

The original funtion f{t) can be obtained by applying the IMT to (/'(^), 

f{t) = — dzt-^cPiz), (1.6) 

ZVrZ Jc—ioo 

where c is a positive constant. In this case I can compute the moments of 
the PDFs (for the two fermion correlation functions) or of the SFs (for the 
four fermion correlation functions) from the relation |13j 

M{N)= [ dxx^~^F{x), (1.7) 
Jo 

where in this case is a positive integer and F{x) is either a PDF or a 
SF. Unfortunately, by using the MITBM for the four fermion correlations 
(Fig. 13. 2p 1 am only able to compute a very limited (three to be exact) of 
moments, and the process of inversion provides a non unique solution (see 
Sec 14.31 for details). Theoretically, one would need an infinite number of 
moments, but it turned out that, for the case of the two fermion correlation 
functions, five moments were sufficient to provide a stable solution. 

To have an idea of the order of accuracy of the model, 1 also computed the 
twist two PDFs. These PDFs are normally twice as big as the one extracted 
from data. As pointed out in Ref. [33], theses PDFs are obtained at the bag 
energy scale (which is of the order of few hundred MeV), which becomes 
a parameter of the model, that is fixed by evolving the PDFs at an higher 
energy scale by using the DGLAP evolution equations, and comparing them 
with the one extracted from experiments at that energy scale (see Sec. I4.2.3p . 
After evolving the PDFs obtained in the MITBM, 1 have a quantitative 
agreement in magnitude, but only a qualitative agreement in shape, with 
the ones extracted from experiments (see Sec. 14. 2p . However, since my 
intent is to provide a semiquantitative estimate of the HT contributions, 1 
consider the accuracy of the model adequate for my purpose. 

In Sec. [33] I use a variation of the MITBM ( modified MITBM or 
MMITBM henceforth) suggested by the authors of Ref. [33] . In this model 
the authors overcome one of the major problem of the MITBM, namely the 
lack of translational invariance of the system, by application of a Peierls- 
Yoccoz projection [3l]. This model has the advantage of reproducing the 
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twist two PDFs quite well, but at the expense of a larger number of parame- 
ters which need to be fitted. Moreover, in this model, it is extremely difficult 
to compute numerically the four fermion correlation functions, which I was 
able to estimate only in the MITBM. A common problem to both methods is 
the absence, at the moment, of the evolution equations (the equivalent of the 
DGLAP equations) for the twist four operators. In this study I compute the 
HT corrections to the SFs at the bag model scale, and compare the results 
with some model independent extractions from experiments [181 119j. It is 
noticed that, if one trusts these extractions, the magnitude of the HT effects 
obtained with my models is too large. This could be an indication that also 
for the twist four case, evolution might play an important role. To have 
an idea on how the twist four contribution might change upon evolution, I 
make the ansatz of using the DGLAP evolution equations to evolve the two 
fermion twist four correlation functions. This procedure is quite arbitrary 
and I do not have any physical justification for applying it, nevertheless I 
noticed that, after evolving the two fermion twist four correlation functions, 
the corresponding results were in reasonable agreement with some model 
independent extractions (see Sec. I4.4.2p . I only applied evolution to the two 
fermion correlation functions, since the software I used [HT], requires that I 
introduce the distribution for each parton. For the four fermion correlation 
functions the partonic concept is lost (see the form of these contributions to 
the SFs in Sec. 14. 3p . and the software cannot be used. It turned out that 
these contributions were small even if not evolved (see Figs. 14.331 and 14.341 
where the relative corrections to A^il for a deuteron target, arising from the 
four fermion correlation function are plotted). 

In Sec. 14.51 I present the results for Afi,l, including all the corrections 
discussed in Chapter [2] and the HT contributions. 

1.3 Bremsstrahlung Contribution to the Resonance 
Electroproduction Cross Section and Aju 

In Chapter[5]I compute the parity violating asymmetry for the bremsstrahlung 
contribution to the proton electro-excitation, i.e. the process e~ P — > e~ N* 
(where P is a proton, e~ is an electron and A'^* is a proton resonance). 
It was speculated that such a process could have been one of the major 
backgrounds in the experiment E-158 performed at SLAG |35j . 

Since the M0ller scattering consists in electron-electron scattering, one, 
normally, does not have to worry about hadronic backgrounds. But in order 
to reach the desired luminosity, atomic electrons in hydrogen atoms were 
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used as a target in the experiment. Therefore, because of the presence 
of the protons, a number of hadronic backgrounds needs to be taken into 
account. Most of these backgrounds were estimated in the proposal [35], but 
the bremsstrahlung emission was thought to be neghgible, and therefore not 
considered. From prehminary experimental results, it was observed that 
the inelastic background constitutes 40% of the total asymmetry measured 
in the M0ller ring, which implies that the prescribed accuracy of 8% could 
only be reached if all the inelastic backgrounds are known to a 20% level or 
better. 

The reason why I decided to compute the bremsstrahlung contribution 
to the proton electro-excitation can be understood after knowing a bit more 
about the detector used in the E-158 experiment. This detector might be 
thought as a series of concentric rings. The innermost is the M0ller ring 
(which is actually made up of three rings), where the M0ller events are 
detected. The outermost ring is the EP ring where most of the elastic 
and quasi-elastic events are restricted. Electrons which undergo inelastic 
scattering, might end up in both detectors. Very roughly, one might assume 
that, for 50 GeV electron beam, the events with energy between 12.5 and 25 
GeV end up in the M0ller ring, while those with energy in the range 25 to 
50 GeV are detected in the EP ring. The electron resulting from the proton 
electro-excitation, are normally confined in the EP ring. But if the electron 
undergoes a bremsstrahlung process, emitting a photon either before or after 
the scattering off the proton, it might lose enough energy so that it ends up 
in the M0ller ring. 

The main problem is that in order to compute the cross section and 
asymmetry for such a process, one needs to know the electro-magnetic and 
electro- weak form factors for the different resonances for a large range of . 
Since we have a good knowledge of the electro- magnetic form factors just 
for a limited number of resonances, and the electro-weak ones are poorly 
known for all of them, I would like to investigate how the uncertainty on the 
resonance form factors impacts the error on the total asymmetry measured 
in the M0ller ring. To this end I consider a simple model in which only 
two resonances contribute to the electro-excitation process, in particular I 
consider the P33, which is a spin | isospin | resonance, and the -D13 which 
is a spin | isospin ^ resonance, and by letting the form factors vary in a 
sensible range of value I investigate the dependence of the asymmetry on 
their uncertainty. 

Chapter [5] is organized as following. In Sec. 15.11 1 describe the E-158 
experiment in some detail, while in Sec. 15.21 1 present the calculation for 
the bremsstrahlung contribution to the e~ P — > e~ A^* process for the D13 
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and the P33 resonances to the total asymmetry. In Sec. 15.31 1 outhne the 
calculation of the loops corrections to the proton electro-excitation, which 
needs to be added to the bremsstrahlung cross section to obtain a finite 
result. Finally, in Sees. 15.41 and 15.51 1 present the numerical results and the 
conclusions. 

1.4 Electroweak Radiative Corrections to PV Elec- 
troexcitation of the A 

In Chapter [6] I present a work done in collaboration with Ramsey-Musolf, 
Zhu, Holstein and Maekawa |36j . In this study, I contributed in the compu- 
tation of the chiral corrections to the axial ^ A electro-excitation ampli- 
tude, where is a nucleon and A is the delta resonance, which could help 
the interpretation of the measurement of the axial — > A matrix element 
planned at Jefferson Laboratory [37]. The goal of the prospective experi- 
ment, is the extraction of such matrix element to a 25% level of accuracy by 
measuring the parity violating asymmetry Ajil in the range of between 
0.1 and 0.6 GeV^. Here, I considered those corrections to Aju which are of 
order 0{a) compared to the leading term, which is of order 0{Gf), which 
arises from the Zq exchange. At first glance, one might think that, in a 
25% determination of the matrix elements, the 0{a) corrections could be 
neglected. On the other hand, has it has been shown in Refs. [38^ [39} BO] 
similar corrections to the electron-proton PV amplitude, were unexpectedly 
large and theoretically uncertain. If the same thing would happen for the 
A'^ ^ A matrix elements, the extraction of the axial form factors from Aju 
would be considerably more complicated than it was originally thought in 
the proposal. 

The PV amplitude for the process eP eA is generated by the diagrams 
in Figs ll.2b -e. The tree-level SM amplitude is 



iM 



PV 




(1.8) 



where 





l^^A\Jx\N) 



(1.9) 



arise from the diagram in Fig |1.2b and 



iM, 



rPV 
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l\A\Jx5\N) 



(1.10) 
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Figure 1.2: Feynman diagrams describing resonant pion electroproduction. 
The dark circle indicates a parity violating coupling. Fig. I1.2l i gives tran- 
sition anapole and Siegert's term contributions. Fig. 11.2b leads to the PV 
d-wave vrA^A contribution. 
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from the diagram in Fig |1.2b . The quantities J\ ( Jas) and l\ (/as) denote the 
vector (axial vector) weak neutral currents of the quarks and the electron, 
respectively [2]. The diagram in Fig |1.2l i contains two contributions, the so 
called Siegert term, and the anapole term whose amplitudes are 

iM^^egert = -i^^^e-7^eA4(M - Ma)^'^^ - ^ V] , (LH) 



iM^apoie = -i^^^^n,eA,N , (1.12) 

X 

where ca (called anapole) and are low energy constants, which consist 
in a calculable part, and the counter term part which needs to be fixed from 
experiments (see Sec. l6.3l for more details), M and Ma are the masses of the 
nucleon and the A respectively, is the nucleon spinor, is the A spinor, 
and ^ 1 Gev is the chiral symmetry breaking scale. The amplitude 
arising from Fig ll.2b is the d-wave amplitude that also contributes to the 
radiative corrections. Its form is much more lengthy than the previous two, 
but can be found in Ref. One thing worth noticing, is that, because 

of the factor ^ in the Siegert PV amplitude, this term is highly enhanced 
in the low region accessed by the Jefferson Lab experiment, compared 
to the other amplitudes. This peculiarity of the Siegert amplitude, might 
actually be used to extract the counter term part of the low energy constant 
d/\, since, as it has been shown in Sec. 16.31 for reasonable values of ^a, 
the Siegert term dominates in the low region compared to the other 
contributions to Arl. 

Chapter [6] is organized as following. After a general introduction in 
Sec. 16.11 about the motivation and the general background of the study, I 
introduce, in Sec. 16.21 the definition of the various kinematical variables 
and the different quantities that appear in the following sections. Sec. 16.31 is 
devoted to presenting the results of the one loop and 1/Mn chiral corrections 
to the low energy constants oa and d/\, and in Sec. 16.41 to the counter term 
parts of OA and dA is estimated, and the corrections due to the Siegert, 
anapole and d-wave contributions to Arl, are computed. Finally, in Sec. 
16.51 1 summarize the results of the study. 



13 



14 



Chapter 2 



Complete Treatment of the 
Parity Violating Asymmetry 
at Leading Twist 



In this chapter I present the calculation for the parity violating (PV) asym- 
metry AjiL at leading twist. Recently, some experiments have been proposed 
in which a PV deep inelastic scattering (DIS) asymmetry, using deuteron 
target, will be measured to a very high degree of accuracy [HI |9] to look for 
possible extension of the standard model. The intent is to extract sin^ 9\y 
at a 0.5%, which translates into a 1% accuracy in measuring A'^jf . A prob- 
lem that needs to be considered is if our knowledge of the nucleon structure 
is good enough to allow such a high accuracy to be achieved. In the ap- 
proximate case in which sea quarks are neglected, the mass of the target M 
is negligible compared to the momentum transfer, and perturbative quan- 
tum chromo dynamics (PQCD) corrections are not included, A'^jf for an 
isoscalar target, like the deuteron, assumes a particularly simple form [TT] 



A 



RL 



Gpq^ 9 r, 20 . 2^ , , ^ . . ( I - {I - y) 



sin^ ew) + {l- 4sin^ Ow) 



2V27ralOL^ 9 ' ' \l + {I - 

(2.1 



As can be seen, under the previous assumptions, the asymmetry does not 
depend on the structure of the target. However a more complete treatment, 
which includes the hadronic corrections, is necessary if one wants to reach 
a precision of 1% in the interpretation of the measurement of A^jf . 

In Sec. 12.11 1 present an overview of the DIS experiments since their 
advent in the late 1960's, to provide a better understanding of the fun- 
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Figure 2.1: Feynman diagram of a DIS process, where an electron scatters 
off a nucleon target. 

damental role they played in expanding our knowledge of the nuclear and 
nucleon structure. I also briefly introduce the quark parton model (QPM). 
In Sec. 12.21 1 present a complete calculation of Ajil at leading twist, in which 
sea quarks, target mass (TM), and PQCD corrections are included. To do 
so I use two sets of PDFs parameterizations provided the CTEQ [1^ and 
the MRST [T7] groups. The idea is to use the difference in the prediction 
for AjiL obtained by using the two parameterizations as an estimate of the 
theoretical error due to the hadronic structure. 

2.1 Deep Inelastic Scattering and the Quark Par- 
ton Model 

The scattering of electrons off atomic nuclei has proved to be one of the most 
effective ways to probe the nucleon structure. The process normally involves 
an electron of four-momentum k = (E, k) scattering off a nucleus target at 
rest, and measuring the energy and the angle of the scattered electron, whose 
four-momentum now is k' = (£'',k') (see Fig. 12. ip . In this process a virtual 
photon with four momentum q = k — k' is transferred to the target. Since 
the wave length of the virtual photon is proportional to the inverse of the 
momentum transfer, the higher the momentum the higher is the resolving 
power that can be achieved. By increasing the Q^, different layers of nuclear 
matter can be revealed, starting from the nucleus, going to the nucleons, and 
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reaching the quarks. For instance, one could consider the scattering of an 
electron off a proton with mass M. From energy conservation, the Bjorken 
variable |41j xb = 2^1771 where P is the proton momentum and v = E — E', 
is equal to one for elastic scattering, and the form factors (FFs) and cross 
section is proportional to a delta function 5{xb — 1) due to the kinematical 
constraint. The FFs could also have a dynamical dependence due to the 
finite size of the target, which implies a different response depending on how 
deep inside we are probing the nuclear structure. This dependence for 
proton elastic FFs is known to be approximately of type 



where A^^^^^^ ~ .7 GeV^. In elastic scattering for Q'^ <C ^Iroton^ fiQ) ~ 1 
and the proton is seen as point like, and the internal structure is not re- 
solved (no dependence). On the other hand, for » ^proton, the FFs 
decrease with Q^, showing a dependence, and therefore revealing the 
proton structure. Similarly, one can imagine to describe the dependence 
of the cross section of electrons scattering off a nucleus, using a similar de- 
pendence as in Eq. (j2.2p (which is not necessarily true, but it is useful 
to describe the different layers of matter) introducing a new scale Anudeus- 
If one considers the scattering of electrons off a nucleus of A atoms and 
mass Ma, when <C Anudeus the electrons scatter off the nucleus coher- 
ently, seeing the target as a whole. By increasing Q^, the internal structure, 
i.e protons and neutrons, is revealed, and the most probable process is an 
incoherent elastic scattering of the electrons off the constituent nucleons. 
Considering the behavior of the nucleus cross section as a function of the 
variable xa = 2Ma(e-E') ' peaked in the kinematical region where the 
incoherent elastic scattering of the nucleons happens. The elastic scattering 
of the nucleus, implies the kinematical constraint ~ 1, and therefore one 
has XA — ]^ = \, where A is the number of constituents. 

In the late 1960 the first DIS experiments on proton and deuterium 
targets were performed at the Stanford Linear Accelerator Center (SLAC). 
Since the in these kinds of experiments is high, one can probe the inter- 
nal structure on the target. It was surprising to notice that the electrons 
scattered at large angles, thirty time more copiously than expected [62], re- 
vealing the composite structure of the proton which appeared to be made of 
hard point like particles, named partons. Also in this case the cross section 
seemed to be dominated by the incoherent elastic scattering of the electrons 
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on the constituents of the target, rather than on the target as a whole. The 
cross section showed a peak for a value of the Bjorken variable xb ~ ^ 
suggesting that the number of partons forming the proton was three. 

Another important feature which these experiments showed was that 
the nucleon structure functions (SFs), normally function of two kinematical 
variables chosen to be and xb, seemed to depend only on xb, and very 
mildly on Q^. This phenomenon, called scaling, was predicted by Bjorken, 
who noticed that the SFs would exhibit this property in the limit and v 
oo but keeping the ratio Q'^ jv constant j4l]- Moreover it was observed that 
the longitudinal cross section (the cross section of a longitudinally polarized 
virtual photon scattering off the nucleon) was practically zero, leading to the 
conclusion that the partons should be spin one-half particles. In order to 
explain the scaling behavior, Feynman [63] gave an intuitive picture in which 
the proton is regarded as a collection of free particles off which the virtual 
photon scatters elastically, and the total cross section is the incoherent sum 
of the individual cross section weighted by the probability q{xB) (the parton 
distribution function) of a parton of type q to carry a fraction of the proton 
momentum between xb and xb + dxB- 

It was a quite remarkable property that those partons would behave as 
free particles, even though nobody was ever able to observe one outside the 
hadrons, implying that they would be tightly bind inside them. This prop- 
erty lead to the conclusion that the interaction among the partons should 
exhibit the property of asymptotic freedom, i.e. the fact that the coupling 
constant should get smaller as the momentum transfer used to probe the 
nucleon increases. Immediately the search for a gauge theory which sat- 
isfied this property began, until 't Hoff [l2]. Gross and Wilczek |43j and 
Politzer [H] realized that non-Abelian theories possessed this feature. It 
became then evident that Quantum Chromo Dynamics (QCD) would be 
the most promising candidate to describe the strong interaction. 

Since then many DIS experiments have been performed, allowing one 
to better investigate the structure of the nucleon and, ultimately, the stan- 
dard model (SM). In the late 1970's a parity violating DIS experiment was 
performed at SLAG [1] to acquire a better knowledge of the PV neutral cur- 
rent. In the recent years the technological progress allows to measure SM 
parameters to a precision that was unimaginable just a decade ago. Recently 
proposed PV DIS experiments aim to measure the weak mixing angle to a 
very high degree of accuracy, allowing to test possible scenarios of physics 
beyond the SM. 



18 



2.1.1 DIS Kinematics and Variables 

In this section I introduce the basic concepts and formahsm of DIS. Most of 
the material can in found in Ref. [32j. 

Consider a process where a lepton beam scatters inelastically on a nu- 
cleon target. In the deep inelastic regime, the momentum transfer is so 
high that the target is broken apart, and a multitude of hadrons are pro- 
duced. We are interested in inclusive experiments, namely, those in which 
only the outgoing electrons are observed, while the produced hadrons are 
ignored [55]. The scattering amplitude for the process is depicted in Fig. 12.11 
where, for simplicity, the lepton is taken to be an electron and the target a 
proton. From Fig. 12. li the following kinematical quantities in the laboratory 



frame could be defined: 
k • k' 

= cos^"*^ (-———) = scattering angle (2-3) 
|k||k'| 

k = (E, 0, 0, E) incoming electron four momentum (2-4) 

k' = (E' , E' sm9,0, E' cos9) outgoing electron four momentum (2-5) 

P = (M, 0, 0, 0) proton four momentum (2-6) 

q = k — k' momentum transfer (2-7) 

Px = P + q final hadronic system four-momentum (2-8) 

q2 ^ -q^ c^AEE'sin^ - {me <€. E) (2.9) 

= {P + qf = Af^ + 2M{E - E') - where W is the invariant mass 

of the recoiling nucleus (2.10) 

P ■ a 

u = -j^ = ( in lab frame ) = E - E' (2.11) 
g2 

xb = '2M~ ^ ''"^ Bj or ken variable (2.12) 

(2.13) 



where the range of x_b can be easily shown by noticing that for inelastic 
processes W"^ > M^, from which the constraint < 1 follows immedi- 
ately. In the previous relations the electrons are considered massless, an 
approximation which is well justified for DIS processes in which > 2 
GeV2. 

Considering the Lagrangian describing the electron-proton electromag- 
netic interaction 

jCin=[- eij{x)-f^i;{x) + eJf,{x)] A^ix) , (2.14) 
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where ip{x)j^ip{x) is the electron current, J^{x) is the proton current and 
A^{x) the photon field, and applying the Feynman rules to the process in 
Fig- one finds, in the first order approximation, the scattering amplitude 
to be 

M = ^u{k',m')j^u{k,m){X\J^{0)\PS) , (2.15) 

where m and m' are the spin components of the incoming and outgoing 
electrons respectively, \PS) is the proton state with spin S and \X) is the 
hadronic final state. The differential inclusive unpolarized cross section is 
given by [32] 



^ ' m',m,S X 



^]^l^uW^^'' (2.16) 



k-P 

2 

where a = |^ is the fine structure constant, l^y is the leptonic tensor 



and W^y is the hadronic tensor 
and 



T X] ("m'(fc')7M^m(fc)) {um'{k')-1 (k) 



(2.17) 



^/^^ = i E(2^)''^'(^^ - ^ - J^(0)|X)(X| J,(0)|P5) . 

X s 

(2.18) 

From the relation 

''^^u{k,m)u{k,m) = Ijt + me (2-19) 

m 

where trie is the electron mass, one finds for the leptonic tensor 

= k'^k, + k^kl + ^g^,. (2.20) 
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For the hadronic tensor one has 

X s 
= ^ I d'xe'^-^P\J^{x)MO)\P) , (2.21) 

where in the last step I have used translational invariance which imphes 
(P|e*^-V^(0)e-*^^-nX) = (P|e^P-V^(0)e-*^-^|X) = {P\Mx)\X) 

(2.22) 

and, since the states Xs are not observed, I have also used the completeness 
relation ~ ^- Moreover the average over the target spin is 

understood, so that for an operator O I have 

{P\d\P)^^Y.(PS\0\PS). (2.23) 

s 

Even though the hadronic tensor cannot be directly computed - since the 
proton wave function is not known - one can still use Lorentz, parity, time- 
reversal invariance and current conservation to constraint its form. From 
current conservation, the hadronic tensor has to satisfy q^W^'^ = q^W^^^ = 
(which follows from the continuity equation d^J^'' = 0). Moreover, as can 
be seen from Eq. (|2.21|) . the hadronic tensor has to be a second rank tensor 
function of the momentum transfer and the target momentum; therefore, it 
has to be built from the metric tensor g^y and products of q^ and P^. The 
most general form for W^y satisfying the previous constraints is 

W^, = MWi{xB,Q^)(-9^u + ^ 



, W2{XB,Q') ( _P-q \ ( P-q, 



MI A* 9 ^A' ; I 9 

V Q / \ Q 



(2.24) 



where the quantities Wi and W2 are the structure functions (SFs), which 
contain the information about the nucleon structure. This functions cannot 
yet be computed from first principles and they need to be measured from 
experiments. Substituting Eqs. (I2.20|) and (|2.24p in (j2.16p the unpolarized 
cross section becomes 

dd IcP'pj'^ yd 9' 
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It can be shown that, in the case in which the electrons scatter on a 
point-hke object, the SFs assume a particular simple form. If, for instance, 
one considers electron-muon scattering the result is 



Analogously, if the electrons scatters off a parton of charge and mass niq, 
one would have 



Following the parton picture suggested by Feynman, one could write the 
SFs as an incoherent sum over the elastic contributions from the scattering 
of the electrons off the constituent partons, weighted by the probability 
distribution q{x) for a parton to carry a fraction of the nucleon momentum 
between x and x + dx, with < x < 1. By using Eq. (j2.26p and the relation 
TUq = xM one finds 




(2.26) 



(2.27) 




(2.28) 



(2.29) 



W[{xb) 




2Mx'v 



) 




(2.30) 



1 



Similarly for W2 one gets 




(2.31) 



By introducing two new dimensionless SFs 




(2.32) 



and 




(2.33) 
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one sees how the parton naturally explains the phenomenon of scaling, since 
the new SFs, in this approximation, depend only on the Bjorken variable 
xb and not on Q^. Another important consequence that could be inferred 
from the previous equation is the relation [64] 

F2{xB) = 2xBFi{xB), (2.34) 

a very important result known as the Callan-Gross relation. 



2.2 Parity Violating Asymmetry in Deuteron 

Parity violating experiments have played a crucial role in testing the SM 
of electro- weak interactions since the late 1970's [1]. In the recent years, 
thanks to a phenomenal technological progress, DIS PV asymmetries can 
be measured to a level of 1%, providing opportunity, in principle, either to 
put more stringent constraints on the SM parameters, or eventually to look 
for possible physics beyond the SM. In this section I compute Abl consid- 
ering all the possible leading twist corrections such as sea quark, TM and 
PQCD corrections, and check if their uncertainties might vitiate the inter- 
pretation of the measurement to a 1% level. To estimate the uncertainty 
due these hadronic corrections I adopt two different sets of PDFs param- 
eterizations |16t I17j. The difference between the two results is taken as a 
measure of the uncertainty. 

In order to isolate the PV part in a scattering process of an electron off 
a nucleon, one defines the polarized asymmetry 

Abl ^ (2.35) 

CTB + CTL 

where o'b(l) stands for the cross section of a right handed (left handed) 
polarized electron beam scattering on a nucleon target. The cross section 
for such a process, up to order where Mz = 91.1882(22) GeV is the 
mass of the neutral Zq boson, is represented schematically in Fig. (|2.2p . 
Since the electromagnetic part is parity conserving, the cross section is the 
same for right and left handed electrons. On the other hand, the interference 
term changes sign upon flipping the helicity of the electron; therefore, the 
asymmetry isolates the PV part. In analogy with what has been shown 
before for the electromagnetic case (see Eqs. ()2.20p and (|2.24p ). one could 
define a set of tensors for each of the processes in Fig. 12.21 in the case of 
polarized electrons. 
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Figure 2.2: Polarized electron-nucleon cross section to first order in the 
Fermi constant Gp- 

The polarized electromagnetic tensor is |46] 

rn' 

1 (7^ 

= -(/c;A;^ + /cX + y9M^ (2-36) 

where A = ±1 is the helicity of the initial electron. In deriving the previous 
relation I have used the relation 

u{k,X)u{k,\) = {}^ + me)^-^, (2.37) 

In the case of relativistic electrons the spin is related to the momentum by 
the relation Su = 

f-^ Tile 

The electron interacts weakly with the Zq boson through a current 

j^{x) = V^(x)7^(5^ + 5a75)^(x) , (2.38) 



where gy = (—1 + 4sin^ 9\y) and = 
of the electron to the Zq respectively, 
the interference term becomes 



1 are the vector and axial coupling 
In this case the leptonic tensor for 



lfu=i9v + mi;u- (2.39) 
One may also want an interference hadronic tensor, given by the product of 
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the hadronic electromagnetic current with the neutral electroweak one. 



F2{XB,Q') _P-q^ \ (p _P-q^ 



+ Fs{xB,Q')'-^P''q^ , (2.40) 

where the Fi{xB, Q^) are the electro- weak interference SFs, analogous of Fi^2 
defined above. Note that the SF F3, multiplies the pseudo-tensor e^j^yap- The 
polarized cross section, neglecting terms of order is BTl 



dQdE' " Q4 



I + (1 - y - ^)i^2(x, Q^)] + 



M| 8 sin2 ^,1. cos2 ^"^^'^^ ^ 
+ (l-y-^)F2(x,Q2)±(y-^):rF3(x,Q2)]} (2.41) 

where L, i? refers to a left handed or right handed electron (corresponding 
to the plus and minus sign respectively), s = {P + k)"^ and y = Using 
Eqs. (j2.35p and (j2.4ip the asymmetry can be written as 

^HL = <^p^ (2.42) 

where 

and Gi? = 1.6639(1) x 10"^ GeV~^ is the Fermi constant measured in muon 
decay, defined as 

Gf = . (2.44) 

V2M| sin2 Ow cos2 Ow 

The two quantities W^^^^ and W^^^^^ are respectively defined as 
W^^-^ = ^{2.i[W^ + F.((l-.)-^)" 

+ g<^F^XB{l-{l-yf)\ (2.45) 
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xsMy 



(2.46) 



In order to write the SFs in terms of the quark PDFs, it is convenient to 
express the hadronic electromagnetic and electroweak neutral currents as 
sums over the different quark fields [5] 

q 

J^i^ = E aU^^^I + E 9\ql^.lm (2.47) 

where the sum runs over the u, d, s and c quark, eg is the electric charge of 
a quark of type q, and Qy and 5^ are the vector and axial-vector couplings 
of a quark q to the neutral boson Zq. In the SM one has 

8 4 

9v = 5y = 1 - 3 sin^ 0iy , 5^ = = -1 + - sin^ 

gl = g'X = -l , g'{=gl, = l. (2.48) 

Following Ref. [2], one could write the different hadronic currents in terms 
of their isospin content 

jEMT=0^^yic) jEMT=l^yi3) (^.49) 



where 



so that 



Vi'^ = \{ui^u-d-i^d) (2.50) 



jEM ^ jEM T=0 ^ jEM T=l (2.51) 



The superscript T indicates the isospin, i.e. isoscalar for T = and isovector 
for T = 1. The neutral current could then be decomposed as 

+ d='4'^+eAS7/.75S + ^AC7/.75C (2.52) 
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where the A/l can be obtained from the definitions of the by replacing 
7/x ~^ 7^75 ^'^d the couphngs are defined as 

' =9^-9'v eA=' =9\-9i (2-53) 

il=' = V^{9^+9i) eA=' = V2.{gl+9i) 

iv=9v + 9v+ 9v eA = 9l+9i + 9% 

ev = 9v- 25y - 2<7f ^A = 9l- '2g\ - 2gi . 



Notice that, in this, as opposed as in Ref. [2], I also considered the presence 
of charm quarks, since in the proposed experiment of Ref. [9], the is 
much larger than the charm quark threshold production. 

In order to compute the SFs one has to consider the following matrix 
elements 

{p\j^{x)Mo)\P) ■ 

If now we consider an isoscalar target as it has been proposed in [HI [9], one 
would need to retain only those product of currents which transform as an 
isoscalar. Therefore, for the product of two elecromagnetic currents, which 
is needed to compute the electromagnetic SFs, I havq^ 

(2.54) 

where the symbol ~ indicates that this is not an identity, but I only retained 
those terms which give a nonzero contribution once they are sandwiched 
between isoscalar states. In a similar way, the only products needed for the 



^Remember that the product of two isoscalars is still an isoscalar, the product of an 
isoscalar with an isovector is an isovector while the product of two isovectors gives raise 
to an isoscalar, an isovector and an isotensor of rank two 
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interference SFs are 

(2.55) 

The number of terms in the previous two relations could be further sim- 
plified by limiting ourselves to the case of leading twist contributions. In 
this approximation the product of two quark currents of different flavor 
can be neglected, since they give rise to twist four contributions (see Sec. 
14. 3p . The contributions to the SFs from a product of a unit charge currents 
such as q{x)^^q{x)q' {0)^^(1 + 75)g'(0) have been computed by the authors 
of Refs. [la [Ml 1^. In Refs. [HI [H] also the TM corrections have been 
included. The result is 

/,""'(«,«) ^ 4^ (2.50) 

/,"™(.B.g) = ±(||,.K) + ^^^>,.(,)) (2.58) 

where in (|2.58|) the plus sign refers to quarks and the minus to antiquarks 
since quarks couple to the Zq trough a vector minus axial current {V — A) 
while antiquarks couple through aV + A current [U] . The small case letters 
have been used to emphasize the fact that they refer to a single quark 
(antiquark) SF and TM refers to the fact that the target mass corrections 
have been included. The qi{x) is the PDF for the i quark, and I have 
introduced the quantity 

k=(l + 4xl—] (2.59) 



+ 
+ 
+ 
+ 
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and the Nachtmann variable ^ 



2XB 

1 + k 



(2.60) 



which reduces to the Bjorken variable in the limit of zero target mass. 

In addition to the TM corrections, there are other leading twist correc- 
tions to the SFs which are of first order in the strong coupling, which I call 
PQCD corrections. Such contributions to the SFs have been computed by 
the authors of Refs. [IS HH]. As opposed to the leading logarithmic cor- 
rections (LLC) to the SFs, which do not depend on the process and the 
particular SF one is considering, and therefore could be absorbed in the 
redefinition of the PDFs, these corrections are process dependent and they 
also differ depending on the specific SF. Moreover they are also renormal- 
ization scheme dependent, an important point that I discuss later. 
The contributions to the SFs Fi, F2 and F^ in the modified minimal sub- 
traction scheme (MS) are 



Qi PQCD 



asiQ) XB 



+ Tr{Fg{ 



y 



« y 

-4^(1 

y 



my) 



4-)g.(y) 
y 



(2.61) 



\XB) 



4tt k'^S^ y . 



CFFg{^)q^{y) + 2TRFG{^)g{y) (2.62) 



rqi PQCD 
J3 



(xb) 



as{Q) CpXB dy 



4tt k"^^ y 



Fg{^) - 2 - 2^ q^{y) , (2.63) 



where g{xB) is the gluon PDF, and where Fq[xB) and Fg{xb) are defined 
as 



Fq{XB) 



3 l + x 



B 



2 (1 - x)^ 



+ l{9 + bxB)-2 



+ 2(1 + 4) 



1- XB 



log XB 

I- XB 



5{l-XB)i9 + -7T^) (2.64) 



and 



Fg{xb) = (1 - 2xB + 2x|) log 



1- XB 

XB 



1 + 8xb{1-xb), (2.65) 
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where Cp = | and Tp = ^ and g{xB) is the gluon PDF. Notice that the 
function Fg{xb) appears multiphed by a factor ^ compared to way it ap- 
pears in [39], the reason being that in their expressions for the SFs, the sum 
runs only over the quarks, while in my case it runs over the quarks and 
antiquarks. Once again I have used small case letters to indicate that the 
contribution refers to a single quark with unit coupling, and the TM correc- 
tions, which affect also the PQCD corrections, have been included |52j . 

Since the SFs are physical quantities extracted from experiments, they 
cannot be scheme dependent. On the other hand, since the quantities 
/J' ^^'~'^{xb) depend on the renormalization scheme adopted, also the quan- 
tities fj^^'^ixs), and therefore PDFs qjixs), have to, so that the depen- 
dence cancels out in the SFs. There are two different approaches normally 
used to define the PDFs, which are [47] : 

1. One specifies the renormalization scheme used to compute the func- 
tions fi'' ^^'^^ (xb) and then extract the PDFs in these scheme, which 
now carry a label indicating the scheme used (such as MS, MS). In 
this case, as long as we are consistent in using the same scheme, the 
PDFs obtained are universal and can be used for any kind of process 
besides DIS. 

2. Another approach is to consider a SF that can be precisely measured, 
usually taken to be F2, and absorb all the PQCD corrections for that 
SF in the PDFs definition, namely, to any order in QCD F2 is defined 

as 

F2{xb,Q^) = xbY,4'1(''b,Q^). (2.66) 

Q 

In this case the PDFs obtained are physical quantities and do not de- 
pend on the renormalization scheme utilized. On the other hand, they 
are not universal, and if we are going to compute a process different 
from DIS, they would need some compensating factor. This scheme is 
called DIS scheme and in this case the PDFs carry the label DIS. 

I now proceed in writing the two electromagnetic and the three interference 
SFs in the case of a deuteron target. I assume isospin symmetry, which 
implies that the up (down) quark distribution (for both quark and antiquark) 
in the proton (neutron) is the same as the down (up) distribution in the 
neutron (proton), i.e. = [dF = u^). I also assume that the remaining 
sea PDFs are the same for proton and neutron, and assume symmetric 
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distribution for the quark and antiquark, i.e. s = s and c = c|. This 
translates in the following relations 

sP = s^ = sP = -s^ , c'' = c^ = cP = c'' (2.67) 

Moreover one defines the PDFs in the deuteron as 

= (2.68) 

Since, because of isospin symmetry, all the PDFs can be expressed in terms 
of the ones in the proton, I omit the suffix P from now on. The result for 
the five SFs is 

FlixB^Q') = ^(/r™(xB)+/f™(:rB) + /r^(xB)+/i™(xB) 

+ i/9(/f ™(xb) + /f ™(xb) + /r''^^''(-B) + fr'^^'^ixB) 
+ 4/9(/f ™(xb) + /r™(xB) + /r^'^^^(xB) + /r^''''(^B) 

(2.69) 

+ i/9(/r^^(xB) + /r^^(xB) + /^^''''(-i^) + /r^^''(xB) 
+ 4/9(/i™(xb) + /i™(xb) + /r'^^^(xB) + /r^''''(xB 

(2.70) 

i^i(-i.,Q^) = ^ (er^+l^) (/r^"'(-B)+/r*^(xB)+/i^™(xB)+/r^^(xB) 

, ruPQCD, N , nuPQCD, N I rdPQCD, ^ , ndPQCDf n\ 

+ /i (a^i?) + A (a^B) + /i (a;^) + A ^ [xb] j 

- wv{ft^'\xB) + fr\xB) + /r^^^^(xB) + /r^''''(xB) 

+ 2/3g^y (/f ™(xb) + /f ™(XB) + /r^«^^(xB) + /r^''''(xB) 

■^Recently, in order to reconcile the discrepancy between the SM value of the weak 
angle and the value measured in the NuTeV experiment [53], it has been suggested that 
the previous assumption, isospin symmetry and symmetric sea quark distribution, might 
not hold [541 155| . I do not consider this possibility here 



(2.71) 
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3^3 

+ 2/3<7J/|™(x^) 



(2.72) 



(2.73) 



i^3(xB,g^) = ^d=^ + |^)(/r*^(xB)-/|™(x^)+/|™(xB)-/3™(xB) 
+ /3"^^^^(xi.) - /f^'^^^(xB) + /3^^«^^(xb) - /3'^'^^^(xb) 



Let us now introduce some new quantities which allows us to rewrite Arl 
in a form more similar to the one in Refs. [U [9] 

= 2.^Fi(x^,Q2)-^ ('-^'^ 

^""'^^ 2xbFi(xb,Q2) ^ ^ 

which are both equal zero in the naive parton model because of the Callan- 
Gross relation 



Q^{XB,Q') = -[/f™(xB,Q^) + /f^^^^(xB,Q^)l (2.76) 
xb l J 

Q,.{XB,Q') = /|^™(xB,Q^)-/|^'^'' + /3*''''''''(xB,g^)-/|^^^^'^(xB,Q^) 

(2.77) 

with Qi = u, u, d, d, s, s,c,c, qj = u, d, Qi = U,U, V, V, S, S, C, C and Qj = 

p. . C(XB,Q2)+C>B,Q2) 

Hc{XB,Q ) = 2- 



K{XB, Q2) + UixB, Q2) + V{XB, Q2) + V{XB, Q^) 



RAXB,Q') = 2 SixB,Q')+SixB,Q^ 



RvixB,Q 



U{XB, Q2) + U{XB, Q2) + V{XB, Q2) + P(xb, Q2) 

2^ _ U^{xB,Q^)+VvixB,Q^) 

U{XB, Q2) + li{xB, Q2) + V{XB, g2) + ^(3:5, Q2) 

(2.78) 
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^ e u 

-i^QaQv 


Ci,d = 


-\9a9v 


C2,u = 


-\9v9a 


C2,d = 


^ p (1 
-1^9v9a 



(2.79) 



Y = ^ ^^^/^l,, (2.80) 

J- -I- l-L Lj) y TTJiH—nT) E 



1+R(XB, 




E 




}^) 


jjxbM 


1+R(XB, 




E 


R{xb,Q 


P) 


yxBM 


1+R{XB, 




E 



y = m TT (2-81) 

i y) y i+R{xB,Q 

After some tedious algebra it can be shown that 

3G Q2 ((^Ciu[l + Rc{xB,Q)]-C^d[l+Rs{xB,Q)]^Y 

" 2Wl| b + Rs{xB,Q)+^Rc{xB,Q) 

R,{xB,Q){2C2u-C2d)Y \ 
h + Rs{xB,Q)+^Rc{xB.Q)\ ^ ■ ' 

It is worth to notice that one could obtain Eq. (12. ip by keeping only the 
valence quark distributions and neglecting the TM, PQCD and the sea quark 
corrections. Under this assumption I have R^ = R^ = R = R = 0, Rv = 
y = 1 and Eq. dH} follows right away from I^M)- 



2.3 Numerical Results 

In this section I compute Ajil in the case of a deuteron target, at the kine- 
matics proposed in Refs. [51 [U]. In Ref. [S], the goal is a 1.1% measurement of 
AjiL from which sin^ Ow should be extracted with an accuracy of 0.4%. The 
kinematical region chosen is such that xb > 0.3, to reduce the contributions 
from the sea quark and from the gluons, and an average of 20 GeV^ so 
that HT contributions should be negligible (see Chap. [3]). The scattering 
angle, in the rest frame of the deuteron, is 9 = 12°. The experiment is 
performed at two different beam energies, at 36 and 39 GeV. For simplic- 
ity I choose a beam energy of 37 Gev and compute the asymmetry for a 
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Figure 2.3: Predicted asymmetry at leading twist from the CTEQ (solid red 
line) and MRST (dashed green line) parameterizations for the kinematics 
proposed in Ref. [9]. 

range of outgoing electron energies between 9.8 and 17.2 GeV, as suggested 
in the proposal. Correspondingly, the range for xb is 0.31 to 0.75, and the 
range for is 15.8 to 27.8 GeV^. It is worth noting that my expression for 
the asymmetry differs slightly from the one in the proposal. First, in the 
proposal the factor Y is not present. It turns out that, in the kinematical 
region of the experiment, this function is very close to one (to a 0.1% level 
or better), at least in the case in which HT are not included. Another dif- 
ference is that in tEqs. (j2.80p and (j2.8ip . I have a term {yxsM) / E which is 
absent in their expression. This term is usually not included in the analysis 
of DIS experiments, because of the small ratio M/E. 

Since the PDF cannot be computed from first principles, I need to resort 
to some available parameterization. I use two of the most recent, provided 
by the CTEQ [16] and the MRST [1^ groups, and take the difference of 
the values obtained for Am as an indication of the theoretical uncertainty 
on the PDFs. In Fig. 12.31 1 have plotted the two asymmetries computed 
by using the two different parameterizations. As it can be seen, the two 
predictions mostly differ in the low xb region as shown in Fig. 12. 4^ where I 
plotted the quantity 



(the factor of two come from taking the average of the two asymmetries). 



5Arl 
Arl 



= 2 




(2.83) 
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Figure 2.4: Relative uncertainty in the asymmetry defined as in Eq. (12.830 
for the kinematics proposed in Ref. [9]. 

the uncertainty is, at the most, 0.6%, which is below the aimed accuracy of 
the experiment. 

I proceed at the same way to compute the uncertainty on the asym- 
metry, for the kinematics proposed in Ref. [8]. In this proposal a single 
measurement at beam energy of 11 GeV, a scattering angle in the lab frame 
of 12.5°, and a scattered electron energy of 5.5 GeV, is considered. I actually 
consider, similarly to the previous case, a larger kinematical region, which 
correspond to an outgoing electron energy ranging between 5.5 and 8 GeV. 
The corresponding range of the Bjorken variable is 0.28 to .72 and the 
range for is 2.87 to 4.17 GeV^. By doing so, I can have a better pic- 
ture on the behavior of the twist two corrections over a broader kinematical 
range, allowing one to see where they can be minimized. 

The asymmetry obtained from the two parameterizations is shown in Fig. 
12.51 Also in this case the largest discrepancy appears in the low xb region, 
and , as it can be seen from Fig. 12.61 at an outgoing energy of 5.5 GeV, one 
has a maximum value of nearly 0.8%, which is very close to the prescribed 
accuracy of the experiment of l%.It looks, therefore, that the kinematical 
point chosen for the experiment presents the highest uncertainty, and that 
performing the experiment at higher E' would be more desirable if one's 
intent is to look for physics beyond the SM. 
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Figure 2.5: Predicted asymmetry at leading twist from the CTEQ (solid red 
line) and MRST (dashed green line) parameterizations for the kinematics 
proposed in Ref. [8] . 
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Figure 2.6: Relative uncertainty in the asymmetry defined as in Eq. (j2.83p 
for the kinematics proposed in Ref. [8]. 
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Chapter 3 

Introduction to Higher Twist 



3.1 Heuristic Definition of Higher Twist 

In this chapter I introduce the concept of twist, and write down the set of 
operators that contributes to the hadronic tensor, up to twist four. Most 
of the material can be found in a beautiful paper by Ellis et al. [13], where 
the reader is referred for more details. A somewhat easier introduction may 
also be found in Ref. |48j . 

The original definition of twist of an operator is normally given in terms 
of its mass dimension and its spin. In order for an operator to have a definite 
spin, it has to be traceless, with definite symmetry property in the exchange 
of the Lorentz indices, and local (i.e. it must depend only on one space-time 
coordinate). If these properties are respected, then its spin is given by the 
number of Lorentz indices. 

For instance, the totally symmetric, traceless, and local operator defined as 



has spin equal two. The twist of an operator is then defined has t = d — s, 
where d is the mass dimension of the operator, and s is the spin [21]. Since 
each field ip has mass dimension 3/2 and the derivative has dimension one, 
one infers that the operator in Eq. (j3.ip has twist two. It is straightforward 
to see how a similar operator, but with an arbitrary number of derivatives, 
and appropriately symmetrized, still has twist equal two. Indeed, each new 
derivative counts for a new mass dimension, but also for a new Lorentz 
index, and therefore a new spin, keeping the value of the twist the same. It 
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can be shown that an operator with twist r enters a high energy physical 
process (such as cross section) scaled by a factor (j^^ [5T] . 

As just seen, this definition might require us to consider an infinite set 
of different operators with all the same twist, which is not very practical. 
Normally, a more convenient and less rigorous definition of twist is adopted 
(which is the one I adopt henceforth in this study). In this case, the twist 
of an invariant matrix element of a bifocal operator is defined as the order 
(plus two) in (M/Q) at which it contributes in DIS processes. 

This simpler definition allows one to identify the twist of a matrix el- 
ement by inspection. For instance consider the following matrix element, 
which appears in DIS processes 

{P\^iOh^i;{Xn)\P) , (3.2) 

where \P) is the state of a nucleon with momentum P (I use the normaliza- 
tion condition {P\P') = 2E{2'k)'^ 6'^ (P — P') , where E is the nucleon energy), 
■0 is a quark field, A is a light-cone coordinate, and n is a light-cone vector 
defined as 

= ^(1,0,0,-1), (3.3) 

with P being the boost parameter with dimension of mass. In the limit P — > 
oo one recovers the infinite momentum frame, while P = M/2 corresponds 
to the nucleon rest frame. Introducing a second light-cone vector 

p^ = P(1, 0,0,1), (3.4) 

the nucleon momentum becomes 

P^=P^ + l/2M^n^ , (3.5) 

where I have chosen the system such that the nucleon moves along the z 
axis (notice that P^ = as required). Since the operator in Eq. (j3.2|) can 
only be function of A, and n^, it can be decomposed as 

(P|^«(O)7^0'^(An)|P) = MX)p^ + f2{X)M^n^ , (3.6) 

where 

/i(A) = (P|^^(0)#'?(An)|P), (3.7) 

and 

f2{x) = ^{p\rmr{xn)\p), (3.8) 
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where, in obtaining Eqs. ()3.7p and p.Sp from Eq. (j3.6|) . I used the fact that 
= = and p ■ n = 1. Since the left hand side of Eq. ()3.6p has mass 
dimension one (three from the two quark fields and minus two from the two 
states |P > and < P'\), one infers that the two functions /i(A) and /2(A) are 
dimensionless. At the end of the day, when computing physical processes in- 
volving the matrix element in Eq. (13. 6p . the two functions /i and /2 appear 
in a combination such that the term /2M^ acquires a multiplicative factor 
l/Q^, to make the /2 contribution dimensionless (as /i). Therefore one has 
that the function /2 always appears suppressed by a factor {M/Q)^ with 
respect to /i in high energy processes. One says that the matrix element 
/i has twist two, while the matrix element /2 has twist four. It is worth 
noticing that, as opposed to the previous more formal definition of twist 
for local operator, the matrix element of a bifocal operator does not have a 
definite twist. 

After this heuristic introduction of the concept of twist, I am now going to 
give a more precise and more technical definition of the twist two and four 
matrix elements and their contribution to the different structure functions 
(SFs). I closely follow the work done by Ellis et al. in Ref. [13], where the 
second definition of twist is adopted. 



3.2 Formal Definition of Higher Twist 

As seen previously in Sec. 12.1.11 all the information on the hadronic struc- 
ture is contained in the hadronic tensor. It can be shown that the hadronic 
tensor can be related to the forward virtual Compton amplitude T^y by 
means of the optical theorem j51j . 

W^y{xB, g') = ^Im[r^,(xB, g2)] , (3.9) 

where Im stands for the imaginary part, is related to the momentum 
transfer by = —q^, and the forward virtual Compton amplitude is defined 
as 

i- j d4ze^«^(P|T[j^(0)J,(z)]|P), (3.10) 

where T is the time ordered product, and the hadronic current (I consider, 
for now, the couplings of the quarks to the different bosons to be equal one), 
can be written in terms of the quark field ip'^{x) as J^{x) = ip'^ {x)^ ^ij:'^ [{)) . 



39 



It is assumed that in high energy processes, one can always factorize T^^ in 
two parts, a short distance (or hard part) which can be computed by using 
perturbative QCD (PQCD) and depends on the hard momentum q, and a 
long distance part (or soft part), which cannot be computed perturbatively, 
and contains information on the hadronic structure (factorization theorem). 
Omitting the Lorentz indices from now on, the forward amplitude, up to 
twist four, can be then written as 



T = [S{k)t{k)] + [5^(A:i, fc2)f '^(fci, fca)] + [5^.(fci, ^2, h)t''''ih, k2, k^)] 



where the four terms correspond to the processes depicted in Fig. I3.1b -c 
(two fermion correlation functions) and Fig. 13.21 (four fermion correlation 
function), respectively. The restriction to diagrams involving no more that 
four fermion lines is justified later on. In the previous equation, the square 
bracket indicates a trace over the color and spinor indices, the repeated k 
and ki imply a four momentum integration in these variables, which links 
the upper part of the amplitudes (the S"s) to the lower part (the F's). The 
upper part of the graphs can be computed perturbatively, while the lower 
part is expressed in terms of matrix element of quark and gluon fields. Let 
us delineate, without going too much into details, how such a factorization 
of the forward amplitude could be achieved. Considering the current matrix 
element in Eq. ()3.1ip . and keeping the zero-th order term in the perturbtive 
expansion one gets. 



where now the fields are in the interaction picture (indicated by the sub 
or super script /, which I omit henceforth). The time ordered product of 
the four quark fields gives, according to the Wick theorem (see, for instance 
Ref. [77]), a quark propagator from the contraction of two of the four fields, 
which is the upper part of the graph, and a nucleon matrix of the product 
of the two remaining fields, which represent the lower part, and that cannot 
be computed perturbatively. Broadly speaking, because of asymptotic free- 
dom [42l[43l[44] (i.e. the fact that the strong couphng becomes small at high 
Q^), the quarks inside the nucleon behave like free particles when probed 
by high energy virtual photons, and perturbation theory can be applied to 
compute the upper part of the diagram. 



+ 




(3.11) 
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q+ X .I . p 



q+ X2p 



XI p 




Xlp 




Figure 3.1: Feynman diagrams for DIS processes contributing up to twist 
four involving two fermions. In Fig. a) is shown the separation between the 
short distance part (S) and the long distance part (F). In a high energy 
process, the virtual photon strikes one of the quark inside the hadron. This 
virtual Compton scattering, indicated by can be computed in perturba- 
tion theory, applying the well known Feynman rules. The bottom part of 
the graphs, the F, requires the knowledge of the nucleon wave function to 
be computed, which at the moment is not available. In Fig. b) and c) the 
scattered quark exchanges one and two gluons with the remaining quarks 
inside the hadron. 




q+ 5!2p q+ X3p 

nnnnr\ — 



^'(Xl-X2)p >l(?i2-X3)p 



Figure 3.2: Feynman diagram for DIS process contributing to twist four 
involving four fermions. There are other three similar diagrams obtained by 
crossing the first and/or the second pair of fermion lines.) 



After some algebra, it can then be shown that 



[s{k)t{k)] = / dxSp{x)rp{x) 



where 

Sp{x) = [jpS{xp)], 
where S{xp) is the amputated amplitude in Fig. 13. 3^ and 



(3.13) 



(3.14) 



(3.15) 



One thing to notice is that now the short and long distance part of the 
process are linked by an integral over the light cone fraction x, instead of 
being linked by the quark momentum k. The scattering amplitude of the 
quark S now depends on xp^, which means that the quark has momentum 
parallel to p^,, that, in the massless limit (see Eq. (13. Sp ). coincides with the 
nucleon momentum. This is an assumption valid in a high energy process, 
where the interaction time between the virtual photon and the struck quark 
(which is proportional to the inverse of the energy of the virtual photon) is 
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q+xp 
> 



h xp 



Figure 3.3: Feynman diagram corresponding to the forward Compton am- 
plitude S{xp). This hard part of the process is computed using the normal 
Feynman rules. 



much smaher than the time of interaction among the quarks (or partons) 
in the nucleon. Since the partons inside the nucleon do not have time to 
interact, the system is frozen, and the photon "sees" one of the parton (the 
struck one) carrying a fraction x of the momentum of the nucleon, while the 
remaining 1 — x fraction is carried by the other partons in the nucleon, with 
no momentum in the transverse direction (transverse to direction of motion 
of the nucleon), which would arise from the interaction among the partons. 

One can proceed similarly, considering higher terms in the expansion of 
the exponential in Eq. (j3.12p . which are of the form 

dtHi{t)] , (3.16) 



with n integer, and where 

Hi{t) = J dxJ^{x)A''{x), (3.17) 

where is the quark current, and is the gluon field (color indices have 
been suppressed). In the most general case of a process involving 2F fermion 
lines and B gluon lines, the forward amplitude would then be written as 

T = J2 [ d{x}Sp,...p,^,...^, {{x}, Q')a.;; . . . c^MBrpi-PFK-.M^ d^}^ A) ^ 

F,B ^ ^ 

(3.18) 
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where Sp-^„,pp^^___f^g represents the hard part of the process involving 2F 
fermions and B gluons, rPi---P-FMi---Ms ^j^g ggft part which consists in a 
nucleon matrix element with 2F quark fields and B gluon fields (actually B 
covariant derivatives), and = 5^ — rif^p'^ is a projector that removes the 
collinear component (i.e. the component along p^) from an arbitrary vector. 
The quantity A is the QCD energy scale, that sets the scale at which the 
long and short distance could be separated, and {x} represent a set of light- 
cone fractions (fractions of the momentum carried by the parton). One 
can express the soft part of the forward amplitude, by introducing a set of 
dimensionless quantities Ti{{x}) such that 

TP^-P^'''-^'^{{x},p,n,A) = ^ A^'-^ef ■••''^'^^•■■''^(n,p)ri({x}) , (3.19) 

i 

where e^^'"^^^^'"^^ {n,p) is a polarizer, which could be a function of p^, n^, 
dpu = gfiu - PiJ.nu - rippy and e^^ = e^uaisp'^n'^ , and n is the twist of the 
amplitude Fj. 

Before going any further, let us see the meaning of twist defined in this 
way. By dimensional analysis, one can see that the forward amplitude T is 
dimensionless. This means that when multiplying the soft part by the hard 
part, which depends on Q, the terms Ti{{x}) will be multiplied by a term 
{1/QY^~'^ to compensate the mass term A'^'~^. Therefore one has that, in 
high energy processes, the higher is the twist, the more suppressed is the 
corresponding amplitude Ti{{x}). 

Let us now consider the minimum twist at which an amplitude might con- 
tribute. Since, to an amplitude with 2F fermion lines and B gluon lines 
corresponds a nucleon matrix element containing 2F quark fields and B 
gluon fields (or derivatives), the mass dimension of such matrix element 
is 3F+B-2, which comes from the 2F quark fields, each of dimension 3/2, 
the B gluon fields or derivative, each of dimension one, and the minus two 
comes from the dimension of the two states < P\ and \P >■ Therefore, the 
minimum twist at which the amplitude rPi---PFMi- -Ms ^an contribute is 

r™" = 3F + B- max{i} [Dim(ei)] , (3.20) 

where Dim(ej) is the dimension of the polarizer e^. To compute the maxi- 
mum dimension of a polarizer, the key point to notice is that all the gluon 
indices ni . . . hb of r''i -''^'^i---'^s a,re projected onto the w's. Because of the 
relation 

= (9^ - n^PlP^ = - = , (3.21) 
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where I used the fact that p ■ n = 1, the polarizer cannot have any four 
vector Pf^. in the gluon indices /xi . . . /i_B . Since the available vectors and 
tensor to build up the polarizers are p, n, d^i, and e^^, and since they 
have dimension one, minus one and zero (for both tensors) respectively, the 
maximum dimension is achieved by considering F vectors pp. in the fermion 
indices pi . . . pp, and a combination of the two tensors in the remaining 
gluon indices. In the case in which B is an odd number, it is not possible 
to use only tensors in the gluon indices, because they have an even number 
of indices, so one needs to use a vector n. Keeping in mind this particular 
case, the maximum dimension is 

max|i}[Dim(eO] = F - ^(1 - (-1)^) , (3.22) 

which implies 

Tmin = 2F + B + ^{l-{-lf). (3.23) 

This formula justifies the fact that, when considering processes which con- 
tribute up to twist four, one has to consider only the diagrams in Figs. 13.11 
and l3.2i The diagram in Fig. 13.1b . has 2 fermion lines (F=l) and zero gluon 
lines (B=0), therefore the minimum twist of the corresponding matrix ele- 
ment is Train = 2. For the process in Fig. 13.1b . one has F=l, B=l which 
implies Tmin = 2 + 1 -|- 1/2(1 + 1) =4. Similarly, one can show that the 
remaining processes in Fig. 13.11 and 13.21 correspond also to matrix elements 
whose minimum twist is four. 

It can be shown that the complete set of matrix elements necessary to 
express the three independent SFs Fi, Fi and F^ up to twist four are|^ 

• two fermion twist two correlation function 

flixB) = ll^e^^^^P\^P'^mij^{Xn)\P). (3.24) 

• two fermion twist four correlation functions 

(3.25) 

^Notice that in Ref. [13] the longitudinal SF Fl is defined as Fl = F2/XB — 2Fi, in 
Ref. [48] is defined as Fl = F2/2XB — F\ while normally, in the literature, is defined as 
Fl = F2- 2xbFi [13 [m [67]. I adopt the first definition. 
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X {P\^P'}{0)-f^rJ^-f,Dt^irJn)D!^{r]n)lP''{Xn)\P) {3.26) 
where the covariant derivative is defined as 

D^iz) = id,, - gA^{z) , (3.27) 
and the transverse component of a vector V is 

= (5^^ - pf'n" - n^Y)Vy . (3.28) 

four fermion twist four correlation functions 

TJJr V z\ - [ ^/^ i\x it,{v-x) iu(z-y) 

X {P\i)'i{'d)iie^'i{vn)i)'i' {pn)itaij'^' {\n)\P) (3.29) 
g2 /■ dXd^idv 



X (P|iJ«(0)./75tV(i"i)*''(A'n)i/«»75'/'''(An)|PB.30) 

where g is the strong coupling constant and ta = ^ a = 1 ■ ■ ■ 8 are 
the SU{3) color generators. 

The expressions of the SFs in terms of the different correlation functions 
can be found in Sec. 14.2.11 for the two quark, and in Sec. 14.3.11 for the four 
quark correlation functions. 

Before I end this chapter, let me mention that it can be shown, by using 
the equation of motion, that the operator in Eq. 13.251 can be cast in the 
following form [13l HH] 

A'TiixB) =xlj ^e*^^«(P|V;''(0)^V''(An)|P) , (3.31) 

since I have thereby eliminated the gluon degree of freedom, this form is 
more suitable to be computed in the MIT bag model. Taking the proton at 
rest, one has 

= ^(1,0, 0,-1) and p^ = y(l,0,0,l), (3.32) 
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and the two matrix elements in Eqs. p.24p and (j3.3ip can be expressed in 
terms of the "good" and "bad" components of the spinor -0 defined as 

V^+ = i(l + 7073)V' = ^(l + «3)V' (3.33) 

and 

= ^(1 -7073)V' = ^(1-«3)V', (3.34) 

respectively, so that, recalling that t^, = l/M(7o + 73) and|i = M/2(7o — 73), 
I have 

^'^^^^ = Wlj ^e''''"(^lV'?(0)V';(An)|P) (3.35) 

and 

k^T^{xB) = Mxl I ^e*^^-^(P|^!.^(0)Vl(A?i)|P) ^ xI/Uxb) ■ (3.36) 

This form for the two correlation functions will be useful in Sec. 14.41 where 
they will be computed using a modification of the MIT bag model. 
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Chapter 4 



Estimate of the Twist Four 
Contribution in the MIT 
Bag Model 



When attempting to make precision tests of QCD inelastic processes, one has 
to face the problem of including the higher twist (HT) effects. In particular, 
the twist four effects contribute to DIS processes suppressed by a factor of 



\q) ' ^^^^^ ^ ^® QCD scale, of the order of few hundreds MeV, and 
is (minus) the square of the four momentum transfer. In the attempt of 
exploring all the possible windows in which to look for new physics, some DIS 
experiments have been suggested in which the momentum transfer is of the 
order of few GeV. The idea is to extract the weak mixing angle by measuring 
AfiL to a one percent level of accuracy. Since the momentum transfer is 
not very large compared to the QCD scale, the HT effects might become 
significant, and eventually spoil the interpretation of the experiments, if not 
properly taken into account. In the past, not many attempts to compute the 
HT effects have been done [251 ESI EZl [281 129] , and I would like to provide my 
contribution on the matter. Besides, experimental data on the HT effects 
are almost inexistent, and a model estimate could be useful in planning 
future high precision DIS experiments. In this chapter I compute the twist 
four corrections to the nucleon SFs using two different variations of the MIT 
bag model. In section 14.11 1 briefly describe the MIT bag model as it was 
originally introduced in ref. [30] . In section 14.2.21 1 compute the moments 
of the twist two PDFs and reconstruct them by means of the inverse Mellin 
transform (IMT). In Sees. 14.2.41 and 14.31 1 utilize the same approach used 
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for the twist two PDFs to compute the two fermion and four fermion twist 
four corrections to the SFs. 

4.1 The MIT Bag Model 

In this section I briefly describe the main features of the MITBM in its 
simplest form, as it was initially introduced by the authors of ref [30]. The 
MITBM is a phenomenological model that describes how the quarks are 
confined inside the nucleon. Since no free quarks have been experimentally 
observed, they must be tightly confined inside the hadrons. The main as- 
sumption of the model is that the quarks in the hadron reside in a region 
(the volume of the bag) of true vacuum (or perturbative vacuum), in which 
they behave as free particles. On the surface of the bag the QCD vacuum (a 
continuous creation and annihilation of quark-antiquark pairs and gluons) 
exerts a pressure which keeps the quarks confined in the hadron. The bag 
model can be described by introducing the lagrangian density [50] 

CBag = {CQCD-B)e{\l^\-R{e,<l))) (4.1) 

where R{6,(j)) is the bag radius, B is the bag constant, which physically 
represents the external pressure exerted by the QCD vacuum on the bag 
surface, 9{x) is the step function, and the QCD lagrangian is defined as 

^QCD = '^q{x) {i^ - mg)q{x) . (4.2) 

From now on I assume massless quarks. Since the quarks are confined inside 
the hadrons, one must impose that flux of quarks through the bag surface 
is zero. If one defines the quark current as 

(q^ ^ X ^ \ 
qb{x) (4.3) 

where the indexes r, b and g are the color charges 'red', 'blue' and 'green' 
and A° a = 1...8 are the SU{3) color matrices, the condition translates as 

n^J^\surface = a=1...8 (4.4) 

where = (0, n) is a unit vector normal to the bag surface. The equations 
of motion which follow from the lagrangian in Eq. (j4.1|) for massless quarks 
are 

i^q(x) = (4.5) 
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for |x| < R, and the boundary condition (j4.4p becomes 

in^-i^q{x) = q{x)\surface (4.6) 

One could also compute the outward pressure that the quarks inside the bag 
exert on its surface. Given the energy momentum tensor 



T 



^ [q7^l'^u - q'du'JiJ.q (4.7) 
1 

the pressure is 

Pd = n^'rfT^.u ^\'^'^^^ m\surface (4.8) 

Since this pressure has to be balanced by the exterior pressure due to the 
QCD vacuum, the previous relation also fixes the bag constant B = Pq. 

Solving the equations of motion with the boundary condition (|4.6p one 
gets, for a quark of type q in the lowest energy state, 

V'™(x,t)=A^.f. '"^^'"'I'^™ )9{R-\^\)e-^^' (4.9) 
\ia ■ xji(£;|x|)xm/ 

where S = g and ~ 2.04, x is a unit vector in the x direction, a is 
a vector whose components are the Pauli spin matrices, ji{z) is the i — th 
spherical Bessel function and 

is a normalization factor. 

Since the quarks obey the Fermi-Dirac statistic, the total wave function 
(spatial-spin-flavor-color) has to be totally antisymmetric. Because of con- 
finement, the color interaction is assumed to be restricted to the bag, so 
that from outside the hadron is seen as colorless. For a baryon, which is 
made up of the three quarks, the only colorless object that can be built is 

^col = eabcq'^qS" (4.11) 

where e^fec is the Levi-Civita antisymmetric tensor which makes the color 
wave function totally antisymmetric. As a consequence the spatial-spin- 
flavor part has to be combined in such a way that results totally symmetric. 
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For a proton with spin up the total wave function is [31j 



ti|d2^3 — ujiigdg — u^d^ul — d\u2uJ^ — (fjugUg 



(4.12) 



Eq. ()4.12p is the MITBM wave function, where u and d are the up and down 
quark fields and the arrow indicates the third component of the spin. The 
spatial part of the quark fields is given in Eq. ()4.9I) and for each combination 
uud the color part is given by Eq. (j4.1ip . This wave function has been 
extensively used to compute hadron masses, magnetic moments and charge 
radii and the results are within twenty to forty percent of accuracy |31j . 
Since my intent is just to estimate the HT contributions, such a precision is 
sufficient to reach my goal. 

4.2 Two Fermion Correlation Function 

In this section I consider the contributions of the two fermion twist two and 
twist four correlation functions to the different nucleon SFs. I then outline 
the method I use to compute the moments of the SFs, and the technique I 
use to reconstruct the SFs from them. I first start with the twist two PDFs, 
so that I will be able to compare the results of my calculations with the 
available data, which will allow to test my modus operandi, and then extend 
this approach to the twist four operators. 

4.2.1 Contribution to the Structure Functions of the Two 
Fermion Twist Two and Twist Four Operators 

The contribution to the deep inelastic scattering (DIS) SFs of the two quark 
operators have been previously calculated [131 IMl HH] up to twist four Q| 



^Notice that in Ref. [13] the longitudinal SF Fl is defined as Fl ~ F2/XB — 2Fi, in 
Ref. [48) is defined as Fl = F2/2XB — F\ while normally, in the literature, is defined as 
Fl= F2- 2xbFi [47l[52l|67]. I adopt the first definition. 




(4.13) 



52 



Q 

~ i dxidx2— — — — — — —T^{x2,xi) 

J X2-X1 



(4.14) 



T'^^^ixs,Q') = ^^T.^,9^T?^0 (4.15) 



2 



Q J X2- Xl 

(4.16) 



~ i dxidx2— — — — — — —T^{x2,xi) 

J X2-X1 

(4.17) 

where the F2{xb) SF is related to Fi{xb) and Fi{xb) through the relation 
F2{xB,Q^) = xb[2Fi{xb) + Fl{xb)\. (4.18) 

The superscript EM indicates the electro-magnetic SFs, while the EMNC 
superscript labels the interference SFs, the ones arising from the interference 
of the electro-magnetic with the neutral current. The coupling Cq is the 
charge of a quark of type q while the couplings and 5^ are the vector 
and axial couplings respectively of a quark of type q to the neutral boson Zq. 
For the different couplings, I use the conventions in [2]. I have introduced 
the twist two PDF for a quark of type q as (see Chap. [3]) 

flixB) = \ j '^e^>^-^p\rmr{X)\p) (4.19) 

and the two twist four correlation function for a quark of type q as 

(4.20) 
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Z J ZTT ZTT 

(4.21) 

where the transverse component of the covariant derivative is defined as 

Dl^ = {g"'' - p^n" - n''p'')D^ (4.22) 
where and are the hght cone vectors 

= y(l, 0,0,1) (4.23) 



= 1:(1, 0,0,-1). (4.24) 



q2 

The variable is the Nachtmann variable defined as = where q is 

9 9 O"^ 

the momentum transfer and Q = —q . Moreover, xb = 57:^ is the Bjorken 
variable and P is the nucleon momentum. 

Using the equation of motion Eq. (j4.20p can be rewritten in the following 
form da SH] 

A^T^{xb) = xlj '^e'^'^'>{P\rmrW\P) . (4.25) 

4.2.2 Moments of the Twist Two PDFs in the Bag Model 

In this section I would like to illustrate how the moments of the twist two 
PDFs could be computed it the bag model and then how the PDFs could 
be reconstructed by using the IMT. Computing the twist two PDFs gives 
us also a way to test the validity of my approach, by comparing my result 
with the available data. 

Let up first consider the moments of the twist two PDF in Eq. (j4.19p |13| 



00 

M''(iV) = / dxx^-^fl{x) 







1 /""^ fix 

dxx^-^- / — e*^^(P|Vi^(0)#'?(An)|P) 

2 ./_oo 2vr 



(4.26) 
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By Taylor expanding the term 



oo ^ 

i;i{Xn)=^-{n-dfilj''{0), (4.27) 



P=oP 



I have 



5 ZTT 

1 1 

p=o P' 

oo op 



Finahy, by considering the relation 



dxx^^5{x) = i-irn\6n,p, 



(4.28) 



where in the last step I used 

r/A BP 
^^-e-A^ = (-.)^-.(x). (4.29) 



(4.30) 



I get 

W^iN) = U^~\P\^P''{0)rjl{n • 9)^-^^(0)1^) • (4.31) 

Because of the lack of translational invariance, in order to compute 
local operator matrix elements in the bag model, I need to switch from 
a continuum normalization for the state |P) with normalization {P\P) = 
2M{2Trf5''^{0) to a wave packet normalization for the state \P) = [2M{2Tr f5''^{0)]~'^/^\P),with 
normalization {P\P) = 1 [561 131]. Then, given a local operator of the type 

(P|O(0)|P) (4.32) 
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it is replaced by 

2M j dx(P|0(x)|P) . 
From Eqs. (fi:2i]) and (fOG]) I get 



(4.33) 



/ dx(P|C(x,t)(7o + 



X (ao + a3)^-V^(x,t)|P). 



73 J 



(4.34) 



By using the binomial formula (a + b)" = X^^^q (fc)^" '^^'^ ^^'^ substituting 
the expression for the bag model wave function in Eq. (|4.9p . I have 



N-l 

E 

t=o 



dx(P|V5^(x,t)(7o + 



73 J 



N - 1 



t 



N,Nl 



MN-i 

where I have introduce the quantity 



Ar-4 N~l 



E 

t=0 



N- 1 



dy 



Jo(2/)9*3io(y)+Ji(2/)^^9*, 



ya 



It, 



ji{y) 



(4.35) 



+ h{y)-dl[h{y) 



+ i dy 



,^3 

y 



My)^dl^h{y) 



(4.36) 



and I have made a change of variable y = -^x with y = ^/yf + + yh 
make the integrals dimensionless. One could see, from symmetry considera- 
tion, that the first integral in Eq. (j4.36p contributes only for t even, and the 
second integral contributes only for t odd. In fact, the two spherical Bessel 
functions jo(y) and ji{y) are even functions of y, so that the first integral is 
an even function of y for t even (zero otherwise), while the second integral 



56 



Table 4.1: Moments for the d quark twist two PDF obtained from the 
MITBM. 



Mn of ftixB) 


Ml 


1 


M2 


.167 


M3 


.042 


M4 


.013 


M5 


.004 


Me 


.002 



the opposite is true. The factor Ng comes from the spin-flavor part of the 
wave function so that = 2 and = 1. From the relation = 4 [56], 
Eq. (|4.34p can be rewritten as 

The first six moments of the twist two PDF for a d quark are shown in Table 
KT\ 

4.2.3 Moments Inversion of the Twist Two Correlation Func- 
tions 

Now that I have computed the moments of the PDFs, I attempt to recon- 
struct them by using the IMT defined as [32] 

f!M = ^ / dNx-/M\N) (4.38) 

where N has been promoted to a continuum complex variable. I follow the 
method suggested in [151 [57] in which the moments are fitted to a func- 
tion M{N) for which the inverse Mellin transformation can be computed 
analytically. I use the following parametrization to fit the moments 

M{N) = a{/3[l + c,b + N]+dl3[l + c,b + N + 1] 

+ eP[l + c,b + N + 2]) (4.39) 
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Figure 4.1: Moment function M'^{N) obtained by fitting the first six mo- 
ments in table WA\ The stars are the moments obtained from the MITBM. 



where /9[a,6] is the beta function. The inverse transform for Eq. (|4.39p is 
— / dNx-g^ a(p[l + c,b + N] 



2m 

+ dp[l + c,b + N + l]+gp[l + c,b + N + 2] 
= ax^^{l - xsYil + dxB + ex%) 



(4.40) 



which is the typical parametrization used in the literature for the valence 
quark PDFs, which are the ones one obtains from the bag model since no 
sea quark is included in it. I fit the set of parameters a, 5, c, d, e using the 
first five and six moments. The results for the two sets of parameters are, 

a = 7.582 6 =.215 c = 7.812 d = 7.532 e = 6.061 (4.41) 

for the fit with six moments, and 

a = 7.582 6 =.212 c = 7.778 d = 7.293 e = 6.155 (4.42) 

for the fit with five moments. In Figs l4.ll and l4.2I I have plotted the function 
M'^[N) for the parameterizations in Eqs. (j4.41|) and (j4.42p respectively. The 
stars indicate the value of the moments obtained with the bag model. As it 
can be seen, there is not a relevant difference between the two parameteriza- 
tions (when plotted on the same graph the two curves are indistinguishable) . 
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Figure 4.2: Moment function M'^{N) obtained by fitting the first five mo- 
ments in table 14. li The stars are the momonts obtained from the MITBM. 



One thing to notice is that the PDFs obtained from the bag model are 
independent. As suggested in ref. [53] I consider this PDF as computed 
at the bag scale usag which has to be fixed. The idea is to evolve the PDF 
obtained from the bag model according to the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) [58^ [59l [60] evolution equations, from an initial 
value of the scale Qf^ = lJ.%ag^ *o ^ specific scale Q^, and compare it with 
the PDF obtained from data at the scale Q^. The bag scale ^Bag is that 
value of Q'j^ so that the evolved PDF obtained from the bag model, best 
resembles the PDF from data at the scale Q^. In Fig l4.3l I have plotted 
the PDF XB[uvixB) + dv{xB)] obtained from the bag model and evolved 
from the bag scale fJ-^^g — GeV^, to = 4 GeV^. I also plot the same 
PDF combination not evolved, together with the one obtained from data at 
= A GeV^ using the CTEQ parameterization [16]. The software I used 
to evolve the PDF has been kindly provided to us by Miyama and Kumano. 
More details can be found in [61]. I choose fiBag so that the peaks of the two 
distributions have the same magnitude. As can be seen, the evolved and the 
experimental distributions are not quite the same in this very simple model. 
On the other hand, my main goal is not to reproduce exactly the twist two 
PDFs, but instead to give a semi-quantitative estimate of the twist four 
correlation functions. In this respect, even if they are not perfect, the model 
provides results which are within the scope of my calculation. 
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Figure 4.3: Comparison of the PDF combination xb[uv{xb) + dvixs)] ob- 
tained from the bag model, evolved and not evolved, with model independent 
extraction from CTEQ group. The green dashed line is the prediction ob- 
tained from the MITBM non evolved. The red solid line is the MITBM 
prediction evolved to a scale = 4 GeV^. The blue dot-dashed line repre- 
sents the CTEQ parameterization at the scale = A GeV^. 



4.2.4 Moments of the Two Fermion Twist Four Correlation 
Functions 

As can be seen from Eqs. (I4.14M.1T1) . there are two two-quark twist four 
correlation functions contributing to the different SFs, namely Ti{xb) and 
T2(xi, X2). I use Eq. ()4.25p since it expresses the correlation function Ti{xb) 
in a form that makes straight forward the computation of the moments and 
the process of inversion. In this form I managed to eliminate the gluon 
degrees of freedom in the covariant derivative, which I would not be able to 
include since they are not present in the MITBM. 
Defining the matrix element f^ixs) as 

flixB) ^ I ^e*"^'«(P|V;^(0)^V''(A)|P) , (4.43) 

one can compute its moments by following a similar derivation to that which 
leads to Eq. (|4.37p . The results is 

M!iN> ^ N,M^ (i)" 1' ; ') an . 
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Table 4.2: Moments for the d quark twist four PFD obtained from the 
MITBM. 



Mn of fiixB) 


Ml 


.880 


M2 


.293 


M3 


.110 


M4 


.044 


M5 


.018 


Me 


.008 



where M^{N) indicate the IMT of fi{xB) and where 



Jt 



dy 



Jo(y)9* jo(y) + ji(y) 



(J,, 



h{y) 



- 3i{y)-dy, 



ji{y)— 
y 



+ 



dy jo(2/)5*3 ( ji{y)^ 



h{y)jdl^joiy) 



(4.45) 



The correlation function T^{xb) is then related to flixs) through the rela- 
tion (see Chap. [3]) 



A^T\{xB)=xlfl{xB) 



(4.46) 



In Table IT2] are shown the first six moments of /^{xb) 

I use the same approach I used for the twist two PDF. I fit the moments 
to a function of type (|4.39p and reconstruct the correlation function /f (x^) 
by using Eq. (j4.40p . I report the result from the fit using only five moments 
and not the one with six moments, since either the moment function and 
the IMT, were indistinguishable in the two cases. 



25.778 b = 2.154 



7.445 (i = 57.183 e = 13.848 



(4.47) 



In Fig l4.4l I have plotted the moment function M'^{N) for /^(xb)- In 
contrast with the twist two case, a set of evolution equations, equivalent 
to the DGLAP equations for the twist two, is not available for the twist 
four case. As will be seen in Sec. 14.4.21 the HT contributions to the SFs 
obtained with the MITBM compute at the bag scale, are large compared 
with the model independent extractions of Refs. [18\ I19| . If one considers 
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Figure 4.4: Moment function M'^{N) relative to the twist four distribution 
fiixs) obtained by fitting the first six moments in table [121 The stars are 
the moments obtained from the MITBM. 

reliable these extractions, it might be inferred that evolution could play an 
important role, reducing the magnitude of the HT effects at high Q^. To get 
a sense on how the HT might change upon evolution, I make the ansatz of 
evolving them by using the DGLAP equation. It will seen in Sec. I4.4.2l that. 
although evolving the HT is just a guess, the magnitude I obtain is closer 
to the model independent extractions. I therefore evolve the distribution 
xb fiixs) from the bag scale ^Bag I previously found for the twist two case, 
to a scale = 4 GeV^. This last value is quite arbitrary, as long as it is 
greater than a couple of GeV^. The reason is that by using the DGLAP 
evolution equations, the evolved quantity shows scaling property [HIIBH] and 
therefore it exhibits a smooth dependence for ^ 2 GeV^. In Fig l4.5l 
I show the distribution function xb ff{xB) in the two cases in which it has 
evolved and in the case in which it as not. As opposite as to the matrix 
element T^(x), for which I was able to eliminate the gluon degree of freedom 
by using the equation of motion and compute its moments in the bag model, 
for T2{xi,X2) such elimination is not possible, and its moments cannot be 



computed directly in the bag model. However I can use an integral relation 
which relates T^{xb) to T2{xi,X2) [25] 



The authors of ref. [25] suggested expressing the correlation T^{xb) in terms 




(4.48) 
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Figure 4.5: Twist four correlation XBff{xB)- The blue dot-dashed line 
refers to the correlation function computed at the bag scale. The solid red 
line shows the result of my ansatz, where the correlation function has been 
evolved using the DGLAP equation. 



of the twist two PDF through the relation 



(4.49) 



where A; is a free parameter to be fitted, and then choose T2{x\,X2) in the 
form 



T|(xi,X2) = A5{xi - X2)Tf{xi) 



(4.50) 



so that Eq. (|4.48p is automatically satisfied. The presence of the 5{xi — X2) 
implies that the gluons exchanged in the process depicted in Fig. I3.H -C. 
are soft, i.e, they carry a small fraction of the nucleon momentum. This 
simple form for T^^xb) and T2(xi,X2) represents an educated guess which 
reproduces the N dependence of the moments of the SFs as suggested by 
the authors in ref. [65]. I suggest a similar approach in which the matrix 
element T^{x) is computed in the MITBM using the IMT method, while to 
compute r2'(xi,X2) I use the relation in Eq. (I4.50|) . From Eq. (j4.50p I get 

5{x — X2) — S{x 



dxi dx2 



X2 - Xl 



r|(x2,xi) 



(4.51) 



which is useful to compute the SFs Fi and F^. One last thing to notice is 
that the operator I compute in the bag model already includes the target 
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mass corrections [13)4 To get the pure dynamic twist four contribution I 
need to subtract the target mass corrections. The term I need to subtract 
from A'^Ti{xb) is 

M'xl f %{r^) , (4.52) 

where q{rj) is the twist two PDF for a quark of type q. I subtract the mass 
corrections at a scale = 4 Gev^. Since the mass corrections depend on 
the twist two PDF which exhibit scahng behavior (i.e. the depend on 
and very mildly on Q^), the value of I choose is not critical, as long as it is 
larger than few Gev^ . The plots of the two fermion twist four contributions 
to the different SFs can be found in Figs. (|4.16ti4.2"T]) where I have displayed 
the results from the MITBM and modified MITBM (see Sec. 14. 4p compared 
with some model independent extractions of HT from data (see Sec. I4.4.2p . 



4.3 Four Fermion Twist Four Correlation Func- 
tions 

In this section I proceed as for the case of the two fermion correlation func- 
tion, and estimate the contributions to the SFs stemming from the four 
fermion operators by first computing the moments of such contributions 
and then reconstructing the SFs using the IMT. 

4.3.1 Contribution to the Structure Functions of the Four 
Fermion Twist Four Operators 

The contributions to the DIS SFs of the four fermion twist four operators 
have been previously calculated by the authors of refs. [131 1211 HH] 

Ff^^{xB) = ^^X] / dxdydz eqeq'Ui{x,y, z) 

X (A(x,?/,2:,xb) + A(y - - 

- A{y - x,y,z,XB) - Mx,y,y - z,xb)) 

+ egeg/U2ix,y,z) 

X {A{x,y,z,XB) + A{y - x,y,y - z,xb) 

+ ^{y - x,y,z,XB) + A{x,y,y - z,xb)) (4.53) 



^In Ref. [13] the matrix elements Ti and T2 including the TM are represented by 
bold-face letters 
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pEMNC 



{xb) 



q,q 



X (A(x,y,z,XB) + A(y - - z,xb) 

- ^{y - x,y,z,XB) - A{x,y,y - z,xb)) 

+ egg^,U2{x,y,z) 

X {A{x,y,z) + A{y - x,y,y - z,xb) 

+ A{y - x,y,z,XB) + A{x,y,y - z)) 



(4.54) 



rpEMNC 
^3 



A2 



^^Q2^ I dxdydz {eqgf,)Ui{x , y , 



X (A(x,y,2;, - A(?/ - - z, 

+ ^{x,y,y - z,xb) - A{y - x,y,z,XB)) 

+ egg^,U2{x,y,z) 

X (A(x,y,2;) - A(y - - 

+ A(y - x,y,z,XB) - A{x,y,y - z)) 



(4.55) 



where g and are not necessarily quarks of the same flavor. For the four 
fermion correlation functions the longitudinal SF is zero, therefore I 
have the relation F2{xb) = 2xbFi{xb) for both the EM and the EMNC 
contributions. The functions A(x,y, z, xb) are defined as 



A{x,y,z,XB) 



6{x - Xb) 



+ 



S{y - Xb) 



+ 



5{z - Xb) 



{y-x){z-x) {x-y){z-y) {y - z){x - z) 

(4.56) 



The four fermion operators Ui and U2 are 

-2 



Ui{x,y,z) 



9 
4A2 



2tt 2-k 27r 



X {P\iP'i{f))ie^'i{vn)i^'i {i^m)Tjitaij'' (An)|P) (4.57) 



U2[x,y,z) - J 27r27r27r 

X {P\^P'lm-f5e^^J'^{un)^P'^\^In)rjita-f5^P'l' {Xn)\P) , (4.58) 
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and 



n^ = — (1,0,0,-1) (4.59) 



where g is the strong coupHng constant and ta = ^ a = I ■ ■ ■ 8 are the 
SU{3) color generators. 

As for the case of the two quark correlation functions, I would like to 
compute the moments of the SFs 



M{N) = [ dxBX%-^Fi{xB) i = l,3 (4.60) 

JO 

and reconstruct them by using the inverse Mellin transform 

1 i-c+ioo 

Fi{xB) = ^. dNx-/M{N). (4.61) 

To compute the moments, it is useful to consider following integral 

ldxdyd.ldxBX-BAiy-x,y-.,y,XB) 
^ / ^^^^^''^^''^'"'^^^'^^'"'^/i(^^)/2W/3(An) 

= j dxdydz^iy- zf-^-^ Y.(y- ^y^'V 



j=l k=0 



y jn-drhjO) (n -0)^/2(0) (n • 9)'?/3(0) 
^-^ ml pi ql 

m,p,q=0 

J 2tt 27r 2tt 

where I used the following relation [13] 

. N-i j-i 

dxBX^A{x,y,z,XB) = 5;y^-^-i^x^-'="i/\ (4.63) 

j=l k=0 
and where n ■ d = n^d^. From the fact that 

/ f = ^^-6^) 
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I have for Eq. K62\i 

N^l j-1 



I dxdydz 

j=l k=0 

oo 

E ( 



J 

oo 



, p+g (^^ • g)"^/i(0) {n ■ drf2{0) (n • d)'^MO) 



y ^ • or MO) {n • 9)^/2(0) (n • d)''MO) 



m! o! g! 

m,p,q=0 

Qm Qp Qq 

^{6{z - y))^{5{y - x))^{5{x)) 

dxdydz ( 

j=i t=o ^ ^ ^ 

k=0 s=0 V / 
-i 

m\ p\ q. 

m,p,q=0 

gm QP Qq 

(4.65) 

By repeated use of the relation 

dxx^^{5{x-y)) = {-IT I 

dx'P {n — p)\ \ zero other wise 

(4.66) 

I get for Eq. (1021) 

zvr ZTT zvr 

N-l j-l j-k-lN-j-l s N~2-t-m . ,^ 

EE E E E E ^^-Vir*^^-^-' 

j=l fe=0 s=0 i=0 m=0 p=0 

- k- l\ (s\ /N -2-t-m^ 
X (n • 5)"^/i(0)(n • 9)P/2(0)(n • df-^-'^-^fm 

(4.67) 
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Similar relations for the other three A's functions can be found in Appendix 
Oand are used, along with Eq. (I4.67|) . to calculate the moments of the SFs 
(|4.53II4.55I) . 

I now consider one of the possible contributions to the moment of the 
SFs, the remaining terms can be obtained similarly. For instance consider 
the integral 

(4.68) 

From Eq. (jC.ip I have 



2 Af— 1 J — 1 J— fc— 1 j — 1— m. 

4A2 



5 ,N~2 n' - 1 - "A / i - ^ - 1 

p J \ m 



EEE E 

j=l k=Q m=0 p=0 

(4.69) 

Consider first (omitting color matrices) 
(n • dyqa{x)4j5{n ■ dYq^[x) 

= -j^[in- dnjo{E\^\)e^^'){n ■ dnUE\^\)e-^''')xla,x. 

- i{n . dY{jo{E\^\)e-'^\l){n ■ a)^'(ji(i5|x|)<T • xe^^*x<x) 
+ i{n . a)P(jo(S|x|)e^^*xl)(n • dY{h{E\^\)a ■ xe^^^x.) 

+ Xl{n ■ dnh{E\^\)a ■ xe'^')a,{n • dy{h{E\^\)a ■ xe'^^^Xo 

jSt2 ^ 

= {-ir-'^'^[{n ■ dnME\^\)e^^'){n • dy{ME\^\)e-^^') 

- i{n ■ dy{jo{E\^\)e-'^'){n ■ 9)^'(ii(i?|x|)^e^^*) 

+ i(n • 5)P(io(ii;|x|)e^^*)(n • dy{n{E\^\)^e-^^') 

|x| 

+ (n • 9)^'(ji(£;|x|)^e^^*)(n • a)«(ji(i?|x|)^e-*^*) 

- ix' + y'){n . a)P(M^e^i^*)(n • a)'?(^l^^e-^*)^ 

(4.70) 
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where I have made use of the fact that 

1 , d 



n 



d 



M'dt ^ dz' 



(4.71) 



only acts on the third spatial component. Using Eq. ()4.70p and setting 
q = N — 2 — p — m, Eq. (j4.69p becomes 



2 ,.2 N-1 j-1 j-k-lj-l-m 

^^,..(-1)"-'-'^ E E E E 



N-2 



j=l k=0 m=0 p=Q 

(io(i?|x|)e*^*)(^-5)'"(io(i?|x|)e-^*) 
i{n . 5)'"(jo(ii^|x|)e-^*)(ji(i?|x|)^e'^*) 



j — 1 — ni\ f j — k — 1 



P 



m 



+ i(jo(^|x|)e*^*)(n • dr{ME\^\)f-e-^^') 



|x| 

z 



-iEt\ 



X 



+ (ji(i?|x|)^e^^*)(n • 9)-(ji(i?|x|)— e-^*) 

X X 



2n/Ji(-^|x|) iiJt 



e^"^)(n.a)-(^^i^|^e-^*) 



{n-df{UE\^\)e'n{n-dy{UE\^\)e 
- i{n ■ 9)^(io(^|x|)e-^^*)(n • a)P(ji(i?|x|)^e*^*) 



X 



+ i{n ■ dr{jo{E\^\)e''''){n ■ a)^(ii(i?|x|)^e-^^'*) 



iEt\ 



X 



-lEt 



- {x^ + y^){n-dn 



jliE\^\)jEt 



X 



jl(,E\^\) ^_,Et. 



(4.72) 



where /Cg-o-' is a factor arising from the spin-flavor-color part of the wave 
function. See appendix O for the results arising from the Ui opeator. 

As an example let us consider the contribution to the second moment 
of F^^{xb) due to the operator U2 for a deuteron target. The reason I 
choose such a kind of target is that recent proposed DIS experiments might 
eventually use it [319]. From Eq. (lOSll . for = 2, I have 



j dxB XB A(x, y, z, xb) = 1 



(4.73) 
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for each A, therefore from Eq. ()4.53p I get (only U2 contribution) 



q,q' a^u' 

X /^y[jo^(|y|) + ^^^^^p^^-?(|y|)], (4.74) 

where I made a change of variable y = i?x to make the integral dimension- 
less, and y = y^y^ + y| + (the extra factor of 2M is due to the different 
normalization when I go from continuum normalization to wave packet nor- 
malization. See Sec. I4.2.2p . Since 

1 

Ar2 = _1 M75) 



and as = 1^ and = 4 [56], using the spin-flavor factors from Eq. IB. II 
and the color factors in tables IB.lllB.Sl I have 



M{N = 2) 



500 320 \ as ^^M 



81 81 J 467rQ2 ^5^2 _ ^^^2 
I dy[ji{\y\) + ^^^^^^J?(|y|)] . (4.76) 



The factor ^ ~ ^) comes from the charge-spin-flavor-color term in 

Eq. (|4.74p . which, recalling the form of the operator U2 in Eq. (|4.58p . is 
equal to 



q,q' rT,a' 

q,q' cr,a-' 



where \P)sFC stands for the spin- flavor-color part of the wave function. The 
first term —500/81 comes from the proton contribution, while the second 
term —272/81 comes from the neutron contribution, and the factor 1/2 is 
due to the definition of the deuteron SFs as average of the two. Let us 
compute explicitly, as an example, the proton contribution. The flavor and 
color parts for a proton wave function can be found in appendix [Bj Since I 
am working in the bag model, the sum over the quark flavors is restricted 
to u and d, hence giving the following choices for q and q': 

q = q' = u ; q = q = d ; q = u , q' = d ; q = d , q' = u . (4.78) 



70 



Obviously, the last two cases give the same contribution. The color con- 
tributions are easily computed. Prom tables IB.H IB. 21 and IB. 31 I have that 
the color factor, which is given by the sum of the all entrees in the second 
column, is -8/3 for q = q' = u, 4/3 for q = q' = d and -16/3 for q = u , q' = d 
and q = d , q' = u (remember that the repeated color indices in the SU(3) 
matrices implies a sum over such indices). For the flavor part, whose 
contribution for the different pairs are shown in Eq. (|B.ip . one needs to 
be a little careful because of the factor (—1)°^^°^ where a and a' are the 
third component of the spin of the first and second pair, respectively. If 
the two pairs have same spin component, then the factor give one, if they 
have opposite spin component, the factor is minus one. For q = q' = u the 
flavor factor is 8/3, for q = q' = d the result is 1 and for q = u , q' = d 
and q = d, q' = u I have -4/3. Finally, including the product of the quark 
charges, the factor in Eq. (j4.77p is 

.?)«_22f-15U-lUill = -™. (4.79) 

3/3 9V3yV3y 93 81 ^ ^ 

The contribution from the neutron wave function can be similarly computed, 
except that now, in the flavor and color factors, I need to exchange u with d. 

The value of the strong coupling to be used requires some attention. The 
authors of ref. [23] take it to be cts{fJ''Bag)'^s{Q'^) where ois{^i\ag) 
value of the coupling constant at the bag scale, while asiQ"^) is the value 
of the strong coupling at the I intend to evaluate the HT contributions. 
I take the strong coupling constant to be equal one, therefore my estimate 
provides an upper limit. Notice that according to the value of the bag 
constant I found in Sec. 142:21 ~ .2 GeV^, for a ^ i QeV^ and 

for a value of Aqcd = -25 GeV (value I used for the evolution software 
provided by [HT]), the coupling constant would be of the order of as — -9, 
decreasing to ~ .7 for a = 4 GeV^. Taking a value of = 1, 
performing the integral and taking M = .938 GeV and Q = 1 Gev, I get 
the result in the table [33TT1 The complete set of results for the moments of 
the different SFs are summarized in Tables 14.3.11 14.3.11 All the results are 
for Q = 1 Gev. It must be noticed that, because of the crossing symmetry 
in the virtual forward Compton amplitude, one knows that for Fi only even 
moments contribute while for only odd moments contribute [66] (which 
is true in general, not only for the four fermion twist four SFs). Moreover, 
since Eq. (j4.63p is zero for = 1, we have the exact result that the first 
moment of F^ is zero. 
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Table 4.3: Moment function for the SF Ff obtained from the MITBM. 
Contribution arising from the Ui operator. 



Mn of Ff^^^(xB) from Ui 


M2 





M4 


7.77 X 10-5 


Me 


2.81 X 10-5 



Table 4.4: Moment function for the SF Ff obtained from the MITBM. 
Contribution arising from the U2 operator. 



Mn of F^^^'{xb) from U2 


M2 


-35.85 X 10-4 


M4 


-3.97 X 10-4 


Me 


-0.05 X 10-4 



Table 4.5: Moment function for the SF obtained from the MITBM. 

Contribution arising from the Ui operator. 



Mn oiF^^''^{xB) from Ui 


M2 





M4 


1.97 X 10-4 


Me 


0.71 X 10-4 



Table 4.6: Moment function for the SF Ff obtained from the MITBM. 
Contribution arising from the U2 operator. 



Mn of F^^^^^{xb) from U2 


M2 


-36.15 X 10-4 


M4 


-4.00 X 10-4 


Me 


-.05 X 10-4 
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Table 4.7: Moment function for the SF F^^'^^ obtained from the MITBM. 
Contribution arising from the Ui operator. 



Mn of F^''''^^{xb) from Ui 


Ml 





Ms 


-31.72 X 10^"^ 


Ms 


-6.30 X lO--^ 



Table 4.8: Moment function for the SF F^"^^^^ obtained from the MITBM. 
Contribution arising from the U2 operator. 



Mn of F^''''^'-\xb) from U2 


Ml 





Ms 


18.97 X 10"^ 


Ms 


3.21 X 10"^ 



4.3.2 Moment Inversion of the Four Fermion Twist Four 
Correlation Functions 

To compute the contributions of the four fermion twist four matrix elements 
to the SFs, I use the same approach previously used for the case of the two 
fermion correlation functions. As before I attempt to reconstruct SFs by 
means of the IMT 

1 i-c+ioo 

Hxb) = ^ / dNx-s^M{N) (4.80) 

The problem now is to numerically invert the previous relation by using 
the only three moments that I obtained. Because of the small number of 
moments available, a complete and reliable numerical reconstruction is not 
possible. As before I fit the moments to a function of type 

M(iV) = a{±p[l + c,b + N] + dp[l + c,b + N + l] 

+ ef3[l + c,b + N + 2]) (4.81) 

The plus sign is used to fit those moments for which the lowest one is non 
zero and the minus sign for those for which the lowest moment is zero. In 
this case, since the number of parameters is larger than the number of points 
used for the fit, the solution is not unique. Indeed, it happens that totally 
different sets of parameters fit the moments equally well. I also consider 
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an alternative type of function to fit the moments to test how much the 
inversion processes depends on the function chosen. The function is 



i+2 

^(^) = E(^)' (4.82) 

i=j 

where j is any integer. The corresponding Melhn Transform is 

Because of the presence of the term (logx^)*"^ in the IMT, I restrict the 
calculation to the case j = 1 and j = 2, to avoid unnatural large values of 
the quantities XBFi{xB) in the limit in which xb goes to zero. 

I present the results for the fits for two of the possible contributions to the 
SFs due to the operator U2 (the rest of the results for the other contributions 
of Ui and U2 to the SFs are presented in the appendix). As an example I 
consider the contributions to the moments of F^^\xb) coming from U2 
operator. In this case, for a parametrization of the kind of Eq. (j4.8ip . I 
have found two sets of parameters {a, h, c, d, e} which fit the moments with 
the same degree of accuracy. The first set is 

a = -302.360 x 10""^ b = 1.545 c = 9.362 d = 265.331 e = 3.006 

(4.84) 

The second set is 

a = -77.404 X 10"^ 6 = -.112 c = 7.477 d = 5.447 e = 179.207 

(4.85) 

The set of parameters relative to the parametrization of the type (j4.82p are, 
for J = 1 

ai = 8.709 X 10"^ 02 = -48.282 x 10"^ 03 = -36.746 x 10"^ 

(4.86) 

and for j = 2 

a2 = 59.569 x 10^"^ 03 = -433.076 x 10"^ 04 = 431.403 x 10^^ 

(4.87) 
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Figure 4.6: Comparison of the different fits for the moment function M{N) 
relative to the F^^'^xb) SF {U2 contribution). The sohd Ught blue line is 
relative to the set of parameters in Eq. (j4.84p . the long dashed blue line is 
relative to the set of parameters in Eq. (j4.85p . the dot-dashed red line is 
relative to the set of parameters in Eq. (j4.86p . and the short dashed green 
line is relative to the set of parameters in Eq. (|4.87p . The three moments 
obtained from the bag model are represented by stars. 




0.2 0.4 0.6 0.8 

X 

Figure 4.7: Reconstruction of the SF xbF^^^ {xb) through the use of the 
inverse Mellin transform technique {U2 contribution). The type of line used 
for each parametrization is the same as in Fig l4.61 
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In Fig. 14.61 1 have plotted the different moment functions M[N) relative to 
the different set of parameters (Eqs. (I4.84M.871) ). The stars represent the 
points obtained from the bag model. It can be seen how all the different 
parameterizations fit the data quite well and, in this case, the correspond- 
ing curves are very similar in shape and magnitude. On the other hand, 
after computing the inverse Mellin transform, those different parameteriza- 
tions give rise to different xb dependence for the SFs as one can see from 
Fig l4.7[ Nevertheless, although each parametrization corresponds to a differ- 
ent shape, the overall magnitude is quite similar. I would like to emphasize 
that I am not interested in an exact calculation of the xb dependence, since 
it is already known, from section [4.2.31 that the simple bag model does not 
reproduce such a dependence for the twist two PDFs. What we would like 
to investigate, instead, is if the magnitude of the contribution to the SFs of 
these four fermion twist four corrections is indeed negligible compared to the 
twist two contribution (or other twist four contribution stemming from two 
fermion correlation functions). All the contributions from the four fermion 
correlation functions have not been evolved, since, as will be seen in Sec. 
14.51 their contribution is already small even if they are not evolved (note 
that evolution decreases the magnitude of the function evolved). 

As a second example I consider the contribution to the SF F^, coming 
from the operator U2 ■ In this case I was only able to find one set of param- 
eters relative to the parametrization ()4.8ip (that does not mean that there 
is only one nor that the one I found is the best one). The set of parameters 
for the different parameterizations are 

a = 41.369 X 10"^ 6 = -.821 c = 5.514 d = 30.754 e = 43.527 

(4.88) 

01 = -26.773 X 10"^ 02 = 195.114 x 10"^ 03 = -168.341 x 10"^ 

(4.89) 

02 = -45.854 X 10"^ 03 = 447.451 x 10"'' = -401.604 x 10"^ 

(4.90) 

The different moments functions are plotted in Fig. 14.81 In this case, al- 
though they all fit the points obtained from the bag model calculation, they 
have quite different shape (at least in the range of values between = 1 
and A^ = 3). However, also in this case, the magnitude of the inverse trans- 
form is quite similar (see Fig. 14. 9p . The parameterizations for the remaining 
terms of the SFs are in the appendix lAl 
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Figure 4.8: Comparison of the different fits for the moment function M{N) 
relative to the F^^^'~^ {xb) SF {U2 contribution). The sohd blue line is 
relative to the set of parameters in Eq. (j4.87p . the dashed green line is 
relative to the set of parameters in Eq. (j4.88p and the dot-dashed red line is 
relative to the set of parameters in Eq. (|4.89|) . The three moments obtained 
from the bag model are represented by stars. 




Figure 4.9: Reconstruction of the SF xbF^ {xb) through the use of the 
inverse Mellin transform technique {U2 contribution). The type of line used 
for each parametrization is the same as in Fig l4.81 
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4.4 Two Fermion Twist Two and Four Correlation 
Functions in the Modified MIT Bag Model 

In this section I try to use a modified version of the bag model (MMITBM) 
suggested by the authors of Ref. [33] 5 to compute the twist two and twist 
four contributions to the SFs. In this case the correlation functions (j4.19p 
and (|4.2U|) are directly computed, without calculating the moments and then 
obtaining the SFs through an IMT. The same method applied to the four 
fermion twist four correlation functions, involves the numerical evaluation of 
multidimensional integrals, which I was not able to compute in a reasonable 
amount of time. 

4.4.1 Modified MIT Bag Model 

In the MIT bag model, since the quarks are confined in the bag, translational 
invariance is lost. The authors of [33] suggest using the Peierls-Yoccoz [34] 
projection, to ensure that the proton wave function is a momentum eigen- 
state, thereby making the system translational invariant. They define the 
action of the quark field on a n— quark state as 

'(/'(x)|xiX2 . . . x„) = (5(x — Xi)|x2 . . . x„) + permutations. (4.91) 

A translational invariant n-quark state with momentum p„ is introduced in 
coordinate space as 



where the normalization factor (pipn) is given by the normalization condition 




(4.92) 



(p„|p'J = 2M(2^)35(p„ - p', 



n 



(4.93) 



where M is the target mass, from which one finds 




(4.94) 



Recalling the form of the twist two PDF in Eq. ()3.35p 




(4.95) 



and inserting a complete set of states 




(4.96) 
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q 



^ > 

4 » 

Figure 4.10: Two quark intermidiate state process contributing to the quark 
twist two PDF. 

I get 

m J \ } 

(4.97) 

where the plus component of a four- vector is defined as v'^ = vq + v^, 
I a) is the spin- flavor part of the nucleon wave function, and (a|P^^|a) is 
the projector operator onto flavor q and spin m for a n quark intermediate 
state. In the last step in Eq. (14.970 I have used translational invariance, 
such that 

(p„|V+(An)|p) = (p„|e-^^->+(0)e*^^-"|p) = e"*^'^ (p„|V+(0)|p) (4.98) 

where P is the momentum operator, and where I used the relation v-n = 
Since the quark field operator ^^+(0) can annihilate a quark in the position 
0, but also add an antiquark in the same position, there are two possible 
intermediate states one with two quarks and one with three quarks and 
an antiquark (see Figs. 14.101 and 14. lip . The PDF for an antiquark can be 
obtained from Eq. (|3.35|) by considering the operator instead of '(/'+ [33] > 
which yields 

fti-B) = ^,Y.WJ») j |^5(-B-^^^^)KP.|VJ-(0)|P)|^4.99) 



79 




q 
q 
q 



q 



Figure 4.11: Four quark intermidiate state process contributing to the quark 
twist two PDF. 



q 
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Figure 4.12: Two quark intermidiate state process contributing to the anti- 
quark twist two PDF. 



80 



Notice that, in this case, the operator "0^(0) either creates a quark or an- 
nihilates an antiquark. Since in the MIT proton wave function there are 
no antiquarks to be annihilated, the only contribution comes from adding 
a quark to the proton, giving rise to the four quark intermediate state rep- 
resented in Fig. I4.12[ Inserting a complete set of coordinate states in Eq. 
(14.971) one gets 



X y dxidX2(fx3(p2|^+(0)|xiX2X3)(xiX2X3|p)|^ 



[l j (ixidx2(ix3(p4|^+(0)|xiX2X3)(xiX2X3|p)p] (4.100) 



Let us consider part of the second term in the previous equation. Working 
in the target rest frame so that |p = 0), I get 



(ixi(ix2(ix3(p4|^^J.(0)|xiX2X3)(xiX2X3|p = 0) 



(ixi(ix2(ix3(p4|xiX2X30j)(xiX2X3|p = 0) 

(ixi(ix2(ix3- 



X 



j (iRe~*P4-^V^ (xi - R)^"^ (x2 - R)V'^ (x3 - R) V'? ^ (-R) 
y"(iR>(xi-R')^(x2-R')^(x3-R'), 

(4.101) 



where I used Eq. (|4.92p for a three quark state with momentum p = and 
a four quark state with momentum p4. By making the change of variables 
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R — > — R and Xj + R ^ Xj (z = 1, 2, 3), I have 



dxidx2dX3(p4|V'+(0)|xiX2X3)(xiX2X3|p = 0) 



dxi(ix2(ix3- 



dRe*P4-V(xi)V'^(x2)V^(x3)V'f (R) 



>3(0)<^|(P4 

X y (iRV(xi-R'-R)V'(x2-R'-R)V'(x3-R'-R 



(4.102) 



By shifting the variable R' + R ^ R' and using Eq. (|4.94p I arrive to the 
result 



(ixi(iX2dx3(p4|V'+(0)|xiX2X3)(xiX2X3|p = 0) 



where I have defined the quantity 



(4.103) 



(4.104) 



The spinor 'i/'+(P4) is obtained by considering the antiquark bag wave func- 
tion 



^(^)(P4)= / (iRe^P^-^'0i^^(R). 



<(x,f) = iV„ 



-ia ■ xji{E\x\)xr^ 

Jo(-E|x|)Xm 



e{R-\x\)e 



-iEt 



(4.105) 



A similar relation can be found for the first term in Eq. (|4.100p . so that I 
have 



ftixB) = 2MY,{a\PlM f ^m{l-XB)-pt 



02 (P2) 



+ 2my,{<^\pIm j -^mi.^-XB)-pi 



04 (P4) 



Similarly one can compute the antiquark distribution, obtaining 



ff{xB) = 2Mj2{c^\PlJa) I ^SiM{l-XB)-pt) 



(2^)2 



03(0) 



l^f(P2)|^ 



lA'^(P4)P 

(4.106) 
(4.107) 
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Using the bag model wave function (j4.12p . one finds the fohowing rela- 
tions [331 EE] 



l0n(Pn)P 



Ann"- 



[2(02 _ sin2 0)]" \ n J 
where u = |pn|-R and the function T{v) is 



— — I / dvv 



Tiv) 



1 1 — V I sm2v- 



1 sin 2n 
2^ 20 



'•'~"T"(^;)io(^) (4.108) 



1 sin 2Q sin^ O 



^°^'"+2 + - 20 



02 



(4.109) 



Moreover, one has 



Vr^(p„)=Vr(Pn)=#^(-Pn) 
TTi? 02 



si{u) + 2^siiu)s2iu) + si{u) (4.110) 



2 (02-sin2 0) 
where si{U) and 52(1^) are defined as 



iPnl 



Si{u) = - 
U 



sin(n — 0) sin(u + 0) 
— O + 



(4.111) 



and 



u 



S2{u) = 2jo(0)ji(n) - ^si(u) . 



(4.112) 



In Eqs. (j4.106p and (j4.107p one can perform the angular integration 



//■oo 
dpnS{Mil-XB)-pt)=^T^ / 
J|M2(1- 



'|M2(l-xs)2-A/2|/(2A/(l-xs)) 
where the orthogonal component is understood to be 

pi = 2M(1 - xb)VM^ + pI- M\l - XB? - Ml 



d\vn\ Ip44-113) 



(4.114) 



and Mn is the mass of the n quark intermediate state. 

One thing to notice is the difference in some of my relations with the ones 
in Ref. [33] . For the four fermion intermediate states the ratio in the normal- 
ization factors is reversed, i.e. they have \4'a/'P^\, while I have \(j)?,/(f>i\i but 
the whole calculation that leads to Eq. (|4.103p with the normalization factor 
reversed is not shown. In the rest of the calculation I adopt Eq. (j4.103p . 
As observed in Ref. [33], the Peierls-Yoccoz projection is a non relativistic 
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approximation, which is vahd only for |p„| < M„. As a consequence, the 
integrand in Eq. ()4.113p has an upper hmit cut off equal to M„. 
The valence distributions have to satisfy the conditions 

where 

QvixB) = flixB) - fHxB) ■ (4.116) 

In Ref. [33] the normalization condition cannot be satisfied. The contribu- 
tion stemming from the four fermion intermediate state counts only for a 
nearly 2% of the normalization, while the contribution from the two fermion 
intermediate counts, at the most, for 82%. In order to saturate the condi- 
tion in Eq. (j4.115p they added a term ad hoc of the form (1 — xbY ■ With 
my normalization, I am always able, by varying the parameters, namely the 
radius R, the two masses M2 and M4 of the two and four quark intermediate 
state, and the bag scale HBag (see Sec. 14.2.31 for more details on how to fix 
this last parameter) in reasonable range of values, to fulfill the normaliza- 
tion condition, and obtaining a valence distribution that resembles the one 
extracted from data (see Figs. 14.131 and I4.14p . 

One of the problems of this model is that it yields a u quark valence 
distribution such that = 2d^, which is known not to be the case. The 
model could be improved to accommodate the different xb dependence in the 
two distributions by assuming that the two quark intermediate state could 
form either a spin singlet or a spin triplet state with different masses |70j . 
I am not trying to reproduce perfectly the twist two PDFs, but instead to 
estimate the twist four contributions; therefore, for my purpose, I consider 
the simplified model. In Fig. 14.131 and 14.141 1 have plotted the distribution 
function xb[u^(xb) + dv{xB)] obtained from the MMITBM before and after 
evolution, and compared it with the CTEQ parametrization for two different 
sets of parameters. In Fig. ()4.13p the parameters are: 

R = lim M2 = 700 MeV M4 = 1111.6 MeV HBag = -45 MeV. (4.117) 

In Fig. (I4.14P the parameters are: 

R = 0.9 fm M2 = 750 MeV M4 = 1032.9 MeV ^Bag = -48 MeV (4.118) 

Both sets of parameters are such that the normalization condition (|4.115p is 
satisfied and, as it can be seen from the plots, they both describe the data 
quite well. 
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MMITBM I Ev. 
MMITBM I N-Ev. 
CTEQ 




Figure 4.13: Comparison of the PDF combination xsiuvixB) + ^^{xb)] ob- 
tained from the MMITBM with the CTEQ parameterization (dashed green 
hne), before (dot-dashed blue Une) and after evolution (solid red line). The 
parameters are R = 1 fm, M2 = 700 MeV, M4 = 1111.6 MeV and f^isag = -45 
MeV 
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Figure 4.14: Comparison of the PDF combination xb[uv{xb) + d^ixB)] ob- 
tained from the MMITBM with the CTEQ parameterization (dashed green 
line), before (dot-dashed blue line) and after evolution (solid red line). 
The parameters are R = .9 fm, M2 = 750 MeV, M4 = 1032.9 MeV and 
IJBag = .48 MeV 
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4.4.2 Twist Four Contributions to the Structure Functions 
in the Modified MIT Bag Model 

Having fixed the parameters of my model, I can now proceed in estimating 
the twist four contribution to the SFs. The approach I follow has been 
described in Sec. 14.2.41 I compute the twist four matrix element 



flixB) 



M 



dX 
2^' 



iXxf 



{P\ipl'{0)tl;'L{Xn)\P) 



(4.119) 



which is related to Ti{xb) by Eq. (j3.36|) . As before, since the matrix element 
T2{xi,X2) cannot be directly computed, I estimate it by making the ansatz 



r|(xi,a;2) = 45(xi-X2)rf(xi), 



(4.120) 



and I evolve the twist four function xb/Hxb) using the DGLAP equation 
(refer to Sec. 14.2.41 for details). Following the same steps that lead to Eq. 
(el61), and using the relation 



^?(Pn)=V'i'^(-Pn), 



one has 

ftixB) = 2M3j^(a|P|^|a 

m 

+ 2M3J](a|P, 



i la) 

a,m I"/ 



(2^)3 

dp4 

(2^) 



5(M{1-XB)-Pt) 
^6{M{1-XB)-Pt) 



and 

rd 



dp4 

(2vr) 



:diM{l-XB)-pt) 



MO) 



<A4(P4 



(4.121) 

l#^(P4)P 

(4.122) 

l£^'^(P4)P. 

(4.123) 



In Fig. ()4.15p I have plotted the valence distribution function xsd^^ (xb) 
where 



qrixB) = fUxB) - fHxB) 



(4.124) 



before and after evolution for the two different sets of parameters {R, M2, M4, ^Bag}- 
As it can be seen, shape and order of magnitude are similar to the ones of 
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MMITBM I Ev. 

MMITBM I N~Ev. 

MMITBM II Ev. 

MMITBM II N-Ev. 




Figure 4.15: Correlation function XBltixs) obtained from the MMITBM 
using the two set of parmaters in Eqs. ()4.117p and (14.118P before and after 
evolution. 



the twist two correlation functions, suggesting that the HT might play an 
important role in processes involving a low momentum transfer, where the 
scaling factor is of the order of unity. 

I now proceed as in Sec. 14.2.41 to obtain the twist four corrections to 
the different SFs. In Figs. SUHMIlI plot the HT contribution to SFs F2, 
Fi and xbF^ obtained with the two models I used, namely the MITBM 
and the MMITBM, with some model independent extraction from data. In 
Figs. (I4.16IITT8|) the SFs are computed by using the non evolved correlation 
function XBfi{xB): while, in Figs. (|4. 19114. 2"T|) it has been evolved. I want 
to emphasize the fact that the HT have never been directly measured. Nor- 
mally, what is done, is to introduce a new term in the fitting function for 
the SFs of the type 

F,(XB, Q2) = q2) + :^f^ (4.125) 

where F^'^^^ is the leading twist (LT) SF including the target mass correc- 
tions, and H^'^{xb) represents the HT corrections. It is not easy to judge 
how reliable these extractions could be. First of all the function H^^{xb) 
is taken to be independent, assuming that the scaling as ^ dominates 
with respect to possible log dependence arising from PQCD corrections. 
Moreover it is not always straightforward to discern between the depen- 
dence of the twist four and the \ogQ^ dependence arising from the higher 
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MMITBM I N-Ev. 

MMITBM II N~Ev. 
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Figure 4.16: Twist four contribution to the F2 SF obtained from the MITBM 
(dashed green hne) and the MMITBM (soHd red Hne and dot-dashed blue 
Une) with no evolution, compared to the model indepenent extraction of 

m ■ 



order PQCD corrections to the LT, meaning that sometimes, what is thought 
to be HT contribution, could partially be accounted for, once higher order 
PQCD corrections are taken into account [18]. Nonetheless, this is the only 
attempt done to extract HT from experiment, and the only way I have to 
compare my results with some model independent extractions. In the next 
section I compute the HT corrections to the PV asymmetry for a deuteron 
target, using the results obtained in this section. 
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Figure 4.17: Twist four contribution to the F2 SF obtained from the MITBM 
(dashed green hne) and the MMITBM (sohd red hne and dot-dashed blue 
hne) with evolution, compared to the model indepenent extraction of |19j . 
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Figure 4.18: Twist four contribution to the SF obtained from the MITBM 
(dashed green line) and the MMITBM (solid red line and dot-dashed blue 
line) with no evolution, compared to the model indepenent extraction of 
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Figure 4.19: Twist four contribution to the SF obtained from the MITBM 
(dashed green line) and the MMITBM (solid red line and dot-dashed blue 
line) with evolution, compared to the model indepenent extraction of |19j . 
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Figure 4.20: Twist four contribution to the xbFj, SF obtained from the 
MITBM (dashed green line) and the MMITBM (sohd red line and dot- 
dashed blue line) with no evolution, compared to the model indepenent 
extraction of [18] . 
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Figure 4.21: Twist four contribution to the xbF^ SF obtained from the 
MITBM (dashed green hue) and the MMITBM (soUd red Une and dot- 
dashed blue line) with evolution, compared to the model indepenent extrac- 
tion of [18]. 

4.5 Higher Twist Corrections to Aju for Deuteron 
Target 

Similarly to what has been done in Sec. 12.31 in this section I consider the 
corrections to the LT asymmetry due the the HT contributions for the kine- 
matics proposed in refs. [HI [8]. For simplicity, I computed the LT asymmetry 
using only the CTEQ parameterization. First, I consider the HT contribu- 
tions from the two quark twist four correlation functions. In all of the 
following figures, what is plotted is the quantity 

Arl 

where is the LT asymmetry, whereas the A^[ corresponds to the asym- 
metry in which HT have been included, and the factor of two comes from 
averaging the two asymmetries. I start by considering the kinematical region 
proposed in Ref. [9], which consists in a incoming beam energy E = 37 
GeV, a scattering angle of twelve degrees, and an energy of the outgoing 
electrons ranging between 9.8 to 17.2 GeV. In Fig. 14.221 1 have plotted the 
relative uncertainty in the asymmetry for the three cases in which the HT 
were estimated by using the MITBM and the MMITBM (I considered the 
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Figure 4.22: Relative uncertainty in the asymmetry defined in Eq. (j4.126p 
due to the HT two fermion contributions for the kinematics proposed in 
Ref. [^. The red sohd hne is relative to the HT contributions obtained in 
the MMITBM using the set 1 of parameters, the blue dashed line using 
the set 2 of parameters, and the green dot-dashed line is obtained from the 
MITBM. 



contributions relative to both sets of parameters, the first set (set 1 hence- 
forth) i? = 1 fm. Ma = 700 MeV, M4 = 1111.6 MeV and HBag = 0.45 MeV 
and the second set (set 2) R = 0.9 fm, M2 = 750 MeV, M4 = 1032.9 MeV 
and i-iBag = 0.48 MeV ) and evolved according to the DGLAP equations. 
It can be seen that for all three cases the relative magnitude of the HT ef- 
fects is, at the most 0.15%. In Fig. 14.231 I have instead plotted the relative 
asymmetry obtained with the two MMITBMs in the case in which the SFs 
have not been evolved. In this case the maximum value of the uncertainty 
is a bit more than 0.6%. 

Next I consider the relative asymmetry obtained with the non evolved 
MITBM. I split the energy range in two parts, which are considered in Figs. 
14.241 and 14.251 As it can be seen, the corrected asymmetry becomes unex- 
pectedly large around E' = 15 GeV. This situation is actually unphysical, 
since this singular behavior is due to the fact that the electromagnetic cross 
section, which appears in the denominator on the asymmetry, becomes neg- 
ative around E' = 15.5. 

The main cause for the cross section to become negative is the possibility 
for the HT contribution to the different SFs to be negative, in conjunction 
with the relatively large size of the HT contribution arising from my model. 
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Figure 4.23: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions for the kinematics proposed in Ref. [9]. The red sohd 
hne is relative to the HT contributions obtained in the MMITBM using the 
set 1 of parameters, the blue dashed line using the set 2 of parameters. 



when evolution is not applied. This fact might suggest that my estimate of 
the HT contribution are too large when they are not evolved, and, therefore, 
that evolution could be a determinant ingredient in the study of HT effects. 
Another possibility is that twist higher than four could become important 
and, if added, they might prevent the cross section from becoming negative. 

I now consider the corrections in Am due to the two fermion twist 
four contributions at the kinematics suggested in Ref. [8]. As before, even 
though the proposed experiment consists in only one measurement for a 
beam energy oi E = 11 GeV, a scattering angle of twelve and a half degree, 
and an outgoing electron energy of 5.5 GeV, I actually consider a wider range 
of outgoing electron energies, so that one could see at which kinematics the 
theoretical uncertainties due to the HT effects, could be minimized. The 
range is from 5.5 to 8.0 GeV. 

In this case, the enhancement in the asymmetry due to the smallness 
of the EM cross section, is present in all three cases in which the HT con- 
tribution have not been evolved, since the is much smaller than in the 
previous case. 

In Fig. 14.261 1 have plotted the relative asymmetry for the same three 
cases as in Fig. 14.221 except that now the kinematics is the one in Ref. [8]. 
Also in this case the relative uncertainty is smaller than the prescribed 
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Figure 4.24: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions computed with the MITBM for the kinematics pro- 
posed in Ref. [9] . 
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Figure 4.25: Relative uncertainty in the asymmetry defined in Eq. (j4.126p 
due to the HT two fermion contributions for the kinematics proposed in 
Ref. [8]. The red solid line is relative to the HT contributions obtained in 
the MMITBM using the set 1 of parameters, the blue dashed line using 
the set 2 of parameters, and the green dot-dashed line is obtained from the 
MITBM. 
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Figure 4.26: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions for the kinematics proposed in Ref. [8]. 



accuracy, being always less than 0.5% . In Figs. 14. 27114. 32) 1 have plotted the 
relative asymmetry in the case in which the HT have not been evolved. As 
can be seen in all three models there exists a point in which the EM cross 
section vanishes, causing the asymmetry to diverge. 

As mentioned before the reason for the asymmetry to be divergent is 
that the EM cross section can become zero, if the HT could be large in size 
and negative. Of course the cross section cannot be negative, nevertheless, 
it might happen that it becomes extremely small, causing the asymmetry 
to be large. Let us consider the HT contributions obtained with the evolved 
MITBM, which, as can be seen from Figs. 14.191 and l4. 211 it seems to repro- 
duce the same order of magnitude of the HT effects found in Ref. [19^ [T8] . I 
tried to multiplying all the HT corrections in A^J^ by a factor of two, and I 
notice that the difference between the LT and the HT asymmetry could be 
easily larger than 1.5%, reaching a difference of 20% if the HT SF are scaled 
by a factor of 2.8. In this case the electromagnetic cross section does not 
become negative, but simply small enough so that the HT asymmetry be- 
comes large. I emphasize that all I have done, is to take my model estimate 
of the HT contributions, which seems to agree with some model independent 
extraction, and scale it by a factor between 2 and 2.8. 

Notice that such enhancement in Aju, cannot happen in other models 
such as the one used in Ref. |28^ [29] , where the HT contributions are esti- 
mated by computing the first moment HT correction to the different SFs, 
and where the xb dependence is assumed to be of the same type of the twist 
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Figure 4.27: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions computed with the MITBM for the kinematics pro- 
posed in Ref. [8] . 




two SFs. For instance, the correction i^^J^r t° ^2 SF would be given 

by 



1 







xb{uv{xb) + dy{xB)) , (4.127) 



where the first moment of F2, 

nl 

dxFi^'^ix), (4.128) 



is computed in the MITBM. For such a form of the HT, since all the SFs 
in the asymmetry formula in Eq. (j2.42p . including the HT corrections, will 
have the same xb dependence, such a dependence will cancel out in the 
asymmetry, yielding a xb independent shift in the asymmetry. 

It seems therefore clear that, in order to perform high precision PV DIS 
experiments, it is absolutely necessary to have a very good knowledge, not 
only in the magnitude, but also in the sign and xb dependence of the HT, 
which is the main cause for the electromagnetic cross section to become 
small and give rise to an enhancement in Apti- 

I next consider the relative uncertainty arising from the four fermion 
twist four correlation functions. The results for the kinematics of Ref. [9] are 
presented in Fig. 14.331 As can be seen, for all the different set of parameters 
I use to compute the IMT, the relative uncertainty is of the order of 0.1%. 
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Figure 4.28: Relative uncertainty in the asymmetry due to the HT contri- 
butions computed with the MITBM for the kinematics proposed in Ref. [8]. 



This result could be considered as an upper limit, since evolution, which 
decreases the magnitude of the SFs, has not been applied. Moreover I have 
taken the strong coupling to be equal one, while in reality it should be less 
than this value, and decreasing with increasing 

Q2. In Fig. KM I have 
plotted the relative uncertainty for the four fermion twist four correlation 
functions for the kinematics proposed in Ref. [9], and also in this case the 
uncertainty is, at the most, of the order of 0.1%. 

According to my model it seems that the normal belief that the four 
fermion HT correlation functions should be negligible is confirmed. 
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Figure 4.29: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions computed with the MMITBM set 1 for the kinematics 
proposed in Ref. [8] . 
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Figure 4.30: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions computed with the MMITBM set 1 for the kinematics 
proposed in Ref. [8]. 
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Figure 4.31: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions computed with the MMITBM set 2 for the kinematics 
proposed in Ref. [8] . 
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Figure 4.32: Relative uncertainty in the asymmetry due to the HT two 
fermion contributions computed with the MMITBM set 2 for the kinematics 
proposed in Ref. [8]. 
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Figure 4.33: Relative uncertainty in the asymmetry due to the HT four 
fermion contributions computed with the MITBM for the kinematics pro- 
posed in Ref. [9] . 




Figure 4.34: Relative uncertainty in the asymmetry due to the HT four 
fermion contributions computed with the for the kinematics proposed in 
Ref. [8]. 
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Chapter 5 



Bremsstrahlung 
Contribution to the 
Resonance Electroproduction 
Cross Section and Aju 

In this chapter I present a complete calculation of the bremsstrahlung con- 
tribution to the cross section and the PV asymmetry in the process 
e~ P ^ N* e~ , where N* represents a resonance. This process could be 
a major source of background in the experiment E-158 [35], recently com- 
pleted at the Stanford Linear Accelerator Center (SLAC). In this experiment 
the PV asymmetry Afi^i in the M0ller scattering process e~ e~ e~ e~ is 
measured, and the weak minxing angle sin^ Oyy is extracted. Since the elec- 
tron beam scatters on atomic electrons in a hydrogen target, there are many 
backgrounds that need to be considered due to the scattering of the elec- 
trons on the protons. Most of the backgrounds have been estimated in the 
proposal [35], but from preliminary results it was noticed that the inelas- 
tic background accounted for 40% of the total asymmetry measured. The 
goal of the experiment was to measure the M0ller asymmetry to a 8% level, 
which requires a knowledge of the inelastic background to a 20% of accu- 
racy or better. I study one of the possible inelastic background, namely 
the bremsstrahlung contribution to the A^l in the e~ P — > N* e~ process. 
In order to do so, I resort to a simple model in which only two resonances 
(the A(1232) and the A^(1520)) contribute. I investigate if the theoretical 
uncertainties on the resonances form factors, might produce an uncertainty 
on the asymmetry that is of the order of 8%. 
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Figure 5.1: Neutral current amplitudes contributing to the leading order to 
the M0ner cross section and Aji^. 

The chapter is divided as follow: in Sec. 15.11 1 describe the experiment E- 
158, while in Sec. 15.21 1 present the calculation for the polarized cross section 
for the bremsstrahlung contribution to the proton electro-excitation. In Sec. 
15.31 1 focus on the computation of the one loop corrections to the process 

P ^ e~ N*, which are necessary to make the bremsstrahlung process 
infrared finite. Finally, in Sec. 15.41 1 present the numerical results of the 
calculation and in Sec. 15.51 the conclusions. 



5.1 The E-158 Experiment. 

The E-158 experiment [35] performed at SLAC, is part of an extensive 
program of high precision experiments, with the objective of searching for 
physics beyond the SM. Experiments such as the LEP I and II at CERN, 
and the SLD at SLAC, have measured neutral current observables with an 
extremely high precision at the Zq resonance energy region [7l]. The E-158 
experiment tried to provide an high precision measurement of an electroweak 
observable, such as Ajil, away from the Zq pole, in order to explore a dif- 
ferent kinematical region in which to look for physics beyond the SM. 

At tree level, the processes contributing to Arl are depicted in Fig. 15. H 
and the result is [72l [35] 

Gf IGsin^e 

Arl = -mE— — — TFivi9^^ (^"^ 

V27ra (3 + cos^ 6)^ 

where Gp = 1.16639(1) x 10"^ GeV"^, m and E are the electron mass and 
incoming energy, respectively, is the scattering angle in the center of mass 
frame (CMF), and gee is the pseudo-scalar weak neutral current coupling. 
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Figure 5.2: Dominant contributions of the electroweak radiative correction 
to the polarized asymmetry Arl. 

that at tree level reads 

gee = -^9vgA = \-^™^ ■ (5.2) 

Radiative corrections reduce the magnitude of the tree level value of g^e up 
to 40% |72| . Some of the processes contributing to the radiative corrections 
are shown in Fig. 15.21 In the experiment, electrons with scattering angle in 
the CMF such that —0.5 < cosG < are detected. For an electron beam 
of 50 GeV this translates in the detection of electrons with energy between 
12.5 and 25 GeV. All the events are collected in a cylindrical detector, which 
is made up of concentric rings. The M0ller detector, in which most of the 
M0ller events are detected, consists in three rings: the IN, the MID and 
the OUT ring. Concentric to the M0ller detector there is the DEAD ring, 
whose only goal is to separate the M0ller ring from the outermost EP ring. 
In this last ring, most of the elastic and inelastic events, from the scattering 
of the electrons of the protons, are collected. 

Because the electrons in the target are atomic electrons in hydrogen 
atoms, the presence of the protons is source of a number of backgrounds 
which need to be considered. Most of them were estimated in the exper- 
iment proposal, but preliminary results from the experiment showed that 
the inelastic background was as big as 40% of the total asymmetry. Since 
the stated goal of the experiment was the measurement of Am to a 8% 
level of accuracy, the inelastic backgrounds need to be known to a 20% level 
or better. It was speculated that the bremsstrahlung process, in which a 
hard photon is emitted after the scattering with the protons in the hydro- 
gen atoms, could be a major source of the inelastic backgrounds. I divide 
this background in two major contributions: (i) the bremsstrahlung process 
in which the electron emits an hard real photon either before or after the 
scattering with the proton which gets excited into one of its resonances (Fig. 
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15. 3p . and (ii) the bremsstrahlung process in which the scattering with the 
proton is deeply inelastic, breaking up the target in a multitude of particles. 
In both cases the main concern is that, because of the fact that the electron 
emits an hard photon, it could lose sufficient energy to end up in the M0ller 
instead of in the EP detector as expected. Even though the deep inelastic 
and the electro-excitation cross sections are much smaller than the M0ller 
one, the asymmetry for these two processes is one order of magnitude or 
more larger than the M0ller one, potentially making the total asymmetries 
comparable. Another problem is that only for few resonances the electro- 
magnetic form factors are well known, and even for these, the corresponding 
the PV FFs are usually very uncertain for most of them. These two facts, 
the possible large value of the resonance asymmetry and the uncertainty 
in the resonance form factors, suggest that a complete numerical compu- 
tation of the bremsstrahlung contributions to the proton electro-excitation 
asymmetry be performed. However, a complete calculation is unrealistic. In 
this model I simplify the situation by dividing the events in two kinematical 
regions. I assume that all the events with energy between 30 and 50 GeV 
are detected in the EP ring, while the events with energy between 12.5 to 
30 GeV end up in the M0ller ring. To see how the uncertainty on the form 
factors might impact the measurement of Arl, I compute the contribution 
of two resonances to the bremsstrahlung asymmetry in the EP ring (in this 
energy range the M0ller contribution is negligible) , by assuming some sensi- 
ble value of the form factors. For these conditions I also compute the total 
asymmetry in the M0ller ring, i.e. the asymmetry given by the sum of the 
M0ller and bremsstrahlung process for the two resonances. This is be my 
reference value, which is used for comparisons. I then allow the form factors 
to change in a reasonable range, in such a way that Am in the EP remains 
the same, and then compute the new M0ller asymmetry, and compare it 
with my reference value. In practice, the idea is to see how much a mea- 
surement of the asymmetry in the EP might constrain the asymmetry in 
the M0ller ring. The value of Arl in the EP ring I kept constant might 
be considered as my "measurement". Of course, there could be different 
values of the form factors that might provide such a outcome. The issue I 
am trying to address here is if these different values, all compatible with the 
potential measurement in the EP, will result in a large uncertainty in Arl 
extracted in the M0ller ring. 
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Figure 5.3: Electromagnetic current amplitude for bremsstrahlung process. 




Figure 5.4: Neutral current amplitude for bremsstrahlung process. 



5.2 Cross Section and Arl for Proton Electro-Excitation. 



In this section I focus on the computation of the bremsstrahlung contribution 
in the proton electro-excitation to the cross section and the polarized asym- 
metry. I consider a simple model in which only two resonances contribute, 
the isovector A(1232) and the isoscalar A^(1520). The electromagnetic and 
the neutral current amplitudes, contributing to the process I am consider- 
ing, are depicted in Figs. 15.31 and 15. 4i Let us consider the case in which 
the incoming electrons move along the z axis and scatter in the (x, z) plane. 
The following kinematical quantities may be defined: 
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k = (^,0,0, Ve"^ - m^) 

incoming electron four momentum 
s = lb — (\/ E"^ — m^, 0, 0, E) incoming electron spin four vector 



m 



ki = {E',VE^ — m"^ sin 9, — 771,2 COS 0) outgoing electron 

four momentum 

^2 = ^o(l) sin a cos (p, sin a sin cf), cos a) real photon four momentum 

P = (M, 0,0,0) proton four momentum 

Pn* = (i-W*Oi -P/V1) -P/v*2i -P/v*3) resonance four momentum 

q = k — ki — k2 four momentum transferred 

Q2 = -q\ (5.3) 



With such definition the amplitude in Fig. 15.31 can be computed, and the 
result is 
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where e is the electron charge and £0(^2) is the photon polarization vec- 
tor and i labels th, si and S2 are the spins of the incoming and outgoing 
electrons, respectively, polarization state of the outgoing photon, ai is the 
proton spin and (T2 is the resonance spin. The differential cross section then 
becomes 



dn ~ Tgme ^ ^ 

S1S2 <Ticr2 
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{Ai^fME 



(5.5) 
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where I have introduced the hadronic and leptonic tensors W^"^ and Lem 
respectively defined as 

6 {M%, -{P + k-ki- k2f) = 



1 



A J 2PoiV*(27r)= 



and 
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where I have used the relations 



(^i+m)}, 
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Yus{p)us{p) = ^ + 
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2ki ■ k2 
(5.8) 

(5.9) 
(5.10) 



where in the last relation the ellipsis indicates terms which do not con- 
tribute to the L^^. The leptonic tensor, although tedious and lengthy, can 
be computed. The result is presented in appendix iDl 

Using the delta function to eliminate the integration over Eq. (j5.5p 
may be written as 



da 

dn 
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d^k2 W^T^em E' 
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(5.11) 



where I have defined 
dE' 



2{M + E{1 — cos 9) — ko{l — sin cos (?!) sin a — cos a cos 0)) , 

(5.12) 
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and where E' is now equal to 



E' 



-M%, + M2 + 2m? + 2ME - 2Mko - 2Eko{l - cos a) 



2{M + E{1 -cos( 



ko{l — sin 6* cos (/) sin Q — cos a cos 0)) 

(5.13) 



On general grounds, invoking Lorentz, space-inversion and time-reversal in- 
variance, and current conservation, the hadronic tensor can be casted in the 
following form: 



-9fiu + 



+ P. 



P^q 



P-q 



(5.14) 



and the structure functions (SFs) Wf^'' and VFJ™ contain information on 
the structure of the resonance. Let us now express the SFs in terms of the 
hadronic current form factors. For the two spin 3/2 resonances, I take the 
hadronic current matrix element to be of the following form |73l I74j 
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(5.15) 



For the A(1232) resonance the following assumptions describe the data cor- 
rectly [75]: 

C2 = 0, Cl = -^Cl. (5.16) 
From Eqs. ()5.6p and ()5.15p . and from the relation [76] 
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3 Ml, 



3Mn- 
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-fPN*+ Mn* 



2Mn* 

it is straight forward to compute the hadronic tensor and, after performing 
some Dirac algebra and equating the result to equation (j5.14p . I get, for the 
structure functions Wf"^ and W^^"^, the following results: 



(5.17) 



"6M^" ^3 



(5.18) 
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(5.19) 
(5.20) 



-2 - ^ --3 

and I have introduced the following kinematical factor 

a = (M + MN'f + Q"^ 
and 

c= {M - MN*f + . (5.21) 

In order to simply the problem, I also assume that the relations in Eq. ()5.16p 
are true also for the A^(1520) resonance. 

Next, I want to compute the polarized asymmetry defined as 
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where the subscripts R, L refer to the helicity of the incoming electrons. 
Notice that, when considering different scattering processes, the total asym- 
metry is defined as 
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(5.23) 



where 
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(5.24) 



and A^j^ and ak are the asymmetry and the cross section for each single 
process, respectively. In this case I consider three processes: the M0ller 
scattering, and the bremsstrahlung processes of the two resonances. 

Since the electromagnetic contribution to the cross section is insensitive 
to the polarization of the beam, taking the difference in the asymmetry 
isolates the PV part of the process, which changes sign upon flipping the 
helicity of the electrons. In this case the amplitude of the neutral current 
process, in the limit ^ Q^, is 
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and the numerator of Eq. 15.221 becomes 
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(5.26) 



where Gp = 1.6639(1) x 10 ^ GeV ^ is the Fermi constant. 

In complete analogy with what I have done before, I have defined the 
interference leptonic tensor as 



rnc 



^^Us^iki) 



(7"^^2/cf) _^ 
e„, — , T{gv + 9a15) 



S2 I 



2ki ■ k2 



e ,(7°^2-2A:°) 



X Us^{k) 
= -Tr 



2k ■ k2 

■ , (7^^^2fcf) 
7 — TTT, — Z - ^13 



2ki • A:2 
(7"fc^ + 2fcf ) 



tisi(A;) 

(7/^^2-2fc^) 
2k ■ k2 



Us2{ki 



2ki ■ k2 



i^'^J^i^ 2kf) 



X (^ + m)- 



7' 



2A;i • /C2 
(^i+m)} . 

(5.27) 



Notice that in this case there is no sum over the incoming electron spin, for 
this reason Eq. (|5.9p now reads 

Us{p)us{p) = {j> + m)i^^ . (5.28) 
The full expression of the interference leptonic tensor can also be found in 
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appendix [D1 The hadronic tensor for the interference process is 

{2TTf5^ (Ptv* -{P + k-ki- k2)) 

{p.AjrmN*){N*\JTmp..), 



1 



4 7 2PoN'{2Tif 



(5.29) 



and similarly with what has been done before for the electromagnetic case, 
its most general form compatible Lorentz and time-reversal invariance, and 
current conservation, is 



P-q 



P„ 



P. 



P-q 



(5.30) 



If I now assume a form for the matrix element of the hadronic neutral current 
as follows [H] 



u^{Pn'\ 



+u^{Pn') 



jZ 

^Pn'^Mpu - g\pgpu)q^ + C^A9\^l 
the hadronic form factors could be computed, and what I obtain is 



u{P) 
u{P) (5.31) 



wr 



3MMn* ^2 



(5.32) 

(5.33) 
(5.34) 



where also in this case, for simplicity, I made the approximation C^y 
~ * ^iv electromagnetic case, and b is defined as 



6 = (M + Mn'){M - Mn*) + Q2 . 



(5.35) 
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Following Ref. [2], the hadronic currents can be decomposed in terms of 
their isospin content (see also Sec. 12. 2p 

■/r(0) = -^J='(0) + -^r°(0)' (5-36) 

and 

^r(O) = er'^r'(O) + ^^v^'J^'i^) + M^) > (5-37) 

where 

^^=^ = 2(1-2 sin^^v^) , ^^=0 = -Asin20i^, (5.38) 

and the superscript T = 1 and T = stands for isovector and isoscalar, 
respectively. The vector ^^(0) represents the axial part of the hadronic 
neutral current, which I am not interested in decomposing. Because of the 
fact that the A(1232) resonance is an isospin | particle, while the proton 
has isospin ^, the isoscalar component of the currents does not contribute 
to the matrix elements (P/v*| J^"*'"'^(0)|P). As a consequence, I have that 
the interference vector form factor is related to the electromagnetic one by 
the relation 

Civ = Cv='C2. (5.39) 

In the case of the A^(1520), since it is an isospin ^ particle, both components 
of the hadronic current contribute, so that I have 

= cf=^ + C73^^=° (5.40) 



Civ = i^='cf=' + V-Hl='cf=' . (5.41) 

In order to compute the bremsstrahlung cross section and asymmetry, ac- 
cording to Eqs. ()5.11|) and (j5.26p . I would need now to compute the two 
tensor contractions L^^J^Wem and L^'^Wnc . The results can be found in the 
appendix |Dl 

5.3 Infrared Cancellation 

In this section I show how the bremsstrahlung cross section and asymmetry 
become infrared finite once I had the one loop vertex corrections shown in 
Figs. 15.51 and 15.61 to the electro-excitation process. As will be shown in a 
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Figure 5.5: One loop correction to the EM cross section. 



moment, the integral over the photon energy in the bremsstrahlung cross 
section and Arl, diverges if the lower limit is taken to be zero. In order to 
keep it finite, I introduce a lower limit cut off. Of course, since the cross 
section and Arl are physical quantities, they need to be finite and cut off 
independent. Adding the one loop contribution will indeed make them finite 
and cut off independent. The reason why I compute to soft bremsstrahlung 
contribution is that this process, even though does not contribute to the 
asymmetry in the M0ller ring, it still needs to be considered when computing 
the total asymmetry in the EP ring. 

First of all let us consider the zero photon energy limit in the square 
amphtudes W^^L^^;; and W^^^^L^c- The lengthy relations in Eqs. (fD:3]l and 
()D.4p become 



lim VF^T^^m 



[m 



q2) ^ ^|™(4£;2 _ q2 



2k -ki 



m 



m 



iki-k2){k-k2) {k-k2Y {ki-k2f 



(5.42) 



and 



^2E + E'^ f 2k-ki m2 



M ^ / V(^i-^2)(A:-A:2) (fc-fc2)2 {ki-k2y 

(5.43) 



In order to obtain the cross section and Arl, I need to perform the following 
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Figure 5.6: One loop correction to the polarized cross section ai — ctr. 



integral 



/ 



ko \{ki-k2){k-k2) {k-k2Y {ki-h 



2) 

dko f I 2k-ki vn? m? 

ail 



ko J \{ki-k2){k-k2) {k-k2f {ki-k2f 



(5.44) 



where k2 = kl^/kQ, ^ is the lower limit cut off, and Ei is an upper limit cut 
off for which the soft photon approximation is still valid. The integration 
over the solid angle of the last two terms is straight forward. Taking the 
four-vector A;^ (or k^) along the z axis, I have 



2 f 2 !• 1 

m I m I m 



dn — . — = / dn — . — = / dn 



{ki-k2Y J {k-k2Y J {E-\k\cose) 



\2 



2tt I dy 



rr? ^'Krr? /I 1 



2 

m 



1 ''(^-|k|y)2 |k| V^-|k| 

(5.45) 



47r— ^ = 47r , 
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while for the first term I have 



2k-ki f 2k-ki 



{ki-k2){k-k2) Jo 



(^xki-k2 + (1 — x)k-k2 



2 



J dx J dQ— — = Att [ dx- ^ 



(k2-{-xq + k)) Jo {-xq + kY 

= A.f'dx^^^^ (5 46) 

where I introduced the Feynman relation 

dx 



AB Jq {Ax + B{1 - x))2 ' 



(5.47) 



Putting all together and performing the integral over the photon energy, I 
get for Eqs. ()5.1ip and (j5.26p in the soft photon approximation 



- I = ^ log — X 



dVt J 2ttME[M + £^(1 - cos d)\Q^ ^ n 

^ dx{Q'^ + 2m2) 
{Q^x{l-x)+m^) ~ 



(5.48) 



(daLRY'"" ^ »'E'Gf El 

\ 6n J 8V2tt^ M E[M + E {I - cos e)]Q'^ 



X 



'g-^[Wr{^E'-Q')+^WrQ'] + 2g'vQ'^-^^^W, 



M 

^ dx{Q^ + 2m2) 
(Q2x(l — x) + m?] 



(5.49) 



where I did not consider finite terms in the limit in which /i goes to zer c0. As 
can be seen, in such a limit, both expressions are logarithmically divergent. 
In order to obtain a finite expression I need to add the one loop corrections. 
This contribution to the cross section is lengthy and cumbersome to obtain. 



^In the actual calculation, the sum of the bremsstrahlung and the loop contribution 
has been performed numerically, so the finite terms left out in Eqs. (|5.48|) and (|5.49|l are 
included 
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so I delineate only the main steps, referring the reader to standard books. 
My derivation follow closely the one that is found in [77]. Considering the 
process represented in Fig. 15.51 the cross section is 



.L^^WZ, (5.50) 



dnj ME[M + E{1- cos e)]Q 

where I have introduced a new leptonic tensor relative to the vertex correc- 
tion 



L 



d'^p 



J (2^)4 (p2_^2)((p + ^)2_^2)((fc_p)2_^2) 

(5.51) 

while the hadronic tensor is the same as for the bremsstrahlung case. Notice 
that I have introduced a fictitious mass // in the photon propagator to make 
the integral finite and which, at the end of the day, cancels the cut off 
dependence in the bremsstrahlung cross section. Using standard procedures 
in performing the integral and keeping only terms which are singular in the 
limit fi goes to zero, equation (j5.50p gives 



dn (47r)2M£;[M + E(l- cos 6l)]Q4 



F{y,z,Q^ 



1 pl-z 

dz / dy 

Jo [^„^2(l_2)2_^g2y(l_^_y)_^^22j_ 

(5.52) 

where I have defined the function F{y, z, Q^) as 

F{y,z,Q^) = 16Wr{[{z-^)z + M^ri'' + m^Q^[y^ + {z-l)y 

+ z^ + z-l]-Q\y-l){y + z)) 

- AWr[^{m\z - A)z + 1] + Q\y - l){y + z)}E^ 

- Q^{-m\z{z + 2)-l]+Q\y-l){y + z)}] . (5.53) 



The integral in Eq. (j5.52p in the limit ^ goes to zero, becomes divergent in 
the region z = 1 and y = 0. Since I am concerned in the divergent pieces of 
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the cross section, I take such a hmit where possible. The function F{y, z, Q"^) 
becomes then 



so that equation (|5.52p yields 



(5.54) 



IttMEIM + E{1 - cos6l)]Q4 



dz 



dy 



(2m2 + Q2 



m?{l — z)2 + Q2y(l — z — y) + fi'^z 



5.55) 



Beside being infrared divergent, the one loop cross section is also ultraviolet 
divergent. Making it finite is part of the renormalization program. In this 
case I eliminate the divergent part by doing the following substitution in the 
integral in Eq. (j5.55|) (see [77] chapters 6 and 7 for more details) 



(2m2 + g2) 









C ^ dy 




Jo 








C ^ dy 




10 



m2(l — z)2 + Q'^y{l — z — y) + /i^z ) _ 

(2m2 + Q2) 

r?T,2(l — z)2 + Q2y(l — z — y) + fj.'^z 



2m^ 



m2(l — z)2 + ^2-- 



(5.56) 



where basically I have subtracted the same quantity calculated at = 
(on-shell renormalization scheme). Performing the change of variable z = 
1 — t and y = tx the first integral of the previous relation becomes 



dz 



dy 



(2m2 + Q2 



771,2(1 — 2;)2 + Q^y{l — z — y) + p?z 



dx I dt- 



dx 



t'^[m? + Q'^x{l- x)) + ^i'^{l-t) 
m? + Q'^ ( iv? + Q'^x{l - x) 



2[m2 +Q2^(1 _^)] 



log 



(5.57) 
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The first integral is straight forward and the result is 



2m^ 



1 i-l-z 

dz / dy-r- 



m 
log^ , 



(5.58) 



so that the loop contribution yields 



do_'°''P 

dn 
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'((4^2 _ q2)^2 - 4(m2 - Q'^)Wi^E' 
2ttME[M + E{1 - cos6')]Q4 
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log- 



(5.59) 



Now adding equations (|5.48p and (|5.59p I get 



do_' 
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(5.60) 



which is independent on the cut off fi as expected. At the same way the 
one loop corrections to the polarized cross section for the process in Fig. 15.6 
have been calculated, leading to 
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8V2Tr'^ME[M + E{1 - cos e)]Q'^ 



g'x[W2{4.E^ - Q2) + 4t^iQ2] + 2g,Q 



2{E + E') 



dx 



oi El 
21og — 

m 



(2m2 + Q2 



m? + Q'^x{\ — x) 



log 



M 
El 



yjm? + Q'^x{l 



(5.61) 
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5.4 Numerical Studies 



In this section I present the numerical results of my two resonance model. 
The idea, as I mentioned before, is to compute the asymmetry in the EP 
region, which in my case is the range of energy of the ongoing electrons that 
goes from 30 to 50 GeV. All cross sections and Arl^s are integrated over 
the angular resolution, which corresponds to a range of scattering angle in 
the lab frame between 4.5 and 8.0 mrd. In order to simplify the problem, 
I assume the same dependence in all the form factors. Specifically, I 
assume a dipole dependence of the type 

where the parameter a is taken to be the same for all the form factors. For 
the delta resonance electromagnetic form factor, I have 

and, for the interference one, according to equation (j5.39p . I have 

CiyiQ')=C^='C]{Q'). (5.64) 
For the A'^(1520) resonance I write 



(5.65) 



and, according to equation (j5.4ip 

CiviQ') = ev='cf=\Q') + V2,il='cf=\Q') . (5.66) 

In order to fix the initial parameters I have considered the cross sections 
e~ P — > vr"*" N and e~ P ^ tt^ P. The MAID group provides a computational 
tool to compute the aforementioned cross sections, allowing one to isolate 
the contribution arising from each single resonance (it can be found in their 



web site at http://www.kph.uni-mainz.de/MAID/) [78]. I have then best 



fit the parameters to the cross sections, and what I get is 

a = .54 C3 = 1.3 C31 = .7 C30 = ^ . (5.67) 
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Another approximation I have made, is to consider a zero width resonance, 
which considerably reduces the computational time. At this point I consider 
the total asymmetry in the EP ring obtained for this initial set of values of 
the form factors, and consider it as a possible result from the measurement. 
I also compute the total asymmetry in Arl in the M0ller ring, and consider 
it as a reference value . What I want to see is if this measurement in the 
EP ring could actually constrain the measurement in the M0ller ring. The 
same value of the asymmetry in the EP ring could have been obtained also 
for some different value of the form factors. I change the value of C3 and C31 
and treat them as independent variables, and find the value of C30 such that 
the asymmetry in the EP ring does not change. I then compute the total 
asymmetry in the M0ller ring for these new values of the form factors and 
compare it to the reference value A^j^ to check if the variation in the form 
factor induced a large difference in Am computed in the M0ller ring. The 
contribution to A^ arising from the PV SF W^'^ are not explicitly included, 
since it turned out to be negligible (remember that this term always appears 
multiplied by the vector coupling constant gy ~ .075). The total asymmetry 
in the EP ring is 



^EP 



1 



(^^ - + {< - ) (5.68) 



in which I have neglected the contribution from the M0ller process since 
it is much smaller than the bremsstrahlung. Considering the values of the 
parameters in Eq. (j5.67p I get the following initial values for the EP ring 



= 3.18 X 10"^GeV-2 


(5.69) 


-a^ = 2.61 x 10~^GeV-2 


(5.70) 


= 9.38 X lO-^GeV-2 


(5.71) 


- a;^ = 6.35 x 10-i°GeV-2 , 


(5.72) 



giving a total asymmetry of 

A^^ = -7.88 X 10"^ . (5.73) 

In the M0ller ring I consider both processes, the bremsstrahlung and the 
M0ller, so that the total asymmetry can be written as 



(5.74) 
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where the M0ller differential cross section and polarized cross sections are 
[35] 

da^ _ (3 + cos^9)2 

dn ~ 2mE sin^ G ^ ' ' 

and 

daff da^ _ Gpa gee 

dff" V27rsin2e' ^ ' 

and where I used for the coupling gee its tree level value corrected by a 
factor of 0.4 due to radiative corrections [72]. The total cross sections for 
the M0ller process integrated over the angular resolution are 



M 



4.18 X 10~2Gev-2 (5.77) 



-a^ = 6.23 x IQ-^Gev-^ , (5.73) 
while, for the two resonances I get 

= 2.2 X 10"^GeV~2 (5.79) 

(j'^-a^ = 9.77 X 10"^^GeV~2 (5.80) 

= 6.72 X lO^^GeV"^ (5.81) 

-a^ = 2.44 x 10~"GeV-2 (5.82) 
so that my reference value for the total asymmetry is 

^^^^0 = -1.52 X 10"^ (5.83) 

In Fig. 15.71 1 have plotted the quantity 

r 4 aMTOT _ aMO 

^ ^ A^l^"^ + A'^jII ' (5.84) 

where A^^'^'^ is the value of the total asymmetry for which I allow the 
parameters C3 and C31 to change up to 20% of their initial value. Notice 
that the factor of two in the definition of the relative error, comes from 
taking the average of two asymmetries. It seems quite clear that in this 
simple model, even a very poor knowledge of the resonance form factors 
(which is not the case for the ones I used) does not impact severely the 
extraction of the total asymmetry, since allowing a very "generous" mistake 
on the form factors of 20%, will only reflect in a less than 1% error on the 
asymmetry. Nevertheless, including less known resonances to the model, 
might result in a theoretical uncertainty of the order of the experimental 
error. 
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Figure 5.7: Relative uncertainty in the total asymmetry defined as in Eq. 
(|5.84|) as a function of C3 and C31 with initial values C3 = 1.3, C31 = .7 and 
C30 = I/IOC31. The parameters C3 and C31 have been changed up to 20% of 
their initial value, while C30 has been changed accordingly, in order to keep 
the total asymmetry in the EP ring constant. 
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5.5 Comments and Conclusions 



In this chapter I have investigated the contribution of the bremsstrahlung 
radiation to the proton electro-excitation process, by using a simple two 
resonance model. I have studied how the possible uncertainties on the reso- 
nances form factors might be source of error in the measurement of Arl in 
the E-158 experiment. As it turned out, even knowing the form factors to 
20% level, only affects A^il with an error of the order of less then one per- 
cent, which is much smaller than the prescribed precision of the experiment, 
which is 8%. I tried to see what would be the error on the form factors 
that produces a 8% uncertainty in the asymmetry. What I found is that the 
error should be of the order of 150%. Even though such uncertainty might 
seem large for two resonances such as the Delta{1232) and the A^(1520), 
there are still many resonances that contribute to the bremsstrahlung pro- 
cess, whose FFs are not very well known. A more complete study, including 
other resonances, seems necessary in order to completely exclude the possi- 
bility that the bremsstrahlung contribution to the proton electro-excitation 
might introduce an uncertainty of the order of the prescribed precision of the 
experiment. I also changed the dependence by varying the a coefficient 
in Eq. (j5.62p between 0.4 and 0.7, but the result did not change appreciably. 
On the other hand, in my simple model, I have made some assumptions. I 
assumed that the dipole dependence for the form factors was the same for 
both resonances. Moreover, for the isoscalar resonance I took the depen- 
dence to be the same for the isovector and the isoscalar form factors, which 
is not necessarily true. It would be also interesting to be able to include 
the continuous DIS region in the calculation. The main contribution to the 
bremsstrahlung process arises from the kinematical region in which the pho- 
ton is emitted parallel either to the incoming or the outgoing electron. In 
this region the it is no more than 0.1 GeV^, so, in order to compute the 
DIS contribution, one would need to know the proton SFs at a very small 
Q^, but this kind of knowledge is lacking at the moment. In conclusion, it 
seems that in this two resonance model, the bremsstrahlung contribution to 
the proton electro-excitation should not vitiate the extraction on the M0ller 
asymmetry at a 8% level. However, I also showed that a 150% uncertainty 
in the resonances FFs, might produce an uncertainty in the asymmetry in 
the M0ller ring as big as the prescribed experimental precision. This result 
suggests that a more detailed study be done, including a larger number of 
resonances contributing to the bremsstrahlung process. 
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Chapter 6 



Electroweak Radiative 
Corrections to the 
Parity- Violating 
Electroexcitation of the A 



In this chapter I am going to present the results of the work I have done 
in collaboration with Ramsey- Musolf, Zhu, Holstein and Maekawa [36]. In 
this work, we analyze the degree to which parity violating electro-excitation 
of the resonance A(1232) may be used to extract the weak neutral axial 
vector transition form factors. My contribution consisted in double check- 
ing the computation of the chiral loops and the l/Mjy corrections to the 
parity violating asymmetry arising from the anapole, the Siegert, and the 
d-wave parity violating lagrangian. In the next sections I will summarize 
the calculation and present the main results of this work. More details can 
be found in ref. 1361 . 



6.1 Introduction 

The electroweak form factors associated with the excitation of the A(1232) 
resonance are of considerable interest to hadron structure physicists. In 
the large Nc limit, the {N, A) form a degenerate multiplet under spin-flavor 
SU(4) symmetry [79], and one expects the structure of the lowest-lying spin- 
1/2 and spin-3/2 qqq states to be closely related. The electroweak transition 
form factors may provide important insights into this relationship and shed 
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light on QCD-inspired models of the lowest lying baryons. These form fac- 
tors describe ^ A matrix elements of the vector and axial vector currents 
[801 ED [82]: 



(6T) 

< A+ip')\Al\N >= A+'^(pO{[^7V^/^](<?A5M^-g.5A;.)+C5V+^9M9.}«b) 

(6.2) 

where the baryon spinors are defined in the usual way. The form factors 
and are the N —>■ A analogues of the nucleon's electroweak form factors 
Fi and Ga- At present, there exist considerable data on the vector current 
transition form factors = 3 — 6) obtained with electromagnetic probes. 

A comparison with theoretical predictions points to significant disagreement 
(see ref. [83| for a tabulation of theoretical predictions). For example, lattice 
QCD calculations of the magnetic transition form factor yield a value ~ 30% 
smaller than obtained from experiment [84], and constituent quark models 
based on spin-flavor SU(6) symmetry similarly underpredict the data|85|. 
One hopes that additional input, in tandem with theoretical progress, will 
help identify the origin of these discrepancies. 

The situation involving the axial vector transition form factors Cf- {i = 
3 — 6) is less clear than in the vector case, since existing data - obtained 
from charged current experiments - have considerably larger uncertainties 
than for the vector current channel. While QCD-inspired models tend to 
underpredict the central value for the axial matrix elements by ~ 30% as 
they do for the vector form factors, additional and more precise experimental 
information is needed in order to make the test of theory significant. To that 
end, an extraction of the axial vector — > A matrix element using parity- 
violating electron scattering (PVES) is planned at the Jefferson Laboratory 
[37j . The goal of this measurement is to perform a < 25% determination 
for \q^\ in the range of 0.1 — 0.6 (GeV/c)^. If successful, this experiment 
would considerably sharpen the present state of experimental knowledge of 
the axial vector transition amplitude. 

Here we examine the interpretation of the prospective measurement. In a 
previous work [83] , the impact of non-resonant backgrounds was studied and 
found not to present a serious impediment to the extraction of the Cf. Here, 
we compute the electroweak radiative corrections, which arise from C(aGi?) 
contributions to the PV axial transition amplitude. We correspondingly 
characterize the relative importance of the corrections by discussing the 
ratio of the higher-order to tree-level amplitudes. This ratio is nominally 
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0{a), so that one might naively justify neglecting radiative corrections when 
interpreting a 25% determination of the axial term. However, previous work 
on the axial vector radiative corrections to PV elastic electron-proton 
scattering suggests that the relative importance of such corrections can be 
both unexpectedly large and theoretically uncertain [381 [39l HD] • Moreover, 
results obtained by the SAMPLE collaboration [86] suggest that may 
be substantially larger than given by the best theoretical estimate |38|. The 
origin of this apparent enhancement is presently not understood. Were 
similar uncertainties to occur for PV electroexcitation of the A, the task of 
extracting the desired axial transition form factors from the PV asymmetry 
would become considerably more complicated than assumed in the original 
incarnation of the experimental proposal. 

In studying the axial vector radiative corrections, it is important to dis- 
tinguish two classes of contributions. The first involves electroweak radiative 
corrections to the elementary V{e) x A(q) amplitudes, where q is any one of 
the quarks in the hadron and V (A) denotes a vector (axial vector) current. 
These terms, referred to henceforth as "one-quark" radiative corrections, are 
calculable in the Standard Model. For elastic scattering from the proton, 
they contain little theoretical uncertainty apart from the gentle variation 
with Higgs mass, long-distance QCD effects involving light-quark loops in 
the Z — 7 mixing tensor, and SU(3)-breaking effects in octet axial vector 
matrix elements {p\A^^'^^\p). Such one-quark contributions to R^ and R^ 
can be large, due to the absence in loop terms of the small (1 — 4sin^ 9yy) 
factor appearing in the tree level V{e) coupling and the presence of large 
logarithms of the type ln{mq/ Mz). 

The second class of radiative corrections, which we refer to as "many- 
quark" corrections, involve weak interactions among quarks in the hadron. 
In refs. [381 [39l I40j . the many-quark corrections were shown to generate 
considerable theoretical uncertainty in the PV, axial vector ep amplitude. 
A particularly important subset of these effects are associated with the nu- 
cleon anapole moment (AM), which constitutes the leading-order, PV jNN 
coupling. The result of the SAMPLE measurements, which combine PV 
elastic ep and quasielastic ed scattering to isolate the isovector, axial vec- 
tor ep amplitude, implies that the one-quark/Standard Model plus many- 
quark/anapole contributions significantly underpredict the observed value 
of i?^. 

In what follows, we compute the analogous radiative corrections R^ 
for the axial — > A electroexcitation amplitude. In principle, as in the 
elastic case, the one-quark corrections are determined completely by the 
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Standard Model, although long-distance QCD effects - which are finessed 
for the ep channel using SU(3) symmetry plus nucleon and hyperon /3-decay 
data - are not controlled in the same manner for the N ^ A transition. 
We make no attempt to estimate the size of such effects here. Instead, we 
focus on the many-quark contributions which, as in the elastic case, can 
be systematically organized using chiral perturbation theory (xPT). We 
compute these corrections through 0{p^). We find: 

(i) As in the case of i?^, the correction is both substantial and theo- 
retically uncertain. Thus, a proper interpretation of the PVES ^ A 
measurement must take into account O^aGp) effects. 

(ii) In contrast to the elastic PV asymmetry, the — > A asymmetry 
does not vanish at = 0. This result follows from the presence of 
an O(aGir) contribution - having no analog in the elastic channel - 
generated by a new PV 7A^A electric dipole coupling d^- Specifically, 
we show below that 

3 X 

where Aj^i^^q"^) is the PV asymmetry on the A resonance, A^ = 4:7rF.„- ~ 
1 GeV is the scale of chiral symmetry breaking, ~ 2 is the domi- 
nant A^ — > A vector transition form factor, d/\ is a low-energy constant 
whose scale is set by hadronic weak interactions, and the + • • • denote 
non-resonant, higher order chiral, and 1/Afjv corrections. 

(iii) The experimental observation of surprisingly large SU(3)-violating con- 
tributions to hyperon radiative decays suggests that the effect of c^a 
could be significantly enhanced over its "natural" scale, yielding an 
A^ ^ A asymmetry ~ 10~^ or larger at the photon poinl|j. 

(iv) The presence of the PV coupling implies that the Q'^-dependence 
of the axial vector transition amplitude entering PV electroexcitation 
of the A could differ significantly from the g^-dependence of the corre- 
sponding amplitude probed with neutral current neutrino excitation of 
the A. As we demonstrate below, it may be possible to separate the d^ 
contribution from other effects by exploiting the unique g^-dependence 



^For a PV photoproduction asymmetry of this magnitude, a measurement using po- 
larized photons at Jefferson Lab would be an interesting - and potentially feasible [87] - 
possibility. An analysis of the real 7 asymmetry appears in a separate communication 
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associated with this new term. We ihustrate this possibihty by con- 
sidering a low- 1 (/^ I, forward angle asymmetry measurement. 

(v) An experimental separation of the cIa contribution from the remaining 
terms in the axial vector response would be of interest from at least two 
standpoints. First, it would provide a unique window - in the A5 = 
sector - on the dynamics underlying the poorly understood PV AS = 1 
radiative and nonleptonic decays. Second, it would help to remove a 
significant source of theoretical uncertainty in the interpretation of the 

^ A asymmetry, thereby allowing one to extract the N ^ A axial 
vector form factors with less ambiguity. 

(vi) A comparison of PV electroexcitation of the A with more precise, 
prospective neutrino excitation measurements would be particularly 
interesting, as inelastic neutrino scattering is insensitive to the large 
7-exchange effects arising at O^aGp) which contribute to PV electron 
scattering [39lH0] . 

While the remainder of the paper is devoted to a detailed discussion 
of these points, several aspects deserve further comment here. First, the 
origin of the nonvanishing ^^^(g^ = 0) in Eq. (|6.3p is readily understood 
in terms of Siegert's theorem [89^ I90j . familiar in nonrelativistic nuclear 
physics. For electron scattering processes such as shown in Fig. 16.11 the 
leading PV 7-hadron coupling (Fig. 16. Id ) corresponds to matrix elements of 
the transverse electric multipole operator TjLj^^, and according to Siegert's 
Theorem, matrix elements of this operator can be written in the formn 

{f\ff=ix\i) = -^^{f\ I d'xxY,x{n)p{x)\i) + 0{q^) , (6.4) 

where the u = Ef — Ei. The leading component in Eq. (j6.4p is q^- 
independent and proportional to u) times the electric dipole matrix element. 
Up to overall numerical factors, this El matrix element is simply dj\/A-^. It 
does not contribute to PV elastic electron scattering, for which to = 0. The 
remaining terms of 0{q^) and higher contain matrix elements of the anapole 
operator |9H HQ] , which generally do not vanish for either elastic or inelastic 
scattering. When {f\Tj^-^)^\i) is inserted into the full electron scattering am- 
plitude, the 1/q^ from the photon propagator cancels the leading (f' from the 
anapole term, yielding a g^-independent contact interaction. In contrast, for 
inelastic processes such as electroexcitation of the A, w = rriA — w-jv does not 

■^We adopt the "extended" version of Siegert's theorem derived in ref. |90) . 
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Figure 6.1: Feynman diagrams describing resonant pion electroproduction. 
The dark circle indicates a parity violating coupling. Fig. 16. Id gives tran- 
sition anapole and Siegert's term contributions. Fig. 16.1b leads to the PV 
d-wave vrA^A contribution. 
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vanish, and the dipole matrix element in Eq. (j6.4p generates a contribution 
to the PV scattering amphtude Mpy behaving as for low-|g^|. Since 
the parity-conserving (PC) amphtude Mpc - whose interference with Mpy 
gives rise to A^p - also goes as the inelastic asymmetry does not van- 
ish at the photon point. Henceforth, we refer to the dipole contribution to 

the asymmetry as A^^p , and the corresponding 0{a) correction to the 

0{Gp) Z'^-exchange, axial vector neutral current amplitude as i?^^^^^'"^. 

We note that the importance of Aj^^^^^^ - relative to the anapole and Z^- 
exchange contributions to the asymmetry - increases as one approaches the 
photon point, since the latter vanish for = 0. 

It is straightforward to recast the foregoing discussion in a covariant 
framework using effective chiral Lagrangians. The dipole term in Eq. (j6.4p 
corresponds to the operator 



^Siegert ^ ^^a+tapF^^ + H.c. (6.5) 



while the transition anapole contribution arises from the effective interaction 
^anapole ^ ^^+^^^^Am + H.c. . (6.6) 



A2 

X 



6.2 Electroexcitation: general features 

The amplitudes relevant to PV electroexcitation of the A are shown in Fig. 
16.11 The asymmetry arises from the interference of the PC amplitude of Fig. 
la with the PV amplitudes of Figs. I6.1b -e. In Fig. I6.1b -d. the shaded circle 
denotes an axial gauge boson (V)-fermion (f) coupling, while the remaining 
V-f couplings are vector-like. In Fig. 16.1b . the shaded circle indicates the 
PV NAtt d-wave vertex. All remaining A^Avr vertices in Fig. 16.11 involve 
strong, PC couplings. In general, the interaction vertices of Fig. 16. II contain 
loop effects as well as tree-level contributions. The loops relevant to the 
PV interactions (up to the chiral order of our analysis) are shown in Figs. 
16.2116.51 The formalism for treating the contributions to A^jj from Figs. 
I6.1b -c is discussed in detail in ref. [83]. Here, we review only those elements 
most germane to the discussion of electroweak radiative corrections. We also 
discuss general features of the new contributions from Figs. 16. Id and 16.1b 
not previously analyzed. 
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Figure 6.2: Meson- nucleoli intermediate state contributions to the 
A transition anapole and Siegert couplings a a and d/\, respectively. The 
shaded circle denotes the PV vertex. The single solid, double solid, dashed, 
and curly lines correspond to the N, A, vr, and 7, respectively. 

6.2.1 Kinematics and PV Asymmetry 

We define the appropriate kinematic variables for the reaction 

e~ {k) + N{p)^ e~'{k') + A (pa) ^ e^' {k') + N' {p') + vr (p^) , (6.7) 
In the laboratory frame one has 

s = (A; + p)^, q = pA-p = k-k', PA=p'+Pn, (6.8) 
where p = 0, and 

s = k^ + 2k-p + p^ = m^ + 2Me + M^, (6.9) 

e being the incoming electron energy, m and M = niN being the electron 
and nucleon masses, respectively. One may relate the square of the four 
momentum transfer 

Q2 = |gf _ ql (6.10) 
to s and the electron scattering angle 9 as 

sin2 9/2 = J, ^ , -. (6.11) 

' {s- M^){s- Ml-Q^) ^ ' 
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Figure 6.3: Same is Fig. 16.21 but with A-vr intermediate states. 

The energy available in the nucleon-gauge boson (7 or Z") center of mass 
(CM) frame is W = \Jl^ and the energy of the gauge boson in the CM 
frame is 

?o = ^ . (6.12) 



As shown in ref. [SS], one may distinguish three separate dynamical 
contributions to the PV asymmetry. Denoting these terms by A^^^ (i = 
1, . . . , 3), one has 



iv+ - iv„ -G^ g2 



N+ + 47ra 



(6.13) 



where (N-) is the number of detected, scattered electrons for an incident 
beam of positive (negative) helicity electrons, a is the electromagnetic fine 
structure constant, and is the Fermi constant measured in /i-decay. The 
'^(1 2) contain the vector current response of the target, arising from the 
interference of the amplitudes in Figs. I6.1h .b. while the term A^g^contains 
the axial vector response function, generated by the interference of Figs. 
16. Ih and l6.1b -e. 

The leading term, A^^^^, is nominally independent of the hadronic struc- 
ture - due to cancellations between the numerator and denominator of the 
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Figure 6.4: Same as Fig. I6.2l but involving insertions of the baryon magnetic 
moment operator, denoted by the cross. 



asymmetry - whereas AJ2 3) are sensitive to details of the hadronic transition 
amplitudes. Specifically, one has 

Afi)=5iCr' , (6.14) 

which includes the entire resonant hadronic vector current contribution to 
the asymmetry. Here, g'^ is the axial vector electron coupling to the and 
^y^^ is the isovector hadron-Z'' vector current coupling [2l [93]: 

5^ey=' = -2(Ci„ - Cm) (6.15) 

where the Ciq are the standard A{e) x V{q) couplings in the effective four 
fermion low-energy Lagrangian [S3]. At tree level, g'^Cy-^^ = 2(1— 2sin^ 9^^) ~ 
1. Vector current conservation and the approximate isospin symmetry of the 
light baryon spectrum protects A^^^^ from receiving large and theoretically 
uncertain QCD corrections. In principle, then, isolation of A^^^^ could pro- 
vide a test of fundamental electroweak couplings. As shown in ref. [83], how- 
ever, theoretical uncertainties associated with the non-resonant background 
contribution A^g) axial vector contribution A^^^ would likely render such 
a program unfeasible. The interest for the Jefferson Lab measurement [37] 
lies in the form factor content of the axial vector contribution A^g^ . For our 
purposes, it is useful to distinguish between the various contributions to this 
response according to the amplitudes of Fig. 16.11 From the interference of 
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Figure 6.5: Same as Fig. 2 but with PV electromagnetic insertions, denoted 
by the overlapping cross and shaded circle. 



Figs. 16. Ih and 16.1b we obtain the axial vector neutral current response: 

Al^^{NC)^g'y^l='F{Q^s) , (6.16) 

where 



-2{C2u - C2d) (6.17) 

in the absence of target-dependent, QCD contributions to the one-quark 
electroweak radiative corrections. The C2q are the V{e) x A{q) analogues of 
the Ciq [M], while the function F{Q'^,s) gives the dependence of A^3^(NC) 

on the axial couplings Cf^. Following ref. [83] we obtain 



1 + 



2M2 



+ 



qo + W-MC^^ 



2M 



(6.18) 



MMa [[s - M2) + {s- Ml) - Q2) 



F{Q\s] 

where 
V{Q\s) 



(6.19) 

In arriving at Eqs. (|6. 16116. 19]) we have included only resonant contributions 
from the A. Non-resonant background effects have been analyzed in refs. 
[831 [95]. Note that F{Q'^,s) is a frame-dependent quantity, depending as 



^ (g2 + (Ma + Mf) (q2 + (Ma - Mf) + (s - M^) {s - Ml) - Q^s 
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it does on q^. However, for simplicity of notation, we have suppressed the 
g'^-dependence in the hst of the arguments. 

The interference of Figs. 16. lb and 16. Id generates the transition anapole 
and Siegert contributions associated with the interactions of Eqs. (16.5|6.6p : 

A[3)(Siegert) + A[3) (anapole) , (6.20) 

while the interference of Figs. 16.1b and 16.1b generates the response associ- 
ated with the PV A^Avr d-wave interaction: 

A[3)(d-wave) . (6.21) 

From the total contribution 

A[3)(T0T) = A[3)(NC) + A[3) (Siegert) + A[3) (anapole) + A[3) (d-wave) 

(6.22) 

we may define the overall 0{a) correction R'^ to the 0{Gf) axial response 
via 

A^3)(T0T) = 2(l-4sin2 0O,)(l + ^A)^(Q^^^) (6-23) 
where 6^ is the weak mixing angle at tree-level in the Standard Model: 

sin2 ^0(1 - sin2 00 ) = —J^ , (6.24) 

or 

sin^ 6'[i[, = 0.21215 ± 0.00002 . (6.25) 

One may decompose the 0{a) effects described by according to several 
sources: 

i?A ^ ^ewk ^ ^Siegert ^ ^anapole ^ ^d-wave ^ . . . ^ (g 26) 

where the + • • • indicate possible contributions from other many-quark and 
QCD effects not included here. The quantity R^^ denotes the one-quark 
radiative corrections, 

'-"2« ~ '-'2d 

with the superscript "0" denoting the tree- level values of the C2q- The 
correction R^^ denotes both the effects of 0{a) corrections to the relation 
in Eq. (I6.24P as well as the 0{aGF) contributions to the neutral current 
e-q amplitude. While the tree-level weak mixing angle is renormalization 
scheme-independent, both sin^ 9w and the correction R^^ depend on the 
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choice of renormalization scheme. In what follows, we quote results for 
both the on-shell renormalization (OSR) and MS schemes. Note that our 

(k) 

convention for the R\ differs from the convention adopted in our earlier 
work of ref. [38], where we normalized to the scheme-dependent quantity 
1 - 4sin2 9w 

The remaining corrections are defined by 

^Siegert ^ A(^3)(Siegert)/A^3)(NC)0 (6.28) 

^anapole ^ A|^3)(anapole)/A^3)(NC)0 (6.29) 

^d-wave ^ A[3)(d-wave)/A^3)(NC)° , (6.30) 

where the "0" denotes the value of the NC contribution at tree-level. 



6.3 Electroweak radiative corrections 

The parity violating amplitude for the process ep — > eA is generated by the 
diagrams in Fig. I6.1b -e. At tree-level in the Standard Model, one has 

iM^"" = iM^^ + iMZ , (6.31) 

where 



iM^^ = i-^l^^ < A\Jx\N > (6.32) 
2v 2 



from Fig. 16.1b and 



iMZ = i^l^ < A\Jx5\N > . (6.33) 

from Fig. 16.1b . Here, Jx (Jas) and Ix (/as) denote the vector (axial vector) 
weak neutral currents of the quarks and electron, respectively [2]. Note that 
the vector leptonic weak neutral current contains the factor gy = (— 1 + 
4sin^0^) ~ —0.1, which strongly suppresses the leading order Z'^-exchange 
amplitude of Fig. 16.1b . 

The interactions given in Eqs. ()6.5|6.6p generate additional contributions 
to when a photon is exchanged between the nucleon and the electron 

(Fig. 16. Id ). The corresponding amplitudes are 

^^Siegert = -i^-^^^n,eAA{M - M^)g^'^ - q'^^^] N (QM) 
^Knapole = ^^-^n^eA,N . (6.35) 
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We note that, unhke M^J, the amplitudes in Eqs. (|6.34p and (j6.35p contain 
no (1 — 4sin^ 9^) suppression. Consequently, the relative importance of the 
PV 7-exchange many-quark amplitudes is enhanced by 1/(1— 4sin^ Oy^-) ~ 10 
relative to the leading order neutral current amplitude. 

The constants dj\ and contain contributions from loops (L) gener- 
ated by the PV heavy baryon chiral Lagrangians (see ref. |36j ) and from 
counterterms (CT) in the tree-level effective Lagrangian of Eqs. (j6.5l6.6p : 

dA = di + dl^ (6.36) 



a A + a A 



(6.37) 



In HBxPT, only the parts of the loop amplitudes non-analytic in quark 
masses niq can be unambigously indentified with and a^. Contributions 
analytic in the uiq have the same form as operators appearing the effective 
chiral Lagrangian, and since the latter carry a priori unknown coefficients 
which must be fit to experimental data, one has no way to distinguish their 
effects from loop contributions analytic in niq. Consequently, all remaining 
analytic terms may be incorporated into d^ and a^^. In principle, d^ 
and should be determined from experiment. At present, however, we 
know of no independent determination of these constants to use as input in 
predicting R^, so we rely on model estimates for this purpose. 

The structure arising from the PV d-wave amplitude (Fig. 16.1b ) is con- 
siderably more complex than those associated with Figs. I6.1b -d. and it will 
not be presented here (see ref. [36] for more details). We mention, however, 
that the amplitude of Fig. 16.1b - like its partners in Fig. 16. Id - does 
not contain the 1 — 4sin^ 9^^ suppression factor associated with the 0{Gp) 
amplitude of Fig. 16.1b . 

For future reference, it is useful to give expressions for the various con- 
tributions to A^3^ as well as the corresponding contributions to and the 
total asymmetry Ai^j^. For the response function, we have 

A^3)(Siegert) 
A^g-j (anapole) 



8\/27ra c^a 
8\/27rQ ca 



A^3j (d-wave) 



8\/27ra 



qo + W-Mn 



V{Q\s) 



(6.38) 



(6.39) 



A. 



Ma + Mn 



g-TzNA 



H{Q\s)V {Q\s) 



(6.40) 

The appearance of 'P[Q'^,s) results from the different kinematic depen- 
dences generated by the transverse PC and axial vector PV contributions to 
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the electroexcitation asymmetry [831 [2]. The function H{Q'^,s) is a gently 
varying function of , and numerical calculation over the kinematics of the 
Jefferson Lab maesurement shows that 



\H{Q^,s)\ < 0.1. (6.41) 
The corresponding radiative corrections are 

Siegert _ SV^vra 1 g^ + W-M . , 

" G,Al 1 - 4sin2 e 2C^ 2A^ ^ ' ^^'^^^ 

^anapole SV^-rra 1 ,,^2 ao\ 

= - G,Al l-4sin^e W^^^''^ ^ ^ ^ 

G^A^ 1 - 4sm^ 9^ niA + ^^ata 

(6.44) 



where 



' ^ 2M2 2M ^ ^ 

In order to set the overall scale of i?^^^^^^^, and R'^'^^^^, we 

follow ref. [38] and set c^a ~ ~ /atAtt ~ dw, where 5^ = 3.8 x 10~^ is 
the "natural" scale for charged current hadronic PV effects |96l [97] . Using 
~ 1, C^/C^ ~ 1.6, g^NA ~ 1, fiQ^s) ~ 1 and H{Q\s) ~ 0.1, we 
obtain 



^Siegert _ q ^^^^ (A^/Q^) (5.46) 
^anapole _ _0.0041 (6.47) 
^d-wave ^ _o.0002 . (6.48) 

As we show below, ^^'^^P*-'^^ may be significantly enhanced over this general 
scale. 

Finally, the total contribution to the asymmetry from the various re- 
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sponse functions is given by 



A,«[A^3)(NC)] 
^^«[Ai3)(Siegert)] 



4\/27ra 



2(Ci„ — Cid) 



9 X 10-5[gV(GeV/c)2] 

2(C2«-C2d)F(Q2,s) 



4\/27ra 

-6.3 X 10-''F(Q^s)[gV(GeV/c)2] 



2dA -5 

-2 X 10^ 



dA/S'Ti 



(6.49) 



(6.50) 



(6.51) 



Aifi[A[3)(d-wave)] 



2.8 X 10 



CI 



QttNA 



H{Q\s)V{Q',s) 

.-SifNAw/gn-, 



V{Q\s)[Q'/{GeY/cf] 
2Q2 



3.52) 



A^(?TiA + mj^) 



3.0 X 10-«[^^^^i^^]//(g2,s)P(Q2,s)[QV(GeV/c) 

At this point, using the chiral Lagrangian describing the PV 7A^A tran- 
sition [36] , one could compute the contributions to oa and generated by 
the loops of Figs. (j6. 2116. 5|) . Loop corrections to the PV vrA^A d-wave inter- 
action contribute at higher order than considered here, so we do not discuss 
them exphcitly. The final result we get by summing up all the possible 
contributions from the different diagrams is 



aiiTOT) 



+ 



67r 
V3 



ISvr 



g-KNAhn 
g-KNAhA 
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m. 
A. 

4 



24 niN 
24 ruN 



Go 
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\/6 '± 

+ ^gnrNA{hv + |^y)Go 

1 uA+A° 
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1 / , pA+7r^ 



)Go 



(6.53) 
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Coupling 


Best value 


Reasonable range 


\d'^''{RES)\ 
al^iVMD) 

InAtt 




IOO5. 
(-15 ^ 70)<7. 
^ I65. 



Table 6.1: Best values and reasonable ranges for d'^, . 
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m. 



(6.54) 



As we mention previously, the constants and d^ consist of two pieces, the 
loop contributions a£ and and the counter term pieces and d^ , 
which are unknown a priori and should be determined from experiments. 
Also the low energy constant /ata. is unknown and should be extracted 
from data. If one uses naive dimensional analysis (NDA), it can be shown 
that all three quantities, , d'^ and Jna-k, are 0(5.). Using vector 
meson dominance for af^ and resonance saturation method for d'^ and 
JnAtt, one could find a reasonable range for this quantities which might 
significantly deviate from the value obtained by using NDA [36] . The results 
are summarized in table 16.11 



6.4 The scale of radiative corrections 



In the absence of target-dependent QCD effects, the 0{aGF) contributions 



to A^"^^ are determined entirely by the one-quark corrections R^^^ as defined 
in Eq. ()6.27p . As noted above, incorporates the effects of both the 

0{a) corrections to the definition of the weak mixing angle in Eq. (j6.24p as 
well as the OlaGp) contributions to the elementary e-q neutral current am- 
plitudes. The precise value of R^^ is renormalization scheme-dependent, 
due to the truncation of the perturbation series at 0{aGF). In Table [621 
we give the values of sin^ 6^^, —2{C2u — C2d), and R^^ in the OSR and MS 
schemes. We note that the impact of the 0{a) one-quark corrections to the 
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tree- level amplitude is already significant, decreasing its value by ~ 50%. 
As noted in Sec. 16.11 this sizable suppression results from the absence in 
various loops of the 1 — 4sin^ ^^v- factor appearing at tree-level, the appear- 
ance of large logarithms of the type In m^/M^, and the shift in sin^ 9w from 
its tree-level valu^. 

In discussing the impact of many-quark corrections, it is useful to con- 
sider a number of perspectives. First, we compare the relative importance 
of the one- and many-quark corrections by studying the ratios Ra ■ Using 
the results previously obtained, we derive numerical expressions for these 
ratios in terms of the various low-energy constants. For the relevant input 
parameters we use current amplitude 1 — 4sin^ O^L, qa = 1.267 it 0.004 



g-KNli = 1.05 



1.166 X 10"^ GeV-^ 6 = 0.3 GeV, ^ = 



1.16 
1.39 



GeV, fp = 5.26, ^ = 17 [lOl], = 3.8 x 10-^ = 0.87 and 
[85j . It is worthwhile mentioning that 2C^ is normalized such that this 
factor becomes qa for polarized ep scattering. We find then 

^anapole ^ ^ ^ 1^ ^ {-QMK - 0.07/iy + 0.006/ia - O.18/i0 

+G.nh1^^^ + ^m\h%p + h\Np - /i2jvpl + 0.025|/iU^|} 

(6.55) 



^Siegert ^ q.OI x 



2Ci 



idM(fJ - 0.09/i^ - 0.03/iA 



O.lGeV^ qo + W -M 
^ 0.6GeV 
^d-wave ^ 0.00105 x/;vA.x(C3VC5^)xi/(Q2,s) 



(6.56) 
(6.57) 



where 



hv 



h 



A+AO 



V 



^/3 

pA+TT+TT" 



T pA + TT 7r^ 



(6.58) 

(6.59) 
(6.60) 



and where all PV couplings are in unit of g-,^ and \H{Cp, s)\ < 0.1. The ex- 
pressions in Eqs. (j6.55p illustrate the sensitivity of the radiative corrections 
to the various PV hadronic couplings. As expected on general grounds, the 



^At this order, the scheme-dependence introduces a 10% variation in the amphtude, 
owing to the omission of higher-order (two-loop and beyond) effects. 
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Scheme 


sin^ 6',y 


— 2(C2u — C2d) 


pewk 


Tree Level 
OSR 
MS 


0.21215 ±0.00002 
0.22288 ± 0.00034 
0.23117 ±0.00016 


0.3028 
0.1404 
0.1246 




-0.536 
-0.589 



Table 6.2: Weak mixing angle and one-quark 0{aGF) contributions to 
isovector axial transition current. 

u\ 

overall size of the R\ is about one percent when the PV couplings assume 
their "natural" scale (NDA). The relative importance of the Siegert's term 
correction, however, grows rapidly when falls below ~ 0.1 (GeV/c)^. The 
hadron resonance models we have used to estimate the size of the counter 

(i) 

terms, may yield significant enhancements of the R\ beyond the naive di- 
mensional analysis scale. To obtain a range of values for the corrections, 
we list in Table [6^ the available theoretical estimates for the PV constants, 
including both the estimates given above as well as those appearing in refs. 
|102l [96] . We observe that the couplings hy, d/\ and /i^Afp ^'^^ weighted 
heavily in the expressions of Eqs. ()6.55p . At present, these couplings are 
unconstrained by conventional analyses of hadronic PV and there exist no 
model estimates for /i^ and hy. Consequently, we allow the various combi- 
nations of these constants appearing in Eq. (16.551) to range between I^Qt^ 
and — lO^TT, using g^^ as a reasonable guess for their best values. 

The resulting values for the i?^*^ are shown in Table 16.41 and Fig. 16.61 

For the ratio , we quote results for two overall signs (±) for c^a, 

since at present the overall phase is uncertain. From both Table 16.41 and 
Fig. 16.61 we observe that the importance of the many-quark corrections 
can be significant in comparison to the one-quark effects R^^. Moreover, 
the theoretical uncertainty, resulting from the reasonable ranges for the PV 
parameters in Table 16.31 can be as large as R^^ itself. It is conceivable 
that the total correction R^ could be as much as ±1 near the lower end of 
the kinematic range for the Jefferson Lab N ^ IS. measurement. While this 
result may seem surprising at first glance, one should keep in mind that the 
0{aGF) one-quark effects already yield a 50% reduction in the tree-level 
axial amplitude, while the absence of the leading factor of in the Siegert 
contribution to enhances the effect of the unknown constant d/\ for 
low momentum transfer. If the Siegert operator is enhanced by the same 
mechanism proposed to account for the violation of Hara's [99] theorem in 
A.S = 1 hyperon radiative decays, then the magnitude of the effects shown 
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Figure 6.6: Contributions to the electroweak radiative correction R'^ at 
beam energy 0.424 GeV. The short-dashed hues show the upper and lower 
bounds of the reasonable range for anapole contribution. The solid line is the 
one-quark contribution. The upper (lower) long-dashed line is the Siegert 
term with d/\ = 25^7^ (— 25g7r)- The dotted line is the d-wave contribution. 
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Coupling constants 


Source 


Best values 


Ran 


ge 


K 


[102| (|103|) 


7(7) 


^ 


17 




m im\) 


-20 (-20) 


-51 - 







m im) 


11 (10) 


5 ^ 


17 


'''ANp 


m im\) 


20 (30) 


-54^ 


152 




m im) 


20 (20) 


17^ 


26 




m im) 


0(0) 


-0.5 - 


2 


hy 


133 


1 


-10 - 


^ 10 




this work 


1 


-10 - 


^ 10 




lioij 


1 


-10 - 


^ 10 



Table 6.3: Range and the best values for the available PV coupling constants 
(in units of g^) from refs. [1021 [Ml [Ml [M] and this work. 

in Table 16.41 and Fig. 16.61 is not unreasonable. Conversely, should a future 
measurement imply ~ R^^, then one may have reason to question the 
resonance saturation model for both c^a and the hyperon decays. 

For the purpose of analyzing prospective measurements, it is also use- 
ful to consider the contributions to the total asymmetry generated by the 
various 0{aGF) effects. In Figs. (I6.7l6.8p . we plot the ratios 

where 74^jj(NC-tot) is the total neutral current contribution to the asym- 
metry and i denotes the Siegert, anapole, and d-wave contributions. In Fig. 
16. 7[ we show the band generated by the anapole, where the limits are de- 
termined by the ranges in Table 16. 4i For simplicity, we show the Siegert 
contribution for only the single case: = 2bgTj, where the effective cou- 
pling c^A contains both the counterterm and loop effects, noting that is 
dominated by . In Fig. 16.81 we give the variation of the Siegert contri- 
bution for a range of values, where this range is essentially determined 
by the range for d^ given in Table 16. 1[ 

From the plots in Figs. ()6.7|6.8p . we observe that the uncertainty as- 
sociated with the anapole and d-wave terms can be as much as ~ 25% of 
the nominal axial NC contribution. The uncertainty associated with the 
Siegert contribution is even more pronounced. For <0.1 (GeV/c)^, this 
uncertainty is ±100% of the axial NC contribution, decreasing to < 15% 
at = 0.5 (GeV/c)^. Evidently, in order to perform a meaningful deter- 
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Figure 6.7: Ratio of asymmetry components rj = A'^Lji/A^^tof' where A^j^^^^ 
denotes the total neutral current contribution. The dotted line gives the 
Siegert contribution; the long-dashed line is for the PV d-wave; the short 
dashed lines give our "reasonable range" for the anapole effect; and the solid 
line is for axial vector neutral current contribution. All the other parameters 
are the same as in Figure 
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Figure 6.8: Same as Fig. 16.71 but omitting the anapole and PV d-wave 
curves and showing Siegert contribution for several values of the coupling 
c^A- The dotted, dashed and dashed-dotted lines are for = Ig^^, 25g7r 
and 100(77r respectively. The solid line is for the axial vector neutral current 
contribution. All the other parameters are the same as in Figure [6771 
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Source 


(best) 


i?^ (range) 


One-quark (SM) 


-0.54 




Siegert (+) 


0.21 


0.02 ^ 0.85 


Siegert (— ) 


-0.21 


-0.85 ^ -0.02 


Anapole 


0.04 


-0.09 ^ 0.21 


d-wave 


0.0006 


-0.003 ^ 0.003 


Total (+) 


-0.29 


-0.61 ^ 0.52 


Total (-) 


-0.71 


-1.48 ^ -0.35 



Table 6.4: One-quark Standard Model (SM) and many-quark anapole and 
Siegert's contributions to V{A) x A{N) radiative corrections. Values are 
computed in the on-shell scheme using = 0.1 (GeV/c)^ and qq + W—M = 
0.6 GeV. The plus and minus signs correspond to the positive and negative 
values for d'^. 



Experimental Parameter 


Benchmark Value 


Luminosity C 
Running time T 
Solid angle AO, 
Electron polarization Pg 


2 X 10^^ cm-^s"^ 
1000 hours 
20 msr 
100% 



Table 6.5: Possible experimental conditions for ^4^^ measurement. 

mination of the Cf'{Q^), one must also determine the size of the Siegert 
contribution. Since the variation of the latter can be as large as that 
associated with the Cf{Q'^) for 0.1<Q^<0.5 (GeV/c)^, one may not be 
able to rely solely on the Q^-dependence of the asymmetry in this regime in 
order to disentangle the various effects. 

Rather, in order to separate the Siegert contribution from the other ax- 
ial terms, one would ideally measure in a regime where the Siegert 
term dominates the asymmetry. As shown in Fig. 16.91 the Siegert contri- 
bution can become as large as the leading, A^-j^^ contribution for < 0.05 

(GeV/c)^. To estimate the experimental kinematics optimal for a determi- 
nation of cZa in this regime, we plot in Fig. 16.101 the total asymmetry for 
low-Q^. To set the scale, we use the benchmark feasibility estimates of ref. 
[83], based on the experimental conditions in Table [631 

From the figure of merit computed in ref. [83], one obtains a prospective 
statistical accuracy of ~ 27% at ^ = 400 MeV, 6 = 180° and = o.054 
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Figure 6.9: Asymmetry components as a function of and beam energy 
0.424 GeV. Except for d^, all the parameters are taken from the central 
values of the table ()6.3p . The bold long-dashed (dashed) line is for Alr{AJ^-^) 
(^Li?('^(2)))' "^^^ solid, dashed-dotted, dotted and dashed lines are for 
Alr{AJ„.) at dA = 0,25,75 and 100 
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(GeV/c)^. A measurement with such precision would barely resolve the 
effect of d/^ = itlOO^Tr. Doubhng the beam energy and going to more forward 
angles {e.g. 6 = 20°), while keeping essentially the same, would reduce 
the statistical uncertainty to roughly 5% . At this level, one would be able 
to resolve the effect of having roughly the size of our "best value" . More 
generally, a forward angle {6 < 20°) measurement for ~ 1 GeV appears to 
offer the most promising possibility for determining d^. Such a measurement 
would have two benefits: (a) providing a test in the AS" = channel of 
the mechanism proposed to explain the violation of Hara's theorem in the 
AS* = 1 hyperon radiative decays, and (b) helping constrain the ^A-related 
uncertainty in an extraction of the Cf{Q'^) for > 0.1 (GeV/c)^. 

6.5 Conclusions 

Parity-violation in the weak interaction has become an important tool for 
probing novel aspects of hadron and nuclear structure. At present, an ex- 
tensive program of of PV electron scattering experiments to determine the 
strange-quark vector form factors of the nucleon is underway at MIT-Bates, 
Jefferson Lab, and Mainz [T05] . A measurement of the neutron radius of 
208p|j planned for the future at Jefferson Lab |106j , and measurements of 
non-leptonic PV observables are being developed at Los Alamos, NIST, and 
Jefferson Lab jl07| . In the present study, we have discussed the application 
of PV electron scattering to study the ^ A transition, which holds sig- 
nificant interest for our understanding of the low-lying qqq spectrum. We 
have argued that: 

(i) The 0(aGi?) contributions to the axial vector N ^ A response gener- 
ate a significant contribution to the PV asymmetry. One must, there- 
fore, take these effects into consideration when interpreting any mea- 
surement of the asymmetry. 

(ii) A substantial fraction of the 0{aGF) contributions arise from weak 
interactions among quarks. A particularly interesting "many-quark" 
contribution of this nature involves the PV ^NA electric dipole cou- 
pling, c^A) whose presence leads to a non- vanishing asymmetry at the 
photon point. 

(iii) A determination of c^a via, e.g., a low-Q^ asymmetry measurement, 
would both sharpen the interpretation of a planned Jefferson Lab PV 
A electroexcitation experiment and shed light on the dynamics of 
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Figure 6.10: Total asymmetry at small \q^\ for several d/^. The couplings 
are at central values of table (|6.3p . The lines for c^a = 0, 25, 75 and 100 g-,^ 
are the solid, dashed, dashed-dotted, dotted and long-dashed line. 
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mesonic and radiative hyperon weak decays. Indeed, one may conceiv- 
ably discover whether the anomalously large violation of QCD sym- 
metries observed in the latter are simply a peculiarity of the AS* = 1 
channel or a more general feature of low-energy hadronic weak inter- 
actions. At the same time, knowledge of would allow one to place 
new constraints on the axial transition form factors Cf{Q'^) from PV 
asymmetry measurements taken over a modest kinematic range. 

(iv) Experimental results for the AS" = 1 decays suggest that the PV 
— > A asymmetry generated by c^a could be large, approaching a 
few xlO~^ as — > 0. Measurement of an asymmetry having this 
magnitude using forward angle kinematics at existing medium energy 
facilities appears to lie within the realm of feasibility. 

More generally, the subject of hadronic effects in electroweak radiative 
corrections has taken on added interest recently in light of new measurements 
of the muon anomalous magnetic moment |108j and backward angle PV elas- 
tic ep and quasielastic ed scattering [86]. The results in both cases differ 
from Standard Model predictions, with implications resting on the degree to 
which one can compute hadronic contributions to radiative processes. The 
interpretation of future precision measurements, including determination of 
the asymmetry parameter in neutron /?-decay and the rate for neutrinoless 
/?/3-decay, will demand a similar degree of confidence in theoretical calcula- 
tions of higher-order, hadronic electroweak effects. Thus, any insight which 
one might derive from studies in other contexts would represent a welcome 
contribution. To this end, a comparison of PV electroexcitation of the A 
with the corresponding neutral current z^-induced A-excitation would be 
particularly interesting, as the latter process is free from the large 0{aGF) 
hadronic effects entering PV electroexcitation [39^ [2] . 
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Appendix A 



Fitting Parameters of the 
Four Fermion Twist Four 
Contributions to the SFs 



In this appendix we present the results for the different parametrizations 
used to fit the moments of the SFs. The results for to contibution to 
Fi^{xb) arising from the operator Ui are: 



= 36.4x10"'' 5 = -1.153 c = 2.797 d = 3.478 e = 6.118 

(A.l) 

and 

= 754.472 x 10"^ 6 = .038 c = 4.317 d = 2.815 e = 2.181 

(A.2) 

for a parametrization of type (j4.8ip and 

ai = -9.00 X 10"^ 02 = 100.08 x 10"^ 03 = -164.16 x lO"'' 

(A.3) 

02 = -7.92 X 10"^ 03 = 231.84 x 10"^ 04 = -432.00 x 10"*^ 

(A.4) 

in the case of a parametrization of tvpe l4.82l Notice that also in this case 
we found two sets of parameters {a, h, c, d, e} which equally fit the moments 
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Figure A.l: Comparison of the different fits for the moment function M{N) 
relative to the Ff^^{xB) SF (C/i contribution). The sohd hne is relative to 
the set of parameters in eq. (jA.ip . the long dashed line is relative to the 
set of parameters in eq. ()A.2p . the dot-dashed line is relative to the set of 
parameters in eq. ()A.3p . and the short dashed line is relative to the set of 
parameters in eq. ()A.4p . The three moments obtained from the bag model 
are represented by stars. 



of the SF. In Fig lA.ll we have plotted the momentum function M{N) and 
in Fig lA.2l we have plotted the corresponding contribution to xeF^^^xb)- 
The contribution of the operator Ui to Ff^^^'^{xB) gave as a result the 
following set of parameters 



-9.631x10""* 6 = -1.153 c = 2.797 d = 3.478 e = 6.118 



(A.5) 



a = -199.621 X 10-4 6 = .038 c = 4.317 (i = 2.815 e = 2.181 



ai = 2.381 X 10-4 02 = -26.480 x 10"^ 03 = 43.434 x IQ-^ 



(A.6) 



(A.7) 



02 = 2.096 X 10""* as = -61.341 x IQ-^ 04 = 114.300 x 10-^ 
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Figure A. 2: Reconstruction of the SF xbFi^\xb) through the use of the 
inverse MeUin transform technique {Ui contribution). The type of hue used 
for each parametrization is the same as in Fig lA.li 



The graphs of M{N) and xbF^^^'^{xb) are in Figs lA.3l and IA.4I respec- 
tively. The parameters for the contribution of the operator U2 to the SF 
Ff ^^^(xb) are 

a = 302.360 x 10"^ b = 1.545 c = 9.362 d = 265.331 e = 3.006 

(A.9) 

and 

= 77.404 x 10-"^ 6 = -.112 c = 7.477 d = 5.447 e = 179.207 

(A.IO) 



ai = -8.709 X lO^'' as = 48.282 x 10"^ 03 = 36.746 x lO"'' 



(A.ll) 



02 = -59.569 X 10"^ 03 = 433.076 x 10"^ 04 = -431.403 x 10"^ 

(A.12) 

The graphs relative to M{N) and xbFi^^'~^{xb) are plot in Figs I A. 5 1 and 
IA.6i Finally the parameters relative to the contribution of the operator Ui 
to the SF F^^{xb) are 

a = -15.522 X 10"^ 6 = -.868 c = 4.517 d = 35.106 e = 45.197 

(A.13) 
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Figure A. 3: Comparison of the different fits for the moment function M{N) 
relative to the F^^^^^ {xb) SF {Ui contribution). The sohd hne is relative 
to the set of parameters in eq. (lA.5p . the long dashed line is relative to the 
set of parameters in eq. ()A.6p . the dot-dashed line is relative to the set of 
parameters in eq. (lA.Th . and the short dashed line is relative to the set of 
parameters in eq. (jA.Sp . The three moments obtained from the bag model 
are represented by stars. 



ai = 17.213 X 10" 



a2 



-132.571 X 10"^ ag = 115.358 x 10" 



(A.14) 



02 = 22.344 X 10"^ as = -280.536 x 10"^ 04 = 258.191 x 10"^ 

(A.15) 

and the graphs for M{N) and xbF^'~^ [xb) are plotted in Figs IA.7I and IA.81 
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X 

Figure A. 4: Reconstruction of the SF xbFi^^^'~'{xb) through the use of 
the inverse MeUin transform technique {Ui contribution). The type of hue 
used for each parametrization is the same as in Fig lA.Si 




Figure A. 5: Comparison of the different fits for the moment function M{N) 
relative to the Ff'^^^'-" (xb) SF {U2 contribution). The sohd hne is relative 
to the set of parameters in eq. (|A.9P , the long dashed line is relative to the 
set of parameters in eq. (jA.lOp . the dot-dashed line is relative to the set of 
parameters in eq. (jA.lip . and the short dashed line is relative to the set of 
parameters in eq. (IA.12|) . The three moments obtained from the bag model 
are represented by stars. 
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Figure A. 6: Reconstruction of the SF XfiFf ^^^^(xs) through the use of 
the inverse Melhn transform technique (U2 contribution). The type of hne 
used for each parametrization is the same as in Fig lA.51 




Figure A. 7: Comparison of the different fits for the moment function M{N) 
relative to the F^^^^^ {xb) SF {Ui contribution). The sohd hne is relative 
to the set of parameters in eq. (jA.lSh . the dashed line is relative to the set 
of parameters in eq. (|A.14p and the dot-dashed line is relative to the set of 
parameters in eq. (lA.lSh . The three moments obtained from the bag model 
are represented by stars. 
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Figure A. 8: Reconstruction of the SF xbF^^^^'" {xb) through the use of 
the inverse MeUin transform technique {Ui contribution). The type of hne 
used for each parametrization is the same as in Fig lA.Ti 
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Appendix B 

Spin-Flavor-Color Factors 



In this appendix we are going to present the spin-flavor-color factors for the 
different contributions to SFs. Let's first consider the spin-flavor part. If we 
call Pg = qaQa the projector operator onto flavor q and spin a, by using eq. 
(|iA2]l we have 



iP\P^np~^/'\P) = (p|P7^/vy2|p) = I 

_ 1 

~ 3 



{p\p~'/'py'\p) = {p\py'p^'/'\p) 



iP\P-y^P~y^\P) = {p\p~'/'p-y^\p) 

,1/2 p-l/2| 

p\ _ / PI p 

1/2 pl/2| 

2 

3 '''''' ""^ ' 3 



^p|pi/2p-i/2|p^ = {P\PJ''^P^'^\P) 

{p\py^py^\p) = I {p\p;^^''^p~^'^\p) 



(B.l) 



For the color factors, we have 
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Table B.l: Color factors for the opertor {P\u^ Xauu^ X'^u\P) 



{P\u'f XaUU^ X^'ulP) 



X" 

A2 
A3 
A^ 
A5 
A6 
A7 
A8 



Table B.2: Color factors for the opertor {P\S Xadd^ X^dlP) 



X^ 


(P|(iUaddU"d|P) 


x^ 





A2 





A3 


2 
3 


A^ 





A^ 





A^ 





A^ 





A^ 


2 

3 



Table B.3: Color factors for the opertor {P\d'^ Xadu^ X"'u\P) 



{P\d^ Xadu^ X^ulP) 



X" 

A2 

A3 
A^ 
A^ 
A^ 
A^ 
A^ 
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Appendix C 



Moments Formulae for the 
Four Fermion Contributions 



In this appendix I am going to present the rest of the formulae derived in 
Sec. 14.3.11 to compute the moments of the four fermion operators Ui and 
U2. 



C.l Useful Formulae I 



Proceeding in the same way that lead to Eq. (j4.67|) we find, for the other 
three A functions 



X 



N—1 j-l j-k-1 j — l-m 

EEE E 

j=l k=0 m=0 p=0 



j — 1 — m 
P 



X (-^ ^ ^)(n-9r/i(0)(n-5r/2(0)(n- 5)^-2-^-/3(0) 

(C.l) 
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and 



J dx dy dz j dxBX^A{y — x, z,y,XB) 



dX dfi du 



.g*^a;gj/x(3/-a;)gji/(z- 



27r 27r 2-k 

N—lj-lj-k~l t j—l—m 

E E E E E ^"-^(-D* 

j=l k=0 t=0 m=0 p=0 



^Vl(^)/2(/^)/3(A) 



j-k-1 
t 



^ "'V™-^r/i(0)(^- W2(0)(n-9)^-2-?'— /3(0) 
mj \ p 



(C.2) 



and 



^ ^ *1 iAx iMy-) e-(--J/) /i (z.)/2 (;x)/3 ( A) 
27r27r27r ''^ ^^^f-z^JV j 



j=l fc=0 t=0 m=0 p=0 



(C.3) 



c.2 Useful Formulae II 



In Sec. 14.3.11 we computed the contribution of the operator U2 to the mo- 
ments of the SFs (see Eq. (|4.72p ). The same calculation can be easily 
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extended to the operator Ui by noticing that Eq. (|4.70p would now read 



(n • dfqa{x)ii{n • dYq„{x) 

]\t2 p 

(n • dnjo{E\^\)e'^'){n ■ dr^{jo{E\^\)e-'^') 



M L 



i{n . dr{j,{E\^\)e-''^'){n • dy{h{E\i,\) 



JEt 



X 



+ i{n . 9)^'(jo(i?|x|)e*^*)(n • dr{h{E\^\)fie 



-iEt\ 



X 



X 



+ ix' + y')in-dn 



){n-dY 



X 



(C.4) 



Notice that in this case there is not a spin dependent factor (—1)'^ and notice 
the change of sign in the last factor in the previous relation compared to 
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Eq. (|4.70p . Using Eq. ()C.4p we get the operator Ui 

/<ixBxgA(.,,,,.,«)/d.*d.f/,(.....) 



,^2 N-lj-lj-k-lj-l-m 
jr=l fc=0 m=0 p=0 



N-2 



j — 1 — m\ f j — k — 1 



ix[(io(i?|x|)e^^*)(n.ar(jo(i?|x| 



P 

-iEt\ 



m 



- i{n.dr{jo{E\^\)e-'^'){h{E\^ 



JEt\ 



+ i(jo(i^|x|)e*^*)(n.a)™(ii(i?|x|)f^e 



|X| 

+ (ji(ii;|x|)^e^^*)(n • 5)-(ji(ii;|x|)^e-^*) 

n • 9f (jo(^|x|)e^^*)(^ • a)^(io(i?|x|)e-^^*) 
i(n • dy{jo{E\^\)e-'^'){n • a)P(ji(i?|x|)^e^^*) 



+ i(n • dr{jo{E\^\)e''''){n ■ 9)^(ii(i?|x|) 



+ (n • 5r(ji(^|x|)^e*^*)(n • dnh{E\^\)f^e- 



X 



X 



X 



(C.5) 



The rest of the contributions for the operaors Ui and U2 arising from the 
other A functions can be easily computed by using Eqs. (|C.1IIC.3() and Eqs. 
Km and (lOil) . 
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Appendix D 



Bremsstrahlung Tensors 



In this appendix we presents the results for the leptonic tensors (the elec- 
tromagnetic and the interference) , and also the contractions of the leptonic 
tensors with the hadronic ones. The electromagnetic leptonic tensor has the 
following form 



_ k-k2k-q + ki-k2{m^ + ^■^2) v 1 

kiMk-k2Y ' ' {ki-k2y{k-k2Y 

^^k^'k''({m^ -Q^ + 2k2-q - 2k-k2 - 2k-q){k-k2f 
+ 2ki-k2k-k2{—m? — k2-q + k-k2 + k-q) + (fci • /c2)^m^^ 
+ 4(A:^A;^ + k^k!^)(^{2m^ + k-k2){ki-k2f 

- k-k2ki-k2{2m^ + k2-q - 2k-q) + {k-k2f{k-k2 - k2-q)^ 

- 45^^ [((A;-A:2)2 + {k-k2 + k-q)m^ - k2-q{m^ + fe-fca)) {ki-k2f 
+ k-k2ki-k2(^- k-k2Q'^ + 2{-k2-q + k-k2 + k-q){k-q - m^)) 

+ {k-k2f{{k-k2f + (m^ - 2k-q)k-k2 - k2-q{m'^ + k-k2 - 2k-q) 



+ k-qiwJ^-Q^ 



2k- 



,))]}. 



(D.l) 
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while, for the interference one, we get 



TflV 



+ 



+ 



2g^m(k^s''{ki-k2 - k-k2){ki-k2 + k2-q - k-k2) 



{ki-k2Yk-k2 

- k^s^iiki-kif + {k2-q - 2k-k2)ki-k2 + {k2-q - k-k2)k-k2] 

+ k-k2ki-k2[^{-'n? + k-k2 + k-q) - S/ca-g] + {k-k2f[^k2-q 

- 3k-k2 + 2{m'^ - Q"^ - 2k-q)]J s^'k" + {{k-k2 - 2m^){ki-k2f 

+ k-k2ki-k2{2m^ + k2-q- 2k-q) + {k2-q - k-k2){k-k2f^ s^'k!^ 

- 2g\mg^"' {[-k2-qk2-s + {k-k2 + k-q)k2-s 

+ (m^ - k-k2)q-s]{ki-k2f + k-k2[-k2-siQ^ + k-q) 

+ (-2m^ - 2k2-q + 3k-k2 + 2k-q)q-s]ki-k2 + {k-k2)'^[k2-sk-k2 

- k2-q{k2-s - 2q-s) + (m^ - - 2k-k2 - 2k-q)q-s\ 



4g^m 



ki-k2k2- s)k^k2 + [ki-k2k2-s + k-k2{q-s 



ki-k2{k-k2f 
k2-s)]k>^k'' - {k2-s){ki-k2 - k-k2)kt'k'' 

{ki-k2k2-s + k-k2q-s)k2k2 
2ig^m 



J^jkaqisss 



- [(^1-^2 - k■k2)k^' + 2(A:-fc2)A:^](e°^^^ - ^ i^)k2aq/3Ss 

+ {ki-k2 - k-k2fe''^^^ko,Sp - {{ki-k2f + {k2-q - k-q)ki-k2 

- k-k2{k2-q - Q2 + k-k2 - k-q)^e''^^"'k2asp 



2ig^'m 



k2- s{ki-k2 - k-k2fe''^^''k^c 



{ki-k2f{k-k2f 

- {ki-k2-k-k2){{ki-k2 + 2k-k2)k2-s + k-k2q-sy''^^''k2aqp 
+ {{k-k2 - m^){ki-k2f + k-k2ki-k2{2m^ + k2-q - 3k-k2 

- 2k-q) + {k-k2)^{-m^ + Q^-k2-q + 2k-k2 + 2A:-g)) e^^^^^g^s^ 



(D.2) 
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For both tensors we dropped terms proportional to and q^, because gauge 
invar iance implies that q^W^^ = q^W^'^ = 0. 



The contraction of the electromagnetic leptonic and hadronic tensors is 



' ' fiu em 



+ 



+ 

+ 
+ 
+ 



{ki-k2)Hk-k2 



^ ^8(^[m'^ + 2k-qm^ + {k-k2f 

2k2-q{m^ + k-k2)]{ki-k2f + k-k2[2{k-q - rr?) 

{w? + 2k-q) + k2-q{2w? + k-k2 - ^k-q) 

2k-k2{k-q - Q^)]ki-k2 + {k-k2f[m^ + {k-k2f - Aik-qy 

k2-qk-k2 + Ak2-qk-q - Ak-k2k-q - (m^ + 2k-q)Q^])Wf' 



[ki-k2 - k-k2][{ki-k2 - k-k2)m + 2k-k2{-k2-q 
k-k2 + k-q)] - {k-k2fQ^^E^ - 4ki-k2{ki-k2 

k-k2)kQvn?'E — 2k-k2ko{{ki-k2)'^ — k2-qki-k2 + 2k-qki-k2 
{k-k2f - k2-qk-k2)E + {ki-k2)\2kl - 2k2-q 
k-k2 + k-q)m^ + k-k2{2kQ — k2-q + k-k2)] 
{k-k2f[{k-k2f + {m^ - 2k-q)k-k2 - k2-q{m^ + k-k2 
2k-q) + k-q{m? - 2k-q - Q^)] + ki-k2k-k2[-2{k-k2 



+ k-q)wp- + 2k-q{kl + k-k2 + k-q) - 2k2-q{k-q 



m 



k-k2Q^ 



(D.3) 



while the contraction of the interference leptonic tensor with the hadronic 
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one IS 



Hu nc 



2k2-qk2-s + {k-k2 + 2k-q) 



k2-s + 2{m — k-k2)q- s\{ki-k2) + k-k2[ — ^iQ-sm 



3k2-sk-q — 2k2-sQ^ + k2-q{k2-s — 4g'-s) + 5k-k2q-s 



+ 4k- 



■s^ki-k2 + {k-k2)'^[k2-sk-k2 - k2-q{k2-s - iq-s) 



+ 2{m^ - 2k-k2 - 2k-q - Q'^)q-s 



2ki-k2k2-s{ki-k2 — k-k2)'mE'^ — 2( \^k2- sk^m 



+ \/W^^{k-k2-m^)]{ki-k2f + k-k2[2rr?^/W^^ 



k2-sko'm + komq-s + 2{k2-q — 3k-k2 — 2k-q)\/ E"^ — m?] ki • ^2 
(A;-fc2)V^2 - m2(m2 + 2k2-q - 2k-k2 - 2k-q - Q^))e 



+ ki-k2k-k2( - 2q-sm^ + 2kom'^^/E'^ - m? 



+ k2-s{—k-q — Q^)m + g-sm[2A;g — 2k2-q + 3k-k2 + 2k 

\2 



+ 2ko\/ E"^ — m? [k2 ■ q — k-k2 — k-q]j + (/ci • ^2)^ — k2-qk2-sm 

+ k2-s{2kQ + k-k2 + k-q)m — {k-k2 — m'^){mq-s 

- 2kQ\/ E"^ - m?)^ + {k-k2f(q-sra{m^ -2k-k2-2k-q- Q'^ 

+ k2-sk-k2m + k2- q{2yEP—rn?kQ — k2-sm + 2mq-s) 



2k-k2koVE'^ -m"^ 



+ 85, 



mki-k2k2-s[[E{k-k2 - ki-k2) + {ki-k2 + k-k2)ko]Q'^ 
{E - E')[ki-k2{k2-q - k-q) + k-k2{k2-q + k-q)]J 

{k-k2f{Q^f + {k-k2[{ki-k2 - 2k-k2){ki-k2 + k2-q 



+ VE^ 



- k-k2) - 2{ki-k2 - k-k2)k-q] - {ki-k2 - k-k2Ym']Q' + {ki-k2 
+ k2-q - k-k2)k-k2[ki-k2{k2-q - k-q) + k-k2{k2-q + k-q)]^ 

+ m(^EQ\k-k2f - E'Q\k-k2f + ko[ki-k2Q^ - {ki-k2 4 

- k-k2){ki-k2 + k-k2)]k-k2 + E{ki-k2 - k-k2) 
[{ki-k2 - k-k2)m'^ + k-k2{-k2-q + k-k2 + 2k-q)] 

+ E'{ki-k2 - k-k2){ki-k2{k-k2 - m^) + k-k2[m'^ + 2k2-q 



k2-q 



2{k-k2 + k-q)\})q-s]wr]. 
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(D.4) 



CTEQ 

Bag Model 




CTEQ 

— — Mod. Bag Model 
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the photon sphtting to quark/anti-quark pairs (which then exchange the Pomeron with the 
nucleon) . 

Clearly a better treatment of the electroweak coupling at low is needed. At present 
there is no universally accepted description of the interactions at low (where perturbative 
QCD cannot be readily applied due to values of as((5^) approaching unity), but perhaps 
progress can be made if a complete multipole description is obtained from a phase shift 
analysis of all available electron scattering and photoproduction data. Even then, the 
treatment of the strange quark mass can be problematic. Ideally, measurements of the 
inelastic asymmetry should be made near the E158 kinematics. 

Given these uncertainties, I assume A)LR = (1.0 x 10~^)Q^ at all values of z/, which is 
the value for the A(3, 3) resonance, and about 10% above the average of the DIS, Pomeron, 
and VDM predictions. At present, it is hard to assign an error to this assumption. 

Similar to the case of elastic scattering, a scattering probability Pj is evaluated using 

Pi = {da / dQde')27it{e - Pmin)[cos{6min) - cos{6max)] 

and scattering takes place if a random number R < PiEj, where the enhancement factor Ej 
is by default the same as Em and Ee- The final results do not depend on Ej for Ej < 200. 



0.6 




0.5 1.0 5.0 10.0 50.0 

k (6eV] 

Figure 2: The data points are a compilation of the world data on o"t as measured with real 
photons. The red curve is the simple fit used in the E158 modeling. 
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DISSERTATION PROPOSAL FOR THE Ph.D. DEGREE 

Hadronic Effects in Parity Violating Electron Scattering. 

Gian Franco Sacco 

Abstract: Since their advent, Parity Violating experiments have played a crucial role in 
testing the standard model (SM) of electroweak interactions. In recent years, these experi- 
ments have reached such a high degree of accuracy that they might test for physics beyond 
the SM. When making high precision measurements, a recurrent problem is how to properly 
include hadronic corrections. In this study I consider possible hadronic effects that appear 
in such experiments. 

In the first part, I consider the possible hadronic corrections to the deep inelastic parity 
violating (PV) asymmetry, A^^ [1]? stemming from sea quarks, perturbative QCD (PQCD), 
target mass (TM) and higher twist (HT) corrections. To estimate the first three effects 
I use available parametrizations of the leading twist parton distribution functions (PDF). 
For the HT, which are not well known, I suggest using the MIT bag model to give some 
reasonable estimate. HT corrections could be relevant at the kinematics in the proposed 
experiments [2,3] where the weak mixing angle 9w niay be measured by using Ar^. In this 
part of the study I consider the implications of these corrections for the extraction of sin^ 9^. 

In the second part I study the contribution to Arl of the bremsstrahlung radiation in 
electron-proton scattering. Such a process constitutes a background in the SLAC experi- 
ment E-158. In the experiment, which uses polarized M0ller scattering to measure A^i at 
the 8% level, it was noticed that the inelastic electron-proton (EP) scattering background 
amounted for 40% of the total asymmetry, which implies that, in order to reach the pre- 
scribed precision in the Arl, one has to know the EP background to a 20% level. I suggest 
to use a simple model, consisting in only two resonances, to investigate which fraction of 
the bremmstrahlung radiation in the proton electro-excitation contributes to the total back- 
ground and if the lack of knowledge in the resonance form factors might spoil the prescribed 
accuracy of the experiment. 

In the last part of the study I estimate the contribution to Arl in the electro-excitation 
of the A(1232) resonance in chiral perturbation theory. I specifically focus on the chiral 
corrections of the terms stemming from the electric-dipole, the anapole and the d-wave chi- 
ral effective lagrangian and give some reasonable range of values for their contribution to 
Arl in the low region. In particular the dipole (or Siegert) contribution assumes great 
importance at very low since it does not vanish at the photon point, while all the other 
contributions do. 
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